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PREFACE xili 


Project MAC was organized at the Massachusetts Institute of Technology 
in the spring of 1963 for the purpose of conducting a research and develop- 
ment program on Machine-Aided Cognition ana Multiple-Access Computer 
systems. It operates under contract with the Office of Naval Research, 
acting on behalf of the Advanced Research Projects Agency of the Depart- 


ment of Defense, 


The broad goal of Project MAC is the experimental investigation of new 
ways in which on-line use of computers can aid people in their individual 
intellectual work, whether research, engineering design, management, 

or education. One envisions an intimate collaboration between man and 
computer system in the form of a real-time dialogue where both parties 
contribute their best capabilities. Thus, an essential part of the research 
effort is the evolutionary development of a large, multiple-access computer 
system that is easily and independently accessible to a large number of 
people, and truly responsive to their individual needs. The MAC computer 
system is a first step in this direction and is the result of research initiated 


several years ago at the MIT Computation Center. 


Project MAC was organized in the form of an interdepartmental, inter- 

laboratory "project'' to encourage widespread participation from the MIT 

community. Such widespread participation is essential to the broad, long- 

term project goals for three main reasons: exploring the usefulness of 

on-line use of computers in a variety of fields, providing a realistic com- 

munity of users for evaluating the operation of the MAC computer system, ’ 
and encouraging the development of new programming ‘and other computer 


techniques in an effort to meet specific needs. 


Faculty, research staff, and students from fourteen ac ademic departments 
and four interdepartmental research laboratories are participating in 
Project MAC, For reporting purposes, they are divided into sixteen groups, 
whose names correspond in many cases to those of MIT schools, depart- 
ments and research laboratories. Some of the groups deal with research 


topics that fall under the heading of computer sciences; others with research 
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topics which, while contributing in a substantiv2 way to the goals of Project 
MAC, are primarily motivated by objectives outside the computer field. 
J 

The purpose of this Progress Report is to outline the broad spectrum of 
research being carried out as part of Project MAC, Internal memoranda 
of Project MAC are listed in Appendix A, and MAC-related theses are 
listed a Appendix B, Some of the research is cosponsored by other 
governmental and private agencies, and its results are described in 
journal articles and reports emianating from the various MIT depart- 
ments and laboratories participating in Project MAC. Such publications 
are listed in Appendix C of the report. Project MAC Technical Reports 


are listed in Appendix D. 
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A Visually-Controlled Manipulator - Marvin L. Minsky 


Since the completion of a computer-controlled mechanical hand in 1961 

(H. Ernst, Ph. D. Thesis, M.I.T., 1361) little has been done in the area 

of design of autonomous manipulators. We have been working the develop- : 
ment of a versatile, computer-operated, visually-oriented machine for 

handling objects in complicated spatial situations. A real-time, on-line, 

television-camera interface has been built for the PDP-6 computer and 

has been used to track the motion of simple visual objects, including the 

human eye. This is of interest in studying the mechanism of humai 

visual preception. Also, an AMF Versatran, Model C, mechanical 


remote manipulator has been coupled to the PDP-6, 


Unrecognizable Sets of Numbers - Marvin L. Minsky and Seymour Papert 


The problem here is, given a set A of non-negative integers, is there 

a finite-state device that will discriminate between members and non- 
members of A when the numbers are presented in binary, or some 

other similar number-representation system with a different radix (such 
as ternary of decimal)? A technique for showing that cértain sets are 

not recognizable has been written. The technique deals with the asymptotic 
behavior of the function TA (n), which is the number of integers less than 

n inthe set A; unless T, (n) has a specified behavior, A is not recogniza- 
ble. However, not all sets A with a satisfactory asymptotic behavior of 
Ty (n) are recognizable. Another technique developed by McNaughton 
handles some of these cases, The goal for research in this area is to 
determine a necessary and sufficient condition that A be recognizable 

in r-ary (r 2 2) bya finite state device. A paper will appear soon in the 
Journal of the ACM. 
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4 ARTIFICIAL INTELLIGENCE 


The Algebraic Approach to Finite Automata - Seymour Papert and 
Robert NcNaughton 


The work by Schutzenberger, Rhodes and Krohn, and some of our own 
results have convinced us of the fruitfulness of the algebraic approach 
to finite automata, and, in particular, the concept of the scmi-group of 
the machine, The two problems that have been most impressively dealt 
with are the problem of decomposition of machines, and the problem of 
general star height. It is our intention to continue to focus on this area. 
Research, however, will lean toward application of already existing 
knowledge to problems about machines, rather than advancing knowledge 
in the field of algebra. We have heen concerned with simplifying some 
of the proofs of known theorems (applying algebra to machines) both to 
make fundamental notions clear enough for future research and thinking, 
and also to bring these concepts within the purview and parlance of begin- 


ning students in the field. 


It is felt that use of the computer will be fruitful in this area. Some of So 
the programmers in the Artificial Intelligence Group have been writing 
LISP programs that will do such things as find the semigroup of a given 
finite-state machine. Use of a time-shared computer in this way can be 
very useful in checking out hunches in the search for valid gencralities 


or in providing counter-examples, 


ARTIFICIAL INTELLIGENCE % 


The Design of LISP 2 - Daniel G. Bobrow 


The design of a new programming language, called LISP 2, is now being 
carried out jointly by Systems Development Corporation and Information 
International Incorporated, with aid from members of the Project MAC 
Artificial Intelligence Group. LISP 2 will be an ALGOL-style language and 
will have the desirable ALGOL features of dynamic allocation of storage, 
recursive functions, and dynamic arrays. LISP 2 will have list process- 
ing as in the current LISP 1.5 language; however, the control language 
will be in ALGOL rather than in LISP. In this list processing, useful 

list structures will be compacted to provide room for the allocation of 


space for arrays and for new list-structure elements. 


In addition to the LISP 1.5 function-oriented list processing, LISP 2 will 
embody format-directed string processing similar to that of the COMIT 
programming language. With this processing a user can specify a pattern 
against which a string is matched. If it is matched succes sfully, the ‘ 
string is parsed, that is, subdivided into segments which match.the 
individual patterns in the given pattern. These segments can then be 

used to construct a new segment according toa given format. This 
processing is useful in the manipulation of natural language strings. 

Since the front end of the LISP 2 system will contain a finite-state machine 
as an input-string grouper and a syntax-directed compiler, the user will 
be able to change the language ‘easily, specify new syntactic forms, define 


new data types, and make semantic additions to thé language. 


The compiler for LISP 2, first compiled with a bootstrap compiler written 
in LISP 1.5, will be written almost conmipletely in LISP 2 itself; LISP 2 
will thus be relatively independent of machine types, and we expect to ou, wy 
create versions of it on a number of machines, such as the Q-32 at SUDYC,, 
the IBM 360, the DEC PDP-6, and the GE 645. The compiler will work 
in a two-stage process. First, the external communication language will 
be translated into an intermediate language, which is compatible upward 
from the current LISP 1.5 notation. Then a list-structure assembly- 
language program will be generated and compiled by a simple machine- 
code assembly program. A preliminary version of the LISP 2 system i 


should soon be in operation on the Q-32 at SDC. =, e 
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6 ARTIFICIAL INTELLIGENCE 


Man-Computer Symbiosis - Warren Teitelman 


A LISP programming system is being constructed to facilitate the develop- 
ment of sophisticated problem-solving systems. Emphasis is on easing 
the familiar debugging sequence: write somne code, run the program, 
make some changes, write more code, run the program again, etc. As 

a system becomes more complex, making changes becomes harder and 


harder. 


Our goal-is to make the computer, via a programming system, play an 
active role in this modification process by providing the means whereby 
changes can be effected immediately, in ways that seem natural to the 
user. The flavor of the system is such that the user feels he is giving 
advice, or making suggestions to the computer about the operations of 
his programs; the necessary work to achieve the desired effect is per- 
formed by the symbiotic systein. The system thus acts as an interface 


between the user and his programs. 


A ‘format list-processing language", FLIP, has been developed and 
implemented. FLIP performs operations of the type done by COMIT,, 
but is considerably more sophisticated and operates from within a LISP 
system. This language will be used to process the advice given by the 
user and will enabie the symbiotic system to accept suggestions ina 
very relaxed format. It has already been used to construct a set of 


editing functions for LISP programs. 


Consistent with emphasis on the importance of change and flexibility, the 
system itself may easily be modified by giving it advice about its own 


Operation. 
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Mathematical Assistant - Williani A. Martin 
eae LAE ES SOTO LAI 


I have embedded a picture-display language in the LISP language of the 
PDP-6. Now complex display programs, such as the ARGUS character- 
recognition scheme of W. Teitelman, are easily written. I have also 
developed a LISP program that displays mathematical expressions and 
equations. At present, the displayed equations have substantially the 
same form as those appearing in mathematical publications. This 
mathematical display program in. PDP-6 LISP communicates with the 
rest of the mathematical assistant programs in the time-shared MAC 


LISP through an input/output dataphone system. 


The mathematical assistant system is a new technique for the comparison 
and simplification of symbolic mathematical expressions. It uses finite- 
field arithmetic for hash-coding functions of a complex variable. By con- 
structing a large but finite field (p=8, 589, 949, 373) within which there are 
representations of the integers and, in an appropriate sense, one can 
use exponent arithmetic to map symbolic expressions onto integers. If 
two expressions are equivalent with respect to rational, exponential, and 
some trigonometric identities, they will, with high probability, map onto 
the same integer. In the recognition of common subexpressions, this 
system appears to be faster and more accurate than those which use 


canonical orderings. 
Heuristics of Theorem Proving in Group Theory - Lewis M. Norton 
et LL DLE OGY, 


We are developing a system of heuristics of theorem proving in group 
theory. It is coded in LISP and handles problems that involve the simple 
consequences of definitions. We are now using it to test methods of re- 
ducing irrelevant effort in proofs already produced. Rather than use 
explicit forms of the predicate calculus, it uses forms more similar to 
standard mathematical practice. Thus, the definition of abelian for a 
group G is written: . 

((ABELIAN G G) IMPLIES (AND (MEMBER A1 G) (MEMBER A2 G)) 

(EQUAL(*PROD Al A2 G) (*PROD A2 Al G))). 

Since logical axioms are embedded in the flow of the program, the 


involved theorems can be proven more quickly than has been done by 
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other methods. For example, the system proves: 
Given: G a group, J a subgroupof G, H a normal 
subgroup of G. 
Prove: J 8 His normal in J. 
In this work, we are trying to simulate the analytical and rational 


methods of human problem solving. 
Computer Experiments in Finite Algebra - W. Douglas Maurer 


Any set of mathematical methods and techniques may be mechanized on 
a digital computer. Even infinite mathematical structures may be 
studied, provided that the expressions used to speak about them are 
finite, since it is these symbols which the computer manipulates. The 
first mechanization is for finite algebra. This system for finite algebra 
consists of an English-like language with 80 commands for the manipu- 
lation of finite groups and semigroups, subsets, maps, and labeled 


constants. Rings and fields are to be added later. 


The system accepts an arbitrary semigroup from the console, if the 
subgroup is given by its Cayley table (multiplication table). It can do 


the following: 


1. construct the Cayley tables of four special kinds of 
group and three special kinds of semigroup: 

2. take the direct product of any two Cayley tables and 
produce a third: 

3. venerate the subsemigroup of any semigroup generated 
by a single element or subset of that semigroup and then proceed to 
build the Cayley table of the subsemigroup: 

4. find all the subsemigroups of a semigroup and all the 
normal and maximal subgroups of a group; 

tay 5, find all the left, right, two-sided, or only the maximal 
ideals of a semigroup: 
6. test whether a Cayley table is associative and whether 
a semigroup has an identity, is a group, and is abelian; 
7. tell whether a map from one semigroup is homomorphic 


to, one-to-one, or onto another semigroup; 
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8. print out, erase from its storage area, or rename any 
Cayley table, subset, map, and labeled constant: 

9. adda zero element and a unit element to a semigroup: 

10. produce the natural map of a group onto the set of left, 
or right, cosets modulo a subgroup, or the Cayley table of the factor 


group modulo a subgroup. 


I have recently added a function to the system which determines whether 
any two given Cayley tables are isomorphic and, if so, demonstrates the 
isomorphism. This function runs in less than thirty seconds for two 
groups of twelfth order. I plan to introduce a program which accepts 
mathematical statements in a source language and produces object codes 
to verify or disprove the statements over a restricted class cf test objects, 


such as a collection of counterexamples. 


A Chess-Playing Program - Burton H, Bloom 


An elaborate static-position evaluator has been substantially debugged 
for the CHESS program. It includes material balance, development, 
center-control, and space evaluations. Current plans are to also in- 
clude king-safety and pawn-weakness and some other criteria. Several 
plausible-move generators are being tested. These look for captures 
for gain, checks, threads to capture for gain, retreating attacked units, 
and safe moves that avoid immediate loss of material. Plans are under- 
way to provide for handling defensive moves, interposition, development, 


and center control. 


Strategic heuristics are being debugged for complicating situations to 
exploit possible opponent's errors, and for controlling the use of the 
different generators in different situations. An claborate system of 


Macro-operations is used for defining the evaluators and generators. 


With debugging of current features, the program is expected to match 
low-to-average club players. j 


- 
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MATHLAB: On-Line Symbolic Computation - Carl Engleman (MITRE 


Corporation) 


MATHLAB is a LISP program that is intended to provide the services ~ 
of an on-line, time-shared computer to a mathematical scientist, such 
as a physicist, applied mathematician, or electronics engineer, for the 


daily performance of his tedious, mechanical, mathematical tasks. 


In its present state of development, the program is primarily oriented 
towards symbolic, rather than numerical, computation. Employing 
simple commands, such as ''solve" or ''differentiate', the program has 
facilities for symbolic manipulations, such as addition, multiplication, 
and differention of expressions and equations. It is also capable of the 
simplification of expressions and the substitution of one expression 


within another. 


Finally, it can, using the INTEGRATE program of M. Manove, succeed 
in the symbolic integration of certain expressions and the symbolic solu- 


tion of certain expressions. 
A detailed discussion of the MATHLAB program may be found in MITRE 


Corporation Technical Memorandum TM-04258, from Department 73, 
dated 26 May 1965. 
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INTEGRATE; On-Line Indefinite Integration - Michael Manove 
(MITRE Corporation) 


INTEGRATE is a LISP program designed to symbolically integrate 
rational functions in one variable over the field of rational numbers. 
Any rational function is admissible as input, and such functions may 

be represented as the sums, differences, products, quotients, or 
integral powers of constants, polynomials, and other rational functions. 
The antiderivative of a rational function will be expressed as the sum 
of a rational function with rational coefficients, and the logarithm of 


a rational function possibly with complex coefficients. 


The program will also integrate many rational functions of trigonometric 
functions, of the form f{(sin x, cox x, tan x, cot x, sec x, csc x) where 
"f" is a rational function, but this feature is limited by the fact that 


answers are always expressed in terms of "tan (1/2)x". 


The coding of INTEGRATE is divided into two major sections. The 

first section will find the rational part of the antiderivative and determine 
rational functions whose integrals are the summands of the logarithmic 
part. The second section will find the logarithmic parts of the anti- 
derivatives of a broad range of rational functions. This, of course, is 
subject to the size limitations of computer memory and other hardware 


variables. 


This program has been incorporated in the MATHLAB program of 

C. Engleman. A detailed discussion of INTEGRATE may be found in 
MITRE Corporation Technical Memorandum TM-04204, from Depart- 
ment 73, dated 22 April 1965. 
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Molecular Model Building - Cyrus Levinthal 


Many of the recent developments in molecular biology have depended on 
an understanding of the three-dimensional structure of large molecules. 
Such an understanding is usually reflected in the ability to construct a 
three-dimensional model of a particular molecule. For really large 
molecules, however, such construction is difficult and time consuming 
and, if complete, would entail an impossible enumeration of the many 
small interactions which contribute to the molecular stability. Thus, if 
a structure is known, physical models can be built, but construction 
complexity prevents the use of such models in the examination of a large 


number of different configurations. 


During the last six months we have written a set of programs with which 
we can construct, display, and analyse macromolecules on the MAC 
system's ESL display unit. With ‘computer-controlled display and real- 
time control of the molecule, an observer can obtain a true three-dimen- 
sional visulization of a particular molecule; and by interacting with the 
program which generates the molecule, he can indicate the way in which 


particular parts of the structure are to be altered. 


The first program written in this project calculates the coordinates of the 
atoms in a protein. Only those angles about which rotation is possible 
are entered as input variables in the program: all other rotation angles 
and chemical bond lengths are entered as rigid constraints. The calcula- 
tion treats each chemical bond along the linear peptide backbone as a 
translation and a rotation of a coordinate system attached to individual 
atoms. In this way, the step from one backbone atom to the next involves 
a translation and a rotation-matrix multiplication. The updating of the 
rotation matrix is determined either by the fixed rotation angles, which 
are introduced as constraints in the program, or as input data which are 
added by an investigator. The basic data for these constraints on bond 
angles and lengths have been obtained over a number of years from the 
X-ray crystallographic studies of Pauling, Perutz, Kendrew, and many 


others, and they are now firmly based on experimental information. 
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The three-dimensional structure of the protein is determined entirely by 
the physico-chemical interactions which lead to the lowest energy con- 
figuration. Since the bulk of the free enevgy is contributed by many small 
interactions, which for the most part have extremely short range, we 
should like to determine for each pair of atoms whether the two members 
lie close together in real space and thus contribute substantially to the 
interaction energy. Since even a small protein molecule may contain 1500 
atoms, the total number of pairs to be tested is very large. 

\ 
To avoid the enumeration of all possible pairs, we have written a set of 
programs that first divides space into cubes and then makes a list of all 
atoms in each cube. By searching through the cubes, and listing those 
pairs in which one member is in a cube and a second member is in the 
same cube or in one of the twenty-six surrounding ones, we are able to 
enumerate all pairs of neighboring atoms. This list-processing procedure 
finds the holes in a molecule by locating empty cubes surrounded by filled 
ones; and, in addition, it determines the inside and outside of the molecule 
by finding cubes which have filled neighbors on one side and empty neigh- 
bors on the other side. The knowledge of these topological aspects of a 
folded protein is necessary for the assessment of interactions of particular 
amino acid residues with the water normally surrounding all protein mole- 
cules. These programs used SLIP subroutines written by Professor 


Weizenbaum, whose cooperation was invaluable in our work. 


When an investigator wishes to alter the structure of a molecule, he inserts 
a pseudo-energy between a particular pair of atoms. An energy-minimizing 
routine, which looks at the sum of all interactions, including the pseudo- 
energy, calculates the appropriate change required for the variable angles, 
If the change leads to a violation of van der Waals contacts, another routine 
eliminates this violation. The latter routine, as well as those for energy 
minimization, require the calculation of the partial derivatives of the dis- 
tances between pairs of atoms with respeci to all of the variable angles 
which can effect these distances. We have used two procedures in altering 
the structure and eliminating van der Waals violation: the first uses one 

of several different minimization routines in a sequence of small steps; the 


second uses a linear approximation and calculates appropriate angular 


change by solving linear equations, subject to the additional constraint 


ee 
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that the sum of all angle changes should be a minimum. When we have 
tested and run programs and subroutines together in a smooth system, we 


can attempt to solve for the unknown structure of a real protein. 


By making use of molecular models constructed in this way, we can 
evaluate the 1mportance of electrostatic interactions in the stability of 
protein configurations. These interactions occur because of the existence 
of a large number of small electric dipoles. Although each of the dipole 
interactions is small, they have a longer range of interaction than the 

van der Waal forces and are so numerous that they contribute appreciably 
to overall molecular stability. The results of our calculations indicate 
that dipole interactions produce a stabilizing energy, if the helical regions 
of the proteins are antiparallel to each other, and a destabilizing energy, 


if the helices are parallel to cach other. 


The structure of the myoglobin protein molecule has been completely 
determined by X-ray crystallographic methods. When we used the approxi- 
mate coordinates of this molecule in our model, we found the net stabilizing 
effect of these dipole-dipole interactions to be significant. We have also 
used our model-building program to refine coordinates of the myoglobin BB 
model, by supplementing the X-ray data with data on the chemical bonds 
derived from other sources. Prior to this computer réfinement, our 
attempts to carry oul a calculation of this kind involved the construction 
and measurement of physical models made of brass wires. This calcula- 
tion was both tedious and inaccurate for a molecule as De INCE Oni 
We have written a program that uses Project MAC's PDP-6 as a display 
unit for structures whose coordinates are generated on the 7094. This 
program was designed to test the usefulness, as aids in the visualization 

of a three-dimensional structure, of perspective and modulation of line 
intensity with depth. It is written in such a way that one can test the 
feasibility of a three-dimensional visualization scheme that requires fewer 
rotation calculations than the ESL display. These tests will determine the 


required characteristics of additional off-line display equipment. 


We have written a set of programs which calculate the strain energy of 


displacement of chemical bond angles and lengths from their equilibrium 
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positions. When the rigid’ constraints on bond length and angles in the 
programs described previously are relaxed, the number of interactions 
in a protein becomes prohibitively large. However, studies of small- 
molecule structure can be carried out ina reasonable amount of computer 
time. We have used this procedure to refine models of the Sugars which 
are a part of nucleic acid molecules, and we are studying the effect of 
this on the overall structure of nucleic acids. Since we can carry out 
this refinement for any small organic molecule, we can, in principle, 
display the three-dimensional model of a small molecule when our only 
knowicdge of it is what atoms compose it and which ones are covalently 


bonded together. 


[Editor's Note: Figures 1 and 2 are stereographs of a DNA molecule and 
several protein molecules, that were photographed from the ESL Display 
Console at Project MAC. To ,use the stereo pairs without a viewer, 
ignore the double image print for the moment and focus your eyes ona 
point about 10 feet away. Bring the figure slowly up into your line of 
sight, maintaining the same focus of 10 feet. Four images will be within 
your field of view, and then the two middle images should merge intoa 
central image. You now have three images side by side. The center one 
will be the stereo composite of the two printed images. You may need 


practice to hold this image until your eyes focus on the stereo image. 


If you have difficulty with the multiple images, try holding your hand, 
edge on, in front of your nose, so that the left image cannot be seen by 
your right eye, and the right image cannot be scen by your left eye. This 
will form a crude ster coscope and allow you to concentrate ona single 


image, the stereo composite. 


Figure | shows three views of the same portion of a DNA molecule.’ 

These were produced by photographing stationary images that differ by 
approximately 10 degrees of rotation about the Z axis. Normally, a 
person at the ESL Display Console is able to obtain a solid-appearing view 
of the molecule by real-time rotation of the image about three axes. The 


views shown in Figure 1 are approximately 10 million times life size. 
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Figure 1. Stereo Pairs of a DNA Molecular Structure 
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Deoxyribonucleic Acid (DNA) is the carrier of genetic information in 
living systems, and consists of two phosphate chains that run in opposite 
directions and are coiled into a right-hand helix. Nitrogen bases, held 
together by hydrogen bonds, link the ester chains like steps of a spiral 


staircase and provide the sequence which carries the genetic coding. 


Figure 2 shows three views of complete and abbreviated polypeptide 
protein chains arranged in the a-helix configuration. The a helix, capable 
of forming hydrogen bonds between all of its amide groups, is a substantial 
constituent of protein molecules, which perform vital enzymiic functions, 


such as controlling chemical reactions, in living systems. 


The top view in Figure 2 is a fairly complcte polypeptide structure, 
while the other two views are simplified versions of the same molecule. 
The bottom image is the same molecule as the middle one, showing the 


exaggerated stereo separation obtained with approximately 20 degrees of 


rotation. | 
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( 
Figure 2. Stereo Pairs of a Protein Molecular Structure 
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Computer-Aided Teaching - Daniel Roos 


The Department of Civil Engineering has used the Project MAC time-sharing 
system in its undergraduate and graduate education programs. Three 
different types of educational activities relied heavily upon CTSS: class- 
room demonstrations, homework and term paper assignments, and an 


experimental time-sharing section of course 1.15. 


The time-sharing system has been used very successfully in classroom 
demonstrations to assist in the teaching of engineering. Ina classroom 
equipped with closed-circuit television, an IBM 1620 computer, CalComp 
and Gerber plotters, and an IBM 1050 remote console connected to the 
Project MAC computer, the instructor has all the computing facilities 
available to demonstrate points that he is making. If a student asks the 
perennial question, ''What happens if ?'' the teacher can immediately 


demonstrate the implications using on-line computer facilities. 


Several problem-oriented languages have been developed, with which an 
engineer can communicate with the computer in engineering terminology. 
Two of these languages, STRESS for structural problems and COGO for 
geometric problems are being used in conjunction with time-sharing for 
educational applications by both instructors and students. A student does 
not need intimate computer knowledge to effectively use these languages. 
With minimum instruction, generally one hour or less, he can use the 
problem-oriented language and the time-sharing system to obtain solutions 


to engineering problems. 
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The Structural Design Language - Robert D. Logcher, Gerald M. Sturman, 
Richard V. Goodman, Lewis C. Teague, Salvatore G. Mazzotta, and 


John R. Roy 


Structural design, the determination of structural elements and their layout 
and makeup, is a task that involves innumerable variables and decision 
parameters and requires extensive computation and decision-making. The 
nature of the computations, and the manner in which they are executed, 
makes classical forms of computer use inadequate and unsatisfactory. The 


purpose of the work described here is to determine the form and organiza- 


Structural design requires interactions with other technologies involved in 
the design product, such as mechanical and electrical engineering and 
architecture. The engineer must accept and recognize the conflicting aims 
of these technologies and evaluate compromises analytically. He is there- 
fore engaged in an iterative process in which he must remain flexible. The 
iteration consists usually of a series of triat solutions together with an 


intuitive searching process that leads to an optimum solution. 


Generally the amount of computation performed between decisions is small, 
but varies greatly with the problem, the engineer, and the various parts of 
the design task. Problems vary widely in size and complexity, from a few 
hundred variables and parameters to thousands. Depending on the problem 
and the engineer's sophistication, the design task may be automated with 
all decisions made prior to solution or may be broken into many processes 
with decisions inserted during solution. A process of design such as 
elastic analysis involves a great deal of computation; while others, such as 


checking performance against criteria, require very little. 


A time-sharing environment is required to join the creative design process 
to a computer. The user may thus perform small and large tasks in an 
unpredetermined order and express decisions between tasks. The problem- 
oriented language STRESS represents the first phase in computer-aided 
structural design. This language contains a data-input mechanism and 


programs for linear elastic-frame analysis. Since this analysis is the most 


difficult computational problem, but the best defined mathematically, its 
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solution was the first undertaken. General data-structure capabilities for 
engineering data were developed for STRESS which make its logical inclu- 


sion in a broader language possible. 


Over the past year, the Structural Design Language (STRUDL) has been 
developed in an open-end form. The major commands of this system are: 
1. PROBLEM INPUT 6. SELECTIVE OUTPUT 

2. PRELIMINARY ANALYSIS ?. WEIGHT DETERMINATION 
3. DETERMINATE ANALYSIS 8. CHECK PROCESS 
4. MEMBER SELECTION 9. SPECIFICATION CHECK 
5. STIFFNESS ANALYSIS 10. REVISE TABLES 
Each of these calls a particular Operation commonly used in the design 
process. The engineer can decide which operations are to be used and in 


what order they are to be executed. 


For example, PRELIMINARY ANALYSIS permits him to insert assumptions 
regarding the behavior of an indeterminate structure and thereby to make 

the problem determinate. He may then use DETERMINATE ANALYSIS to 
obtain a complete solution. MEMBER SELECTION, based on previous. 
design-criteria input and permanently-stored tables of standard members, 
may follow, and in turn be followed by STIFFNESS ANALYSIS, which pro- 
vides a rigorous solution for the indeterminate structure. With SELECTIVE 
OUTPUT, the user can obtain only that Output which is of immediate interest, 


rather than voluminous data which is actually generated. 


As modifications to the structure are made, either automatically or at the 
user's judgment, new analyses are executed until the engineer feels that an 
optimum design has been achieved. This type of on-line operation provides 
a new dimension to structural design and gives the engineer an ability for 
Optimization not previously possible. 

Seven of the commands listed have been coded and implemented. The 
remaining three, as well as others which further extend system capability, 


require further work. 
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Dynamic Structural Analysis - J. Melvin Biggs and Daniel Beltran-Maldonado 


In the design of Civil Engineering structures, the problem of analyzing 
those subjected to dynamic actions arises frequently. The laborious nature 
of the analysis on the one hand, and the possibility of more refined methods 
on the other, make the solution of these problems by a computer very 
desirable. Indeed, the technical literature indicates extensive use of 
computers in this connection, mostly through special-purpose programs 
coded to solve individual problems. The development of a programming 
System to analyzea variety of structures for free or forced vibration, such 
as Space frames, space trusses, plane grids, plane frames, and plan 
trusses, is in progress. This system is a problem-oriented language 
called DYNALS, similar in Many respects to STRESS. DYNALS is unique, 
not only because it is applicable to a wide range of structural configurations, 
but also because it encompasses a large part of the total design problem. 
Given the geometry of the structure and the stiffness and mass of individual 
elements, it forms the total mass and stiffness matrices, obtains eigen- 
values and eigenvectors, and computes dynamic response of the structure 
in minute detail for any specified distribution and time variation of applied 


forces. 


Three aspects of the problem make implementation in the time-sharing 
mode highly desirable: 

1. The iterative nature of the design process calls for a simple 
and rapid way of modifying a problem and analyzing the effects of such 
modification on the structural response, 

2. The wide variety and great volume of possible output available 
from a dynamic analysis make the ability to request at will selected forms 
of output very desirable, 

5, Time-sharing offers unique advantages in the debugging and test- 
ing process, which are quite valuable to a user who is working with sucha 


large and complex program block. 


\ DYNALS will be incorporated into STRUDL, which deals with a broader 


portion of the design process for Civil Engineering structures. To make 


the system efficient for both large and small structures, the various data 
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arrays are dynamically allocated so that only needed arrays are required 


in core at any given time. 


DYNALS accepts and decodes data pertaining to the continuous structure 
and performs a ''discretization'! of the stiffness and mass characteristics. 
The discrete system consists of a number of degrees of freedom, equal 

to the number of kinematic degrees for a joint, times the number of joints. 
The analyst, presented with various alternatives regarding the formation 
of the mass matrix, chooses with regard to the nature of the problem at 
hand. He may accordingly reduce the number of degrees of freedom to a 
specified selected set, A wide variety of output, available upon request 
from the console, ranges from eigenvalues and eigenvectors of the free 


vibrating structure to the complete response of one and all of the mechani- 


, cal or geometrical quantities of the structure while it vibrates under a 


er, é 
time-varying load system. 


DYNALS may be subdivided into the following programming blocks: 

1, INPUT phase, decoding of input statements, storage of data 
and process parameters, and consistency checks, 

2. Formation of the logical stiffness matrix, 

3. Formation of the logical mass matrix for the case where member 
masses are specified, 

4. Reduction of the number of degrees of freedom to a selected set 
when specified, 

5, Reduction of loading data, 

6. Solution of the eigenvalue problem, 

7. Using the characteristic shapes obtained in 6. as input joint 
displacements, the preparation of sets of static solutions for forces, 
reactions and.member distortions, (n is the number of modes to be 
combined. ) 

8. Using the results from 5. and 6a, the numerical integration of 
the modal equations of motion, the results of which are used to combine 
the static solutions developed in 7., (Record is kept of maxima and time of 
their occurence for forces, reactions displacements and distortions. ) 

9. OUTPUT phase. 


As of June 1965, the system design and a substantial portion of the program- 


ming is complete. 
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Soil Engineering Problem -Oriented Language - Robert L. Schiffman and 
i tiented Language 


Laurence N. Beckreck 


SEPOL (Soil Engineering Problem -Oriented Language) is a system of 
interrelated Programs which serves a three-fold purpose. First, SEPOL 
is designed as a communication system in which the communication between 
man and machine is in an interrogation form that elicits the program from 
the user and channels the program towards the desired Objective. Second, 
it is a computer-aided design system that so widens the area of choice 

and breadth of calculation that by sheer speed and utility, a greater range 
of calculation is Open to the soil engineer. Third, it performs a basic 
Study of the nature of soil engineering analysis. By subdividing the system 
and delineating the areas of common logic, it identifies the basic hypotheses 
and analytical units. 

SEPOL I, covered in this report, is a system designed to calculate, via 
time-sharing, the magnitude and progress of settlement of an earth mass, 
when the soil surface is subjected to a specified loading. It was designed 
as a prototype to determine the feasibility and usefulness of this type of 


computing system, and also to provide useful information for soil engineers. 


SEPOL I is composed of many small subroutines, each of Which represents 
a unique task that is always performed as a unit and can be performed any- 
where and many times in an analytical sequence. The control function of 
SEPOL I is interrogatory. A set of questions is posed; and, since the 
response to the questions provides the linkage between boxes, the set of 
questions and answers builds, for each Situation, a specific main program, 
SEPOL I was designed as an easily-expandable System. To adda new 
method of analysis requires first the Programming, in FORTRAN, of tasks 
or calculations not already in the system, and then the augmentation of a 
dictionary of variables, Access to the new System is accomplished by 


adding choices at the appropriate decision points. 
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1/0 System Research - Ken Reinschmidt, DaAN Gardner, and John W. Weber 
The hardware which connects the IBM 1620 to the MAC time-sharing 
system consists of two logically-different parts. First, a general I/O 
interface to the 1620 makes it easier to connect non-IBM devices, and 
second, a dataphone interface which goes between the general interface 
and the dataphone which is connected to MAC. Because of certain idio- 
syncrasies in the 1620 1/O control, such as the lack of interrupt, time- 
sharing operation must be under software control. Subroutines, written 
in 1620 assembler language and available to the FORTRAN programmer, 
allow him to use the 1620 typewriter and/or card read/punch as a CTSS 
console, pass numeric data between programs in 7094 and the 1620, pass 
alphanumeric data, and perform alphanumeric comparisons on data in the 
1620. The subroutines were developed in a modular form which facilitates 
modification or expansion of the present system. Mr. Weber also 
developed a supervisor subroutine for the 1620, which gives a more auto- 
matic operation. The main uses of the system have been for obtaining 
plotted results of data developed at the ESL Console and for card I/O of 
lengthy I/O operations. 


Another improvement on the 1620 was made by the addition of the Gerber 
VP-600 plotter as an output device. Tnis plotter was put on-line with the 
1620 in August, 1964, through a modification of the old CalComp Plotter. 
Several plotting subroutines were developed for the Gerber unit: line, 
number and character plotting routines by Mr. Gardner; and then curve 


and circle routines. 


An on-line, date-acquisition system has been developed jointly by members 
of the System's Laboratory and the Structural Model's Laboratory. The 
hardware part of the system, connected to the 1620 through the general 
interface mentioned above, consists of electro-mechanical switches, 
digital logic, signal conditioning and amplification, and an analog-to- 
digital converter. The system allows the 1620 to switch to any one of 50 


strain gauges and read the converted voltage into storage. i 
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A large set of software, which facilitates the use of the system, isa 
problem-oriented language for the data acquisition system. With this 
language, the user can do any of the following: 

1. Balance any or all gauges, 

2. Select how many and which gauge he desires to read, 

3. Read each gauge as many times as desired, 

4. Do simple calculations, 

5. Convert his data to several different forms, 

6. Plot results, 

7. Store intermediate results on the disk, 

8. Change loading and repeat operations. 
The system was developed by Mr. Reinschmidt, has been used in class- 
room demonstrations and on some plate studies, and is now being used 
extensively on models of a dam under construction in Spain. 
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Bridge Design - John F. Brotchie, Charles A. Cornell, Jose M. Roesset, 
and Guillermo Guzman Barron-Torres 

The purpose of this research is to develop an integrated and comprehensive 
system for the directed design and drafting of highway bridge structures. 

It is similar in several respects to the Structural Design Language reported 


elsewhere, but is oriented specifically toward the design of highway bridges. 


The design system proposed consists of a series of logical blocks, self- 
sufficient in themselves. These blocks, each one consisting of a large and 
expandable number of links, are to be interconnected and each is amenable 
to direct input and output. Individual blocks perform the following: 

1, An internal definition of geometry; 

2. A determination of suitable alternate forms, span arrangements 
and materials; 

3. A design of superstructure, piers, abutments, approaches, or 
parts thereof; 

4, An analysis of the structure or components above; 

5. Optimization; 

6. Detailing; 

7. Cost analysis and calculation of quantities to various degrees of 
detail, depending on the stage of the design at which the link was entered; 

8. Documentation, including drafting and specification preparation, 
The point of entry to the system, the flow between the blocks, and the kind 
of exit may be chosen by the designer at the beginning of the process, but 
are preferably chosen during it. Time-sharing facilitates this choice, and 
a large computer facility provides the storage necessary to retain all 
required information throughout the process and communicate it to the 
designer in the most efficient way, whether by plotter, printer, teletype 
or oscilloscope. Thus, the logic of the design process and the basis for 


determining design parameters is controlled by the designer. 


Some phases of the system are already completed and were recently linked 


by Guillermo Guzman to demonstrate a time-sharing version of the program. J any 


These were the input phase, the definition of internal geometry, the selection 
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of alternates for straight bridges, and the design of simple-span, 
composite, precast-prestressed, concrete superstructures together 
with their optimization on the basis of a comprehensive niinimum cost. 


A pilot version of the steel design phase will also be added shortly. 


Optimal Synthesis of Road Networks - Alan M. Hershdorfer 


This study was originally an investigation into the properties of optimal 
traffic flo': patterns on metropolitan road networks. A total-travel-time 
function to be minimized is expressed as a piecewise, linear, convex 
approximation to the experimental travel-time vs. traffic-volume relation- 
ship. The linear programming code RSMFOR was used to determine cost 
flow solutions for a multi-commodity flow system. The model was then 
extended to a mixed-integer programmiing formulation. . The integer 
variables represented decision variables related to the.provision or 
deletion of road capacity. In this model, numerous related problems of 
traffic network synthesis are studied. The function to be minimized is a 


linear function of construction cost and user cost. 


The simplest synthesis problem, that in which construction cost is nil, 
Occurs in the problem of optimally orienting a set of one-way streets. 
Integer-linear constraints on lane allocation require all lanes in a road 
link to be oriented in one of the two directions or to be split equally between 
the two directions. This one-way street problem was studied extensively 
with 4 time-sharing version of the Land and Doig Algorithm for mixed- 
integer programming, which prescribes values for integer variables in an 
ordered sequence of linear programs. With time-sharing, the algorithm 
was monitored during computation and, thus, trial solutions were very 
rapidly evaluated. Since the structure of the flow constraints forced 
certain variables to near integer values after a small number of integer 
decisions had been made, a good primal feasible integer solution was very 
rapidly determined. Time-sharing, used as a monitor on the algorithm, 
has enabled us to develop a heuristic method for quickly obtaining solutions 


to problems that involve the synthesis of flow networks. 
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Sequences and the Four-Color Problem 


A Design Language for Digital Systems 
Man-Machine Communication 
Research on the Theory of Automata 
Program Segmentation 
Waveform Transformation and Graphical Display 
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On-Line Braille 
Analysis of Time-Shared Computer Systems 
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Sequences and the Four-Color Problem - David A. Huffman 


Methcds for the generation of binary-valued pseudo-random sequences, 
which are more general than shift-register sequences, have been studied. 
These sequences have the interesting property that they remain invariant 
to a type of transformation which, in a manner reminiscent of tag 


systems, substitute one string of symbols for another. 


The four-color problem (a famous, as yet unsolved problem in topology) 
has been studied with the view that it furnishes a difficult and realistic 
example of a problem for which a man, aided by appropriate program- 
generated computer displays, may be able to gain insight more rapidly 
than a man working with pencil and paper. The maps associated with this 
problem are nicely suited to representation by certain simple, symmetric 


list structures. 


A Design Language for Digital Systems - Gerald J, Burnett and 
Jack B. Dennis 


This research was promoted by the need for a more consistent and con- 
venient method of designing digital systems and by the desire for computer 
aid in this design process, We have developed a design language to meet 
a number of requirements. Some of the more basic requirements are the 
following: 

1, The language must be complete, that is, it must be able to 
describe the logic and timing of any digital system; 

2, The language must be convenient for the designer and must be 
simple enough to express all of the common hardware logic and timing 
functions: 

3. The language must have features analogous to macros and 
repeats in assembly languages in order to describe readily-iterated logi- 
cal structures; 

4. The language must be able to describe a system in sections or 


components that are designed by different individuals. 
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As an example of the description of a familiar logical structure in the 
design language, we describe the lugic required to perform parallel 
binary addition of a number contained in register B to a number con- 


tained in register A and leave the result in register A. 


The addition method used for the example involves two steps. When the 
pulse named “partial-add''occurs, cach bit of register A is complemented 
if the corresponding bit of register B is one. The pulse named "carry" 
complements a stage of A if there is a carry signal into that stage. The 


description of one stage of the adder is as follows: 


begin partial-add: BL i] tALil; 
carry: if C[i] then TA[i]; 
Alot | a Git | AAI SULA BLil; 


end; 


The combinational logical operators A, V and 7 have their usual signifi- 
cance, while the upward arrow denotes the action of complementing, as a 
unary operator, or of conditional complementing, as a binary operator. 
The first two statements of the example show ways of indicating that cer- 
tain actions are to take place on the occurrence of definite pulses, either 
absolutely as in the first statement or conditioned by a logic level as in 
the second, The third statement illustrates several features of the 
‘language. The where clause indicated that the statement following then 

is to be considered part of the description only if the Boolean expression 

i #0 is truc. In this instance, the following statement represents combina- 
tional’ logic that generates level C[i-1| | according to the expression at the 


right of the equivalence sign, 


The complete dese ription of the adder involyes iteration of the above 
logic for each position in the two registers and is handled by a for 
statement as shown below. Certain declarations are included so that the 


types of named objects are unambiguous. The integer quantity n is the 
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number of bits in registers A and B and is presumed to be assigned a 


definite value outside this portion of the description. 


integer n; 
register A[0, n-1], B[0, n-1]; 
level + C[0, n-1] 


add: partial add; 
delay 500n; 
carry; 

C[n-1] =0; 


for i = 0 through n-1 then 

begin partial add: B[i] f A[I]; 
carry: if C[i] then f A[i]; 
where i #0 then 


Cli-l1] =Cli] A Ali] v n Cli] ~a Bli]; 
end; 


The pulse add performs the partial addition and initiates a 500-nano- 
second delay, At the end of the delay, the carry pulse occurs. The 
meaning of the for statements is intrepreted as if the (compound) state- 
ment foliowing then werc written out successively for integral values of 


i from 0 through n-1, 


Work related to the design language will be devoted next to the develop- 
ment of suitable simulation techniques and means whereby a designer can { 
employ the language on-line in the process of creating digital system 


plans. 
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Man-Machine Communication - Joseph Weizenbaum 
Seen Vommunication 


SLIP 
SLIP, a list-processing system capable of being imbedded in a large 
class of algebraic compiling languages, was imbedded in the MAD system. 
MAD is the compiling system most frequently used in CTSS, and there- 
fore, SLIP has been made available to large scctions of the MAC user's 
community. A number of students and faculty members have made use 


of SLIP for research and thesis purposes, 


OPL 


The on-line SLIP-based Programming system OPL was rewritten in the 
MAD-SLIP system discussed above. It has served as atest bed fora 
number of experiments, particularly text manipulation, and hae been an 
invaluable tool in clarifying issues underlying interactive programming 
systems, their design and implementation. The problem of the program- 


ming system-user interface has been illuminated, but not solved. 


ELIZA 3 

A major effort in the design of a natural-language, on-line, man/machine 
conversational system was begun and carried forward during the report- 
ing period. At this time, ELIZA is capable of conversing with users in 
natural language in certain restricted contexts. We are beginning to 
expand ELIZA so that it can form the basis for both information gathering 
and a retrieval system Operating on a natural-language basis, ELIZA 
was also used for experiments on two-person conversational interaction 
at the Massachusetts General Hospital, The development of ELIZA 
required implementation of powerful text manipulation functions within 


MAD-SLIP. These latter are of general use as external MAD functions. 
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Research on the Theory of Automata - Robert McNaughton 


A, CLASSIFICATION AND COMPLEXITY OF FINITE-STATE OPERATIONS 
The most commonplace distinction in the theory of automata has been 

that between finite-state operations and operations requiring potentially- 
infinite automata. Partial motivation of our research has been the 

contention that finite-state operations should themselves be classified in 
some interesting manner, Stimulated by some work of M. P. Schutzenberger, 
we have proposed a notion that we think gives an important classification, 
namely the notion of a "topological" event (a term which we do not feel is 


optimal). 


Interestingly, the notion can be characterized in several distinct ways: 

1) as an event whose reduced-state graph contains no loop that counts any 
word modulo m, for m>; 2) as an event expressible as a regular expression 
containing free use of all the Boolean operators, but without star (denot- 
ing closure, or iteration); 3) as an event expressible in a certain first- 
order language of symbolic logic; 4) as an event realized by a nerve net 
in which feedback is confined to the case in which output of a neuron 
drives one of its own inputs; and 5) as an event whose semi-group 
contains no non-trivial subgroups. The fact that this class of events can 
be characterized in so many different ways - structurally, behaviorally, 
and algebraically - each interesting in its own right - is evidence enough 
of the importance of the concept. A long paper on this matter is planned, 
which will include proofs establishing the equivalence of the five charac- 


terizations, 


There are two interesting ways of defining the loop complexity of a 
regular event, each in terms of a class of regular expressions. The 
first class is that of restricted regular expressions, i.e. those made up 
from the operators union, concatenation, and star. The second class 
consists of those made up of these operators and, in addition, inter- 
section and complementation; such expressions are called "general 
regular expressions". Thus general regular expressions contain all the 
Boolean operators, while restricted regular expressions contain only 
union, in each case, the loop complexity of an event is defined as the 


minimal star height of all regular expressions in the class that represent 
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the event. (The star height of a regular expression is the maximum 
length of a nested chain of the stars, ) Thus, for example, the loop com- 
plexity of a topological event with respect to general regular expressions 
is zero, But there are topological events with arbitrarily large star 


height with respect to restricted regular expressions, 


It seems that the most useful tool for investigating star height with 
respect to restricted regular expression (which is the problem studies 
originally by L. Eggan) is graph theory, On the other hand, the most 
useful tool for investigating star height with respect to general regular 
expressions is the theory of finite semigroups. Both concepts seem 
interesting and important, but both have unsolved problems. Thus, we 
do not even know whether there exists an event of star height 3 with 
respect to general regular expressions; and, although we do know the 
existence of events with arbitrarily large star height with respect to 
restricted regular expressions, we do not know whether there exists an 


algorithm for determining this star height when a regular event is given. 


B, CONSTRUCTION OF WELL-TIMED NETS FROM BADLY -TIMED 
ELEMENTS 

Research in this area has been a continuation of research at the 

University of Pennsylvania. A seminar was held consisting of some 

members of the Machine Structures Group of Project MAC and some 

design engineers from Honeywell. The number attending was small, but 

the discussions were intensive and fruitful. Some improvements were 


made in the report written at the University of Pennsylvania (Badly-Timed 


Elements and Well-Timed Nets, by R. McNaughton). This report will be 
a A eee INE 


issued as a monograph next year. 


C, INFINITE HISTORY OF FINITE-STATE MACHINES 

Two widely-differing studies in the literature have led into this area. 

J. R. Buchi studied a certain system of symbolic logic, and developed a 
decision procedure by showing, in effect, that every formula in the 
System says something about the ihfinite history of a finite-state machine. 
And D. Muller analyzed the behavior of an asynchronous circuit by regard- 
ing the circuit, which is subject to varying input conditions, and in which 


the relative timing of the parts is unpredictable, as a finite automaton 
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with an infinite duration in time: a buzzer, for example, would be con- 
ceived of as a two-state automaton going from one state to the other ad 


infinitum. 


The proof of an analogue to the so-called Kleene Myhill theorem for the 
infinite-duration problem has been written up, and will be submitted for 
publication. Also, there is an interesting connection between this 
problem area, and the topic of infinite games; i.e., games which two 


players play in turn an infinite number of times. 


It would appear that this problem area is more abstruse than the other 
problem areas. In spite of the fact that the area was approached by two 
widely separated studies, no new application seems to be forthcoming. 
We are far from the position that this area is worthless; however, we 


do think it deserves less attention than the other two. 


Program Segmentation - Richard Y. Kain and David J. Kuck 


A question of major importance is segmentation of users programs to 
decrease memory space and swap times. There are several problems 
which immediately arise. The major one relates to efficiency of the 
algorithm which might be developed to perform such segmentation. In 
particular, is the time spent segmenting a program greater than the 


time saved? 


Four approaches to this problem are: 
‘1, The compiler examines a program's gross structure and 
makes simple decisions about where to divide the program; 

2. The compiler divides a program into small "atoms" and 
then attempts to combine the "atoms" into segments; 

3. The compiler divides programs into small segments, and 
then adds to each segment a small program which determines those 
segments which might be brought into memory at the same time; 

4. The compiler describes possible paths through a program in 
terms of regular expressions, which.are examined when allocation 


decisions are made. 
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The first approach is difficult to implement without some subdivision of 
a program; therefore, it will be considered as a special case of the 


second method. 


The second method "reduces" the segmentation problem te a clumping 
problem: given a set of atoms, which are most closely related to eagh, 
other? A number of approaches to this problem have been considered, 
but all unfortunately assume that there is some way to measure "inter - 
relationships'' between pairs of atoms; thus reducing the interrelation- 
ship to anumber. We will assume that these numbers have been placed 
in a matrix, A, where an is the interrclationship between atoms i and 
j. Some clumping strategies are: 

a) Build clumps by adding ''most suitable" neighbors until 
some criterion of maximum size or minimum suitability is met; 

b) Partition matrix A until the sum of the clements not in 
the matrices along the main diagonal is small; 

c) Take random walks through the ''atoms,"' letting the 
probability of a path be proportional to the interrelationships of the 
elements at the ends of the path, with termination by length of the walk; 

d) Find the largest cigenvector, whose components are ane 
probabilities of visits to the corresponding states in the long run, and 


clump according to the size of the components of the cigenvector, 


The third approach is a crude attempt towards "dynamic segmentation. "' 
By including a simple program which determines neighbors (possibly by 
examining the parameters of the running program, such as loop limits) 
with each atom, only those atoms which are likely to be used during the 
next quantum need be fetched. If the decision is too complex, a standard 


set of atoms would be used. 


The fourth approach is another attempt, less sophisticated than the third, 
to dynamically segment the program. Program data are not checked in 
the choice of atom set and less cconomy is achieved. This method is 
superior to a priori segmentation, since the beginning atoms in a segment 
might not be needed if the previous execution had halted in the middle of a 


segment. 
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Progress has been made in finding interrelationships among methods 2a 
through 2d, though no useful generalizations are immediately possible. 
Unfortunately, evaluation must depend upon detailed experience with 
each scheme or simultations which are necessarily very dependent upon 


the models of user programs. 


Waveform Transformation and Graphical Display - Thomas G. Stockham, Jr. 


A subroutine for the plotting of numerical data on the E.S. L. display was 
completed by Henry Ledgard. The plotting process gives freedom from 
details of scale-factor selection and, at the same time, produces graphs which 
are labeled decimally in convenient round numbers and may have logrithmic 
or linear scales on either axis. Subprograms that perform various sub- 
tasks for the logical specification of completed plots are available as a 
separable package. These routines are free from constraints imposed by 
the specific display unit for which they are intended and, thus, provide a 
general service for any plotting equipment. Their main subtasks are 
generation of numerically-convenient plot boundaries and associated 

scaling parameters, generation of normalized grid-line parameters, 
generation of text for dimensioning, and scaling of data to fit the normal- 
ized plotting space. Finished plots are produced from normalized, 
floating-point array pairs of arbitrary length. The total generation time 


is between one and two seconds on the IBM 7094, 


An algorithm has been developed by W. R. Chiodi with sufficient flexibility to 
implement a variety of linear transformations, such as Fourier, Hilbert, 
and convolution transformations, H. G. Murray, Jr. has realized the 
proposed algorithm in a computer program. This system operates ona 
standard-parts technique by which a required variability is concentrated 
in two building-block functions, two simple sets of rules for combining 
them, a set of parameters, and a simple transformation on these parame- 
ters. Since these variables can be specified through suitable numerical 
tables and brief sets of algebraic statements, the bulk of the computation 
is provided by a static framework into which they are embedded. This 
qwork will become part of a problem-solving system in electrical linear- 


system theory. (See Chiodi, and also Murray, Appendix B.) 
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A Table-Directed Translator - Chung L. Liu 


Research in the table-directed translator system has continued for the 
past year. In order to further investigate and evaluate some of the ideas 
conceived in the course of our research, such a system is being designed 


and implemented. 


It is also hoped that this system, aver completed, will extend a great 
deal of our programming powers. While a table-directed translator 
system provides users with the flexibility of designing their own program- 
ming languages, it is most important from the users! point of view to 
learn to use it with ease. This is what we stress in our design, Clear 
boundary lines between modules of the system and simple input languages 


are emphasized. 


It is expected that most parts of the sys#em will be programmed ina 
higher-level language, so that reprogramming work, due to hardware 
transition, will be minimized. We intend to use the MAD language for 


this purpose. 
On-Line Braille - Edward L. Glaser 


It should be understood that Braille is not a simple one-to-one code with 
normal characters of the English language. Rather, it is reminiscent of 
shorthand in that many standard contractions and abbreviations are used, 
Although these contractions and abbreviations are not phonetic, the rules 
for their use are often context dependent. As a consequence, a computer 
is necessary to convert standard teletype code to a modified teletype code 


that is adequate to drive an electric Brailler. 


Although this could have been done within the IBM 7094 computer at 
Project MAC, it was decided to use a separate satellite computer, so 
that Braille would also be available in the earliest phase of the new 
MULTICS time-sharing system on the General Electric 635 computer. 
As a consequence, a Digital Equipment Corporation PDP-8 computer was 
acquired this year, since it is a small, high-speed, versatile machine 


that is completely adequate for the job. To produce Braille output, we 
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selected an electric Braille Writer that was designed and developed in 


the M.1,T. Mechanical Engineering Department. 


The program for the conversion from teletype code to Braille has been 
written and checked out. Also, the interface between the PDP-8 and a 
Data Set has been checked out so that the Braille system can communi- 
cate with the MAC system as though it were a standard teletype station, 
Still remaining are final check-out of the interface between the PDP-8 
computer and the electric Brailler and check-out of the Brailler itself. 


(See Luconi, Appendix B.) 


Future plans include use of the PDP-8 machine as a base for additional 
studies of man-machine communication. Some of the possible areas of 
investigation are: auditory coding, interaction of auditory and visual 
stimuli, and visual and tactile stimuli. To date, a method for producing 
arbitrary auditory wave forms has been designed and fitted to the machine. 
It consists of two 6-bit digital-to-analog converters that are driven by the 
machine, The output analog wave form is: sent through low-pass filters 

to stereo amplifiers and loud speakers, Programming the equipment for 


this application has not been completed. 


Analysis of Time-Shared Computer Systems - Allan L. Scherr 


Some aspects of the operation of time-shared, interactive computer 
System: were analyzed by making extensive user and performance 
measurements of the Project MAC system. Emphasis was on the 
reaction of hardware systems to the demands that users make, Simply 
stated, the problem was to characterize both time-shared systems and 
their users in order to predict the performance of the two Operating 


together, 


Portions of the problem included Specification and measurement of user 
characteristics, the development and verification of both simulation and 
mathematical models for time-shared systems, and the specification and J > 


measurement of performance metrics for such systems, 
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First, simulation models were used to study the effects of changing small 
details in the operation of CTSS-like systems. Then, a continuous-time 
Markov process model was derived to predict the performance of a broad 
class of systems. Three lengthy programs were written to gather CTSS 
data. The first measured interaction parameters; the second measured 
input, output, and idle times of the system; and the third, statistics for 
interaction parameters of five system commands. This data was gathered 


on the MAC system between October 1964 and March 1965. 


The CTSS data were used as a basis for comparison with model pre- 
dictions, and, in order to take measurements and build models, many 
definitions of commonly used time-shared terminology were made precise. 
This research was used as the basis for a doctoral thesis. (See Appendix 
B and MAC-TR-18. ) 


Queueing Models for File Memory Operation - Peter J, Denning 


A problem which has received much attention at Project MAC is the 
scheduling of user programs. The Scheduling Algorithm assigns each 
waiting program a burst of processor time during which it runs; and each 
program continues to receive bursts, at intervals, until it has run to com- 
pletion. While not running, it is stored on adrum or disk. The process 
of placing a program in core memory or of removing it from core and 
placing it on the drum or disk is called swapping. An optitnum Scheduling 
Algorithm results in the most efficient system operation and yet treats the 
users fairly. Now that computing systems are about to be built in which 
segments of many users’ programs are simultaneously present in core 
and several processors operate simultaneously, another problem arises; 
scheduling of the necessary swaps and of any requests generated by the 
various user processes for file memory use. At present, requests are 
made singly to access file memory; and a request does not occur until the 
preceding cne has been processed. In a system with many simultaneous 
processes, requests for file memory use will be placed ina queue. This 
research attempts tofind a reasonable model of the way in which user pro- 
cesses demand file memory use and a reasonable and optimal method of 


handling the queued requests. 
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It is not possible to construct a detailed model for a computational pro- 
cess in a segmented, multi-programmed System, since no such system 
exists. However, we can make reasonable, probabilistic assumptions: 
for example, that segment lengths are Poisson distributed about some 
mean number of pages; that inter-request times are Poisson distributed, 
and that a given process will be able to generate requests while its write 
requests are in service,. but will be suspended while any read request is 


in service, 


Having constructed a reasonable model for user processes, we had to 
consider its request for file usage. The methods of queueing theory were 


applied, in particular, the following queue disciplines: 


1. first-come-first served, 
2. shortest-job-first, 


3. shortest-access-time-first. 


The access time, defined as the rotational positioning delay of a drum or 
disk file, is measured from the time a channel becomes free until the 
desired starting location of a request is in position; it is meaningless if 
there are no requests waiting in the queue. A shortest-access-time- 
first queue was found to be the most efficient of the three queue disciplines, 
Or any combination of them. Expressions were derived for the average 
wait in queue and the average number of waiting requests, and the niost 
efficient queue discipline was taken to be the one for which over-all idle 
time was minimized. Over-~all idle time is the idle time of all user pro- 
cesses, plus idle time of all processors, plus idle time of all file memory 
channels. A method was assumed to exist for optimal core memory 
allocation, that is, the allocation problem for core is irrelevant to the 
discussion of file memory scheduling, although we considered allocation 
schemes. In the analysis, simulation was found to be a most useful tool, 
because it showed which assumptions were valid and which were not. It 
also showed the form of probability distributions for waiting times in 
queue, average number in queue, and so on, that should be expected. 
Finally, simulation was used to verify the predictions of the mathemati- 


cal model. (See Denning, Appendix B.) 
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Input/Output Subsystems - Arthur A. Smith 


The execution of input/output instructions in a multi-access computer 
system may he considered as either simple ordinary instructions, with 
long execution times when the input/output equipment is part of the central 
processor, or as a form of communication between the program being 
executed by the main arithmetic processor and a program written 
especially for the 1/0 subsystem involved, when the I/0 equipment is a 
separate entity. The former consideration ignores the possibility and 

the desirability that the central processor execute an entirely separate 
program while the I/0 is occurring. Execution is implicit in the second 


viewpoint, which therefore is considered more fruitful. 


My studies regard a process as able to force its execution by a particular 
processor, in particular by a processor suited to the properties of the I/0 
device, That is, upon execution of a certain instruction, an image of the 
state word of the process is transmitted to the 1/0 processor, which 
modifies the state word by the execution of instructions written in its own 
order code, This I/0 processor, by executing an instruction analogous 

to the one which started its operation, causes resumption of the usual 
sequence of instructions, with the appropriately-modified state word. 
Thus, it executes an instruction which forces che subsequent instructions 
to be executed by a main (arithmetic) processor. I plan to publish a 


' . ! report on possible realizations of this system and alternatives to it, 


\ Non~Repeatability of Multi-Process Computations - Earl C. Van Horn 


A MAC system consists of several processors, al main memory, several 
file-storage devices, and several input/output devices. The supervisor 
. program of such a system schedules available processors among the 
various jobs to be done. We have defined the notion of a process as a 
precise molding of the job to be done. A process is a locus of control or 


activity within an instruction sequence. We have been studying the 


a 


relationships among processes and various other system entities and have 
concentrated on those relationships that are invariant with respect to the 
way the supervisor schedules available processors among the processes, 


A collection of processes working together for a single user ona 
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particular problem forms a separate computation. All of the processes 
of a separate computation may reference the same set of system entities. 
This set constitutes the sphere of protection associated with the computa- 


tion. 


The possibility exists in a MAC system that different runs of the same 
computation from the same initial state with the same input data will yield 
different output data. Such a computation is said to be non-repeatable. 
Computations running in the present MAC system, CTSS, are repeatable 
because they always consist of a single process. An attempt to adapt 
CYTSS for multi-process computations would encounter the problem of 
non-repeatability. Because in practice no program is ever debugged, 
repeatability cannot be guaranteed to a user unless sufficient conditions 
for repeatability are enforced by the hardware and supervisor of the 
system. These sufficient conditions may be published as a set of rules 
for users to. follow in the writing of programs or they may be followed by 
a compiler ia the generation of multi-process programs. In such cases, 
the burden of establishing repeatability lies on the person debugging the 
program‘or compiler, respectively. In order to study the interactions 

of processes in multi-process computations, we have defined a Multi- 
Process Automaton of Type A as a model of a multi-process computation, 
We have postulated one set of conditions sufficient to make a run of such 
an automaton repeatable. A rigorous proof of this hypothesis is forth- 
coming. Later, we hope to formulate and prove the necessary conditions 


for repeatability in such an automaton. 


Semantics for Multiprogrammed Computations - Jack B. Dennis and 
ce ESS 


Earl C. Van Horn 


The description of computations performed by a multiprogrammed computer 
system requires functions not available in existing programmiing languages. 
In particular, functions related to parallel programming, protection of data 
and procedure information, and communication among independent compu-~ 


tations are usually absent. We have studied the problem posed by these 
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semantic aspects of computations within time-shared computer systems 

and have proposed a set of ineta-instructions through which these essential 
facilities may be realized. Some of the meta-instructions defined already 
exist in multiprogrammed computer systems as monitor call operations; 
however, many of the meta-instructions are yet to be made part of any 
computer system. The operations performed by meta-instructions are 
based on the program structure developed in our previous work. A com- 
putation is defined by a collection of capabilities represented by entries in 

a C-list. The process active in a computation may invoke capabilities = % 
represented in its C-list, for instance, the addressing of a procedure or 


data segment or the executing of an I/0 function. 


A group of the proposed meta-instructions relates to the control by one 
computation of other inferior computations. This hierarchical relation- 
ship of computations is essential to program debugging within time-sharing 
systems and to the proper handling of faults or exceptional conditions 
encountered by a computation process. Other meta-instructions relate to 
the accessing and control of a directory structure for retained objects, 
principally segments of procedure and data. We suggest a hierarchical 
directory structure associated with each intelligence that uses the computer 
system in order to facilitate the choice of-unique identifiers for retained 


objects and make the sharing of objects among computations feasible. 


Optimal Allocation of System Resources - Robert L. Potter 


The problem of optimal memory and processor allocation in a multi- 
processor, time-shared computer is being studied. It is formulated asa 
two-space mapping problem with emphasis on the use of probabilistic 
methods. A computer has a certain number of memory and of processor 
units available. A two-dimensional space can be visualized. I constructed 
two'spaces with the same normal axes, a memory unit axis and a processor 
unit axis; one is called ''resource'' space, the other, ''needs'’ space. The 
resource needs of users are distributed in needs space. An allocation 
algorithm assigns each user to a portion of resource space. This isa 


mapping from needs space onto resource space, and the mapping function 


must assign resources equitably and efficiently. 
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The probability densities of memory and processor needs per user will 
not be known until measurements can be made on an appropriate system. 
We expect the densities, which are time-dependent, to have the same 
general appearance; they are zero at zero need, rise to some maximum, 


and then fall off toward zero as the need becomes infinite. 


The mapping function is not easy to determine; it must include cost, 
efficiency, and user-satisfaction factors reflected in the way distribution 
in needs spaces adapts to the system. Cost factors may be defined in 
many ways: as the allocation scheme which minimizes cost, in one case; 
or wasted time or resource units in another. User-satisfaction factors 


are elusive of definition until a required system is built. 


The mapping must be time-dependent, because the present MAC system 
considers allocation primarily in terms of program length, which dis- 
criminates against long programs. We need a more sophisticated approach 
in which the running status of a user depends not only on program length, 
but also on the nature of the resources employed since the beginning of 
computation. Thus, a mapping function must consider a user's position 

in both resource space and needs space for the previous few moments and 


the present moment. (See Potter, Appendix B.) 


Automatic Flowcharting .- Daniel U. Wilde 


The objective of this research ty development of a tool to aid advanced 
programmers in analyzing and debuggin digital computer programs 
written for general-purpose, single-address machines. Input to the 
analysis program is a BCD output-listing tape from the 7094 FAP 
assembler. Output from the analysis program will be a flow chart in 
symbolic form showing the effects of control flow and data flow. The 


following development work has been completed: 


A. DATA GATHERING 


On the first pass, the analysis program reads a BCD assembly tape and } 1 
produces the following tables as a function of the instruction location: j 
a. Entry and Exit transfer tables, 4. 
b. Active and Passive reference tables, \ 
ee OE | 
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c. Data and Storage pseudo-operation tables, 

d. A Symbol table, 

e. An Entry-Point or Starting-Location table, and 

f. A First-Instruction and Last-Instruction location table. 

B. DATA REDUCTION 
The second pass of the analysis program processes the tables constructed 
during the first pass. Examples of the processing are: 

a. The Entry and Exii tables are used to divide the program into 
"blocks''. A "block'' is a set of instructions between a transfer 
entry point and the next transfer exit point. Thus, a "'block!' is 
cornpletely processed if its first mmember is executed. 

b. The Active and Passive reference tables are used to construct 
the Constant and the Result tables. 

c. The Data and Storage tables are used, along with the Constant 
and Result tables, to determine if the input program changes 
or modifies any of its own instructions. 

C. MICRO DATA COMPRESSION 
The third pass of the analysis program processes the tables of the second 
pass. Examples of the processing are: 

a. The Block table is used to determine if all non-data blocks can 
be reached from either the starting location or the various 
entry locations. If a block cannot be reached, thenthe topology 
of the surrounding blocks is studied to determine what con- 
nection assumptions should be made. 

b. The connected Block table and the Active and Passive tables 
are then used to construct the Latest reference table as a 
function of program topology for each passive reference. 

hy D. MACRO DATA COMPRESSION 
\ a. Each active reference represents data changes that are made 
by the program being analyzed. By using the Latest table entry 
| ¥ 
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for each Passive table entry, a symbolic expression for each 
Active table entry can be constructed as a function of program 


topology. 


b. Cross reference information showing Latest values of used 


registers, ¢e.g., index registers, is also provided. 


The final flow-chart output can be set to show either a topology-only level 
or topology-and-data-processing level. Further refinement of these 
output formats is necessary, and this requires further consideration of 


possible data-flow representations. 
Visual Information Processing - Herbert M. Teager 


Major effort, under ONR Contract Nonr-1841(69) and JSEP Contract 
DA36-039-AMC-03200(E), has been devoted to the area of human sensory 
processing, primarily visual. The object of this effort has been three- 
fold. First, froma pragmatic viewpoint, to provide for handwritten and 
drawn symbolic input to digital systems with a latitude and scope approach- 
ing normal human practice. Second, to develop a testable theory of the 
visual processing that occurs in living'systems, since our theories with 
respect to pattern recognition are generally inadequate for almost any of 
the commonly encountered facets of human visual capability, such as 
visual memory, depth perception, face and object recognition, et al. 
Third, to understand the "what" and "how"! with respect to sensory pro- 
cessing, which is a fascinating and potentially fruitful field of scientific 
research in its own right. 

The work that was accomplished has been presented as a paper, entitled 
Multidimensional Visual Information Processing, to the New York Academy 
of Sciences. There is not the space available to duplicate the material of 
that paper here. However, the results with respect to development of a 


useful system are summarized below. 


It is possible to build a simple hardware system to accomplish real-time 
character recognition, when connected to a 'Teager table/Rand tablet", 


for a known individual's writing with the following constraints and 
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capabilities: Printed-writing script must be disconnected and at the 
general level of complexity of Roman character fonts, with no limitation 
On writing speed. The system is invariant to size, location, rotation 
skew, and all other affine transformations generated by oblique projec- 
tion. Within the limits of character complexity and delay-line storage 
Capacity, the system is unlimited with respect to character font. In 
practice, a font of no more than 256 different characters is optimal with 
respect to the hardware. These conclusions have been verified, with 
simulated hardware, by an extended period of experimentation on writing 
from a varied group of subjects, during which the very low error rate 
that occurred (less than 3%) could be traced almost exclusively to the use 


of symbols that were completely new to the system. 


Further work is being expended to build the hardware, as well as extend 
its capability to non-real-time recognition, unknown source, connected 
writing, anda character font at the level of complexity of Chinese and 


shorthand, 


With regard to the development of a tentative theory for visual processing 
of two-dimensional information, conclusions have been reached with 
respect to the form of measurement and processing which may well be 
Occurring in living systems if the known experimental, psychological, and 
artistic evidence of line, shape, and form incoding is considered. These 
measurements, in fact, form a basis for the processing system that has 
been developed. Evidence for such measurements and analogous ones in 
other sensory modalities are being sought, and consideration is being 
given to means of artificially accomplishing such processing with con- 


temporary processing and sensing hardware. 


One tentative theoretical conclusion with respect to visual processing is 
that optic nerve connections are far too few in number to be carrying in- 
formation to the brain in the assumed one-way, binary fashion. Obviously, 
analogous speculations are for future exploration in conjunction with 


appropriate life and medical scientists. 


At the operational level, a multiplexing system of hard-copy graphical out- 
puts and its attendent programming was turned over to the ESL group; and 
Eric Westerfeld developed a program to translate hand-drawn flow charts 


into machine-compilable MAD programs. (See Westerfeld, Appendix B) 
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CTSS Operation and Equipment - M. V. Solomita = 


i During this reporting period the IBM 7094 was scheduled for operation 
twenty four hours a day, seven days a week with the exception of a shut- 
down for the Christmas and New Years Holidays. The computer was also 
not in operation during the installation of the IBM 7289 channel and 7320A 
high speed drums in September, and the replacing of the IBM 1301 Disk 
storage file with the IBM 2302 Disk storage file in January. These changes 
have helped increase the speed and storage capacity of the Compatible 


Time-Sharing System, 
In June of 1965 the GE Model 635 began arriving and by June 15th the 
installation of the equipment began. We expect operation of the computer 


to start during the third quarter of 1965, 


The IBM 7094 Data Processing System configuration now consists of the 


following: 
Type Description 
" 7109 Control Processing Unit 
7110 Control Processing Unit 
oe 7302 Magnetic Core Storage (32, 768 words each) 
leo Console Control Unit 
3 7607 Data Channels 
3 7617 Data Channel Consoles 
2 ‘7909 Data Channels 
7606 Multiplexor 
7608 Power Converter 
7618 Power Control 
7631 File Control Unit 
2302 Disk Storage Unit (234 million characters) 
7320 Drum Storage Unit (1,118,400 characters) 
7750 Data Communication Channel 
7320A Drum Storage Units (1,118,400 characters each) 
12 729VI Magnetic Tape Units ‘ ~ 
716 Alphabetic Printer 
711 Punched Card Reader ' 
ical Card Punch if 
7289 Data Channel : 
j 
} 
| 
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Peripheral (Off-Line) Equipment: 


Type 
1401 Central Processing Unit (4,000 characters) 
1402 Card Read-Punch , 
1402 Printers (132 character/line each) 
; 4 729-V Magnetic Tape Units 


| 
| 
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Research on the Compatible Time-Sharing System 
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MULTICS'! Hardware System Design 
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Research on the Compatible Time-Sharing System - Fernando J. Corbaté 
SU AOI ALD LORIN CS ONAN PES] 


The system programming group has continued developing the Compatible 
Time-Sharing System (CTSS) for the IBM 7094 computer of Project MAC. 
This evolution has allowed the constant modification and testing of new 
ideas in a framework of man-machine interaction. In fact, the ability c/ 
the system to evolve with time has been a basic design goal. Hardware 
requirements in the past year have also caused further programming re- 
quirements. TwoIBM 7320A high-speed drums, attached to the MAC 


computer in August, have significantly improved response times. 


Thus, time-sharing at Project MAC has been: 1) a service facility which 
has allowed individual users to explore the consequences of man-machine 
interaction, and 2) a laboratory in which Systein programmers can develop 
the tools needed for effective man-machine interaction. The 7094 system 
is still in the process of being tuned, though it is incapable of being tuned 

to any high degree because the basic multi-programming ability and 

storage allocation techniques are very difficult. For this reason, a second- 
generation time-sharing system called MULTICS (for Multiplexed Informa- 
tion and Computing Service) will be implemented initially on a GE 645 asa 
development project with the Bell Telephone Laboratories and the General 


. 


Electric Computer Department. 


The Multics Project was started this year, and papers describing the 
system are to be presented at the Fall Joint Computer Conference, 1965. 
The major design decisions which have been made are in the area of:hard- 
ware. A separate report by E. L. Glaser covers hardware design. Soft- 


ware plans have been harder to firmly establish. 


Tuning ana improvement of the 7094 system has been done with two motives; 
there has been the desire to develop an accurate pilot model of possible 
future systems, and there has been great interest to learn from mistakes. 
On both counts it is felt that considerable success has been achieved during 
the past year. A detailed description of the 7094 system may be found in 
the Compatible Time-Sharing System, 2nd edition, M.I. T. Press, 1965. 
What follows is a brief review of the major elements of the evolvement of | 


the 7094 system. = 
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It has been recognized that the disk file system is a major and central 
part of a time-sharing system. Already CTSS has one of the most 
elaborate file systems developed. Nevertheless, it is also recognized 
that the first version of the file system is inadequate for a more highly- 
tuned time-sharing system. Some of the changes that have been improve- 
ments are the salvager and an improved input/output editor for placing or 
removing material from the disk file. These changes have improved the 
reliability and smoothness of system operation inasmuch ag the features 
involved occupy pivotal software roles and are expected to always function 


without mishap. 


However, of major importance has been the design of the second-generation 
file system. This major redesign is expected to correct four areas of 
difficulty. One, it should be possible to share files among many users, 
thereby allowing multiple reading and suitable interlocking upon writing. 
Two, multiprogramming should be possible such that more than one user 
can be using the file system simultaneously. Three, it should be possible 
to use magnetic tapes from within the time-sharing system. Four, it 
should be possible to incrementally dump only newly-created material, on 
a continuous basis, rather than have a massive dump periodically (e. gs, 
once a day) of the entire contents of the disk. This process is being 
developed under the program name DAEMON. Basic design was concluded 
in the fall of 1964, Implementation was begun and is now approaching the 
terminal stages. The implementation of this new file system is considered 
of prime importance, as it gives the programming staff an opportunity to 


shake-down the multitude of ideas which are involved in anticipation of 


.MULTICS. 


Other improvements to CTSS have been the introduction of a class of new 
editing programs for typing in programs, text, etc.” Of notable interest 
has been the development of the TYPSET and RUNOFF commands for 
English text and the corresponding ED and EDL commands for program 
text. These programs are based on a synthesis of many editing techniques 
by J. Saltzer and have mostly replaced the need for the previous editors 
such as MEMO, MODIFY, DITTO, INPUT, EDIT, and FILE, The new 
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editors are distinguished by a high degree of sophistication for man- 
machine interaction as well as context editing techniques which avoid 


fixed format requirements. 


In addition, nearly all commands of the system have had some corrections 
or improvements. Particularly notable are: inclusion of a new version of 
the MAD compiler; cleaning up and fixing up of the FAP command to allow 
use of more compacted input files; improvement of the LOAD commaiid to 
allow three different program libraries and thus speed up library search- 
ing; and reprogramming of FILE, RUNCOM, SPLIT, COMBIN, SNOBAL 
and other commands. Further new features have been added in the form 
of interconsole message facilities and the ability to ''attach" other consoles. 
In addition, new debugging techniques have been introduced, the most 
notable being the MADBUG program of R. Fabry for the FAPDBG-like 
debugging of MADBUG programs. The MADBUG program is considered 
only a prototype of future debugging techniques, as the present CTSS im- 
plementation forces it to be relatively clumsy. Nevertheless, the design 
goal, that the MAD language user need not know machine language, has 


largely been met. 


Finally, there have been numerous attempts to find solutions for the diffi- 
cult problem area of system documentation. There have been many at- 
tempts to compact disc space requirements and allow larger utilization of 
the available secondary storage Capacity. In particular, the ARCHIVE and 
CRUNCH commands have been developed. In addition, the LOG command 
has been developed to notify system users and allow user teams to keep 
posted all system changes of current interest. Moreover, the system 
programmers themselves have developed conventions and procedures in an 
attempt to systematize record keeping for a large system. In addition, 
facilities have been developed for semi-automatic "carry'' of programs back 
and forth between the two 7094 computers to allow smoother and easier 


maintenance of two slightly different CTSS systems. 


Of major importance in the orderly development of the new system has been 


a complete revision of the present CTSS Programmer's Guide to include 


information from many CTSS Bulletins and important MAC memos des- 


cribing changes in the features of the system. This revised manual will be 
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maintained on-line in the computer and is available to any user at any. hour 
of the day or night whether he be 20 feet or 2000 miles away from the 
Computer. In this Way, users may keep abreast of the latest system 
changes. In particular, a user is given the ability to read a table of con- 
tents of the changes listed in reverse chronological order, It is expected 
that the inanual, which can be maintained by the system programniing 
staff in as dynamic a state as the system itself, will be a Significant part 
of the answer to documentation difficulties of future remote-console 


systems. 


MULTICS Hardware System Design - Edward L. Glaser 


During the reporting period, the hardware design for the new Multiplexed 
Information and Computing Service (MULTICS) time-sharing system was 
formulated and frozen. This design consists of rather extensive modifi- 
cations to a General Electric 635 computer. The resulting computing 
System is to be known as the General Electric 645 computer, Extensive 
modifications were carried into a number of different areas. Hirst of all, 
the addressing logic was enhanced to incorporate the concepts of segmenta- 
tion and paging,* Additionally, Some new instructions were added to the 
instruction repertoire, and the interrupt logic was enhanced, making it 
Possible to interrupt an instruction part way through execution, and sub- 
Sequently continue with the execution of this instruction, These explicit 
modifications wil] enable the system to handle software envisioned for the 


new time ~Sharing system. 


This I/O control unit incorporated functions of both a Standard I/O contro] 
channel and a comnuiunications computer. Thus, the Same equivalent logic 
could be used to communicate with both standard computer Peripherals, 
such as magnetic tape, disc, and card equipment, and with varying com- 


munication lines from teletype speed to very high speed for displays, 


-_—_—— ss 


* See (1) Dennis, J, B, » "Program Structure in a Multi-Access Computer,!! 
Project MAC Technical Report MAC-TR-11; and (2) Dennis, J. B, and 
Glaser, E. L., "The Structure of On-Line Information Processing 
Systems,'! information Systern Sciences, Proceedings of the Second 
Congress, Spartan Books, Inc. 
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The generality of design is intended to make modular software in the new 
time-sharing System not only feasible but efficient. Further, specific 
features are incorporated into the design of the 1/0 equipment control unit 
to facilitate processing of many interrupts in a short period of time, which 
is necessary on any system connected to a large number of remote ter- 


minals. 


The third major area of design is that of 2 new high-speed drum controller 

to be used in conjunction with a large-capacity drum, The availability of 

this new drum controller with its associated drum will facilitate time-sharing 
for two reasons: first, the transfer rate between the drum and core 

memory is approximately fifty percent of memory speed; second, because 

of the way the drum controller operates, it is possible on a statistical 


basis to ignore any latency due to the rotary access to the drum. 


Details on all of these design areas will be available in papers to be 
presented at the Fall Joint Computer Conference, 1965, 
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Introduction - J. Francis Reintjes 


Several research groups of the Electronic Systems Laboratory are parti- 
cipating in Project MAC. Display systems research has been conducted 
in support of the Laboratory effort on computer-aided design and for the 
purpose of enhancing the output-display capabilities of the MAC time- 


computer system. 


The Computer Applicaticns Group continued its research in fundamentals 
of man-machine problem-solving, with emphasis upon the development of 
generalized computing techniques which are expected to be applicable to a 
wide variety of engineering design problems. Emphasis in this project is 


on the use of computers as an aid in mechanical design. 


Computer-aided electronic circuit design was also investigated under a 
grant from the National Aeronautics and Space Administration. The main 
efforts of the group were devoted to: analysis techniques for linear and 
nonlinear networks, logic design and synthesis through use of the MAC 
computer, design of threshold logic, and the development of a low-cost 


teletype-operated graphical display. 


The Project MAC time-sharing system has a'so been very useful for in- 
vestigating simulation of force-rebalance accelerometer Output data. The 
objective of this work is to obtain maximum Output information by using 
specialized data processing. In another aspect of the same research pro- 
gram, the group investigati ng computations peculiar to strapped-down navi- 


gation systems found the Project MAC PDP-6 to be highly useful. 


A total of 16 staff members and 12 students in the Electronic Systems 

. U 
Laboratory have benefited in their research through the use of Project 
MAC facilities. 


Dp lll 


eer 


72 ELECTRONIC SYSTEMS LABORATORY 


Display Systems Research - John E. Ward and Robert H. Stotz 


A. ESL DISPLAY CONSOLE 

The ESL Display Console was designed to satisfy needs of the M.I.T. 
Computer-Aided Design Project, and its construction Was supported by 
the Air Force Materials Laboratory, Wright-Patterson Air Force Base, 
under Contract AF-33(657)-10954. The unusual features of this display 
system, which have permitted its routine use in the Project MAC time- 
sharing system since January, 1964, are: its ability to time share the 
7094 memory through the Direct Data Channel; its incremental digital 

line generation, scaling, and rotation system; and its automatic -hardware 
pen-tracking system. During the past year, the System capabilities have 
been extended, with Project MAC support, by adding a second operator 
station (consisting of a slave display scope and associated manual inputs) 
which time-shares the display generation electronics, thus allowing two 
Operators to work simultaneously on different problems. Typical current 
usages, which are reported in other sections of this progress report, 
include: three-dimensional drawing programs, parametric surface studies, 
and portrayal of stress contours by the Computer-Aided Design Project; 
three-dimensional models of protein molecules by the Biology Department; 
modeling of the human vocal tract by the RLE Speech Group; design of 
highway interchanges by the Civil Engineering Department; and general 
graphical output by a number of other groups who use the display in an 
interactive way. Usage currently averages about 16 hours per day for each 


of the two console Stations. 


The M.I. T. Computation Center has recently decided to install a complete 
duplicate of the ESL Console on its time-shared 7094 to extend the avail- 
ability of display facilities to a larger segment of the M.I. T. community. 
The systems will be compatible and Project MAC display software can be 
used directly. Digital Equipment Corporation is building the electronic 
portions of this second system from M.I.T. drawings, and the manual in- 
put devices are being fabricated in the Electronic Systems Laboratory. 


This second console, which is supported by an NSF Grant, is scheduled for 


installation in November, 1965. 
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B. STORAGE-TUBE DISPLAYS 

Experiments with driving ''one-shot"' storage tube displays from analog 
deflection voltages of the ESL Display Console have been quite successful, 
although the quality of the present storage tubes leaves something to be 
desired. Two such displays, connected to the ESL console by three 
coaxial cables, have been in experimental use during the past few months. 
Programming and control circuit modifications to permit the addressing 
of up to ten remote storage scopes for transmission of ''one shot" pictures 
have been accomplished. Figure 3 shows typical storage. displays. These 
require only about one half second of console time to create, but are re- 
tained for periods up to one hour, unless erased manually or by the com- 
puter. Experiments are continuing with the use of analog data sets on 
three telephone circuits (for X and Y deflection voltages and intensification 
pulses) for remote operation of storage scopes in this mode at greater 
distances. Noise and linearity characteristics of the telephone circuits 
appear satisfactory, though the difference in time delays on the three 


circuits is proving troublesome. 


C, DISPLAYS FOR THE GE 645 COMPUTER 

Study of future display requirements for the Project MAC GE 645 Computer 
has resulted in identification of two goals: 1) a modest number (10-15) of 
remote sophisticated displays, with capabilities similar to the present ESL 
Console; and 2) a large number (perhaps 500) of low-cost, limited- 


capability displays which can be located at every MAC terminal. 


Planning for low-cost remote display scopes currently centers around a 
digitally-controlled unit (rather than the analog transmission system dis- 
cussed in B, previously), with line and character generation being performed 
locally. A set of desirable characteristics for a unit which could operate 
from standard 2000-bit-per-second Dataphones has been prepared and is 
being discussed with the manufacturers of display equipment. Also, 
research work is being conducted on circuits for use in such a device. One 
master's thesis project by Mr. T. B. Cheek concerns new approaches to 
character generation, with the objective of exploiting low-speed requirements 
to achieve very low cost. A low-cost, hybrid digital/analog, line-generation 
system is also being designed to operate from a Dataphone. This approach 


combines a binary rate multiplier with an analog integrator, thus 
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Figure 3. Sample Storage Tube Displays 
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eliminating the need for digital-to-analog converters. A prototype unit, 
With parts costing about $1,000, should he ready for initial tests by 


September, 1965. 


Operating more sophisticated remote displays (like the ESL Console) 
directly from the central computer (in the manner of the present ESL 
Console) will present definite communications problems, unless a flexible, 
low-cost, multi-megacycle, two-way link becomes available. Also, the 
load on the central computer must be considered. Our experience indicates 
that the best approach is to use a small satellite computer with each dis- 
play (or group of displays) to buffer display data transmitted over medium- 
speed telephone circuits (4 to 10x 107 bits per second). The satellite 
computer would also perform AH atin display processing, suchas pen 
tracking, rotation calculations, and interrupt analysis. One such com- 
puter-buffered display (the DEC 338) will be available at Project MAC 
during the coming year for experimental evaluation and study of pro- 
gramming problems. (Note that the small computer program can be 


considered part of the supervisor prograim of the central machine.) 


D. TV SCAN CONVERSION 

The display group has completed the design for an image-maintaining 

display system based on a TV scan converter for the U. S. Naval Underwater 
Ordnance Station under Contract N140(122)-761-18B., Coniputer-generated 
pictures are drawn on the storage surface of ¢ scan-conversion tube by 
slaving its deflection inputs to an existing DEC Type 340 display system. 
Then, by switching the deflection ir puts toa TV raster, the stored I 
picture may be viewed for periods up to five minutes on standard TV 
monitors. A 1203-line scan is used to provide high resolution. Selective 
erasure or writing of up to three characters will be possible during each 

TV vertical retrace time, so that stored pictures may be altered and updated 
without requiring a complete regeneration by the computer. This system 
will be experimentally tested on Project MAC's PDP-6 computer before 


being delivered to the Navy later in the year. 
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E. IMPROVED DISPLAY TECHNOLOGY 

A number of other aspects of display technology are being studied as 
thesis projects. Mr. W. D. Stratton is working on a new high-speed 
pen-tracking system, based on the amplitude comparison of light-pen 
responses to displayed points ina tracking pattern, similar to radar 
conical-tracking schemes. By replacing the present digital point-counting 
schemes, it is hoped that pen-tracking time can be reduced by a factor of 
perhaps 20 to 1. The new electrostatic "beam pen" being developed by the 
group is being used in this work, although the technique could also be 


applied to light pens. 


The possibilities of infrared (IR) data links for remote Operation of dis- 
plays is being investigated in a thesis project by Mr. George Ling. We 
have obtained an experimental IR link, on loan from the Lincoln Laboratory, 
that is capable of 10? bits per second over a two-mile range, and have 
evaluated its performance in transmitting digitally-coded messages. The 
link is now being set up to transmit display information from the ESL 
Console in Technology Square to a storage scope in the Electronic Systems 


Laboratory, a distance of about 600 yards. 


As part of the continuing program to improve the ESL Display Console, Mr, 
Endre Guttman completed his thesis study of a digital differential-analyzer 
(DDA) rotation matrix to replace the present binary rate-multiplier (BRM) 
rotation matrix. In the BRM matrix, shown in Figure 4, the 10-bit signed 
numbers iv is etc., are set by the computer and are normally the vector 
cosines specifying rotation between the x, y, z and h, v, d coordinate sys- 
tems (depth coordinate, d, is not computed). Scaling of these numbers 
controls the size of the displayed image. The boxes labeled "increment 
logic'’ combine and buffer the BRM outputs, releasing them to the scope de- 
flection registers in binary weights of 1, 2, 4, or 8 steps, which the regis- 
ters (up-down counters) can accept. 

Although the BRM matrix is generally satisfactory, it has significant 
round-off errors which create picture distortions. Also, the buffering 


action of the increment logic adds some roughness to displayed lines. 


SR entre tog tet 


“e 


pS a 


— 


a a nr a te ae pe ee re ee an a 


ELECTRONIC SYSTEMS LABORATOR Y 


78 


XTTjeIN uotjeqoy yqqd pasodoig ‘¢ aInsI gy 


| GALLINO 91907 NOIS :3LON | 
SIXV-4 OL 


TVOILNSGI 


2V 


SuY31S193uy y 
NOIL937434a 7 HOLVYSN39 
3d090S a} 3NIn 
O1 dav N i WOU4 
OL Vv SLNdNI 
SLNdLNO W Nivel 
= 3S 1nd 
uy 
XV. 
uy. MOTEUSAO ea 
ee ee ee | —— i 
YSLNdWOD Ag LAS Sy aNv'sf ‘sy 


eet 


nn ss 


ELECTRONIC S¥STEMS LABORATORY 79 


Figure 5 shows the Proposed matrix, which would have two identical three- 
input parallel DDA's, one for h, and one for v. The advantage of the DDA 
Over the BRM is that it has a remainder register to 'remember' round-off 
in a particular step of the calculation, and apply it to the subsequent step. 
Also, the h and v outputs would be obtained directly, eliminating part of 
the function of the increment logic. The increment logic could then be 
eliminated completely by converting the first few stages of the deflection 


registers from counters to adders. 


A DDA for this purpose, however, must be able to add four, ten-bit, signed 
numbers ina total elapsed time of only 1.5 microsecond. It was shown that 
multi-input threshold-logic elements offer the only feasible way of meeting 
such a severe speed requirement at reasonable cost. A potentially suitable 
six-input threshold-adder stage (four register inputs, plus two carries) 
requiring only three transistors was designed and breadboarded. It was 
found to have a transition time of 80 nanoseconds; fast enough for this 
purpose. Further development and testing of threshold circuits will be 


necessary, however, before the DDA rotation matrix can be constructed, 


EF. FUTURE PLANS 

All reported activities will continue during the coming year, with particular 
emphasis on developing the required display technology for future graphical 
man-machine communication and studying the hardware/software problems 
in interfacing large numbers of displays with the next generation time- 


sharing computers, 
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Computer-Aided Design - Douglas T. Ross and Clarence G. Feldmann 


The Computer Applications Group, along with faculty and staff of the 
Mechanical Engineering Department, has continued its Computer-Aided 
Design work, with emphasis on problems of mechanical design. This 
program is being sponsored by the Air Force Materials Laboratory, 
Wright-Patterson Air Force Base, under Contract AF-33(657)-10954. The 
work is an evolution of earlier work, conducted at ESL, to develop the 
Automatically Programmed Tool (APT) System now being used by industry 
to program numerically-controlled machine tools. The approach being 
taken to computer-aided design is sufficiently general tu make it applicable 


to a wide variety of design problems in engineering and science. 


The current research effort has four phases. The first phase involves 
research in fundamentals of man-machine problem-solving, the structure 
of language and language processing, the modeling of problems, and the 
development of generalized computing techniques. The second phase 
concerns application of the theoretical foundations resulting from this 
research to the design and construction of a generalized man-machine 
problem-solving system, including development of a powerful language for 
Writing, Operating, and debugging programs, and facilities for efficient 
construction of specialized problem-oriented systems. The third phase 
is actual development of such systems for specific design areas, and the 
remaining phase is adaptation and application of special systems by 


designers to meet their individual needs. 


In the realm of generalized languages and systems, the AED-0 system 
(standing for Automated Engineering Design) was the first system completed 
by the group. It has been in Operation for over a year and is being employed 
extensively by members of the Design Group and by other groups at Project 
MAC. More recently, through the efforts of Mr. Feldmann and his associates, 
AED-0 has been made available as a batch-processing system and distributed 
to interested industrial users. Mr. Ross presented a series of lectures on 

the salient features of AED-0 language, and his notes now constitute the 


basic AED-0 Programmers! Guide. 
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During the year, many.substantial additions were made to AED-0 in order 
to increase its power for writing advanced programming systems. The 
language is based on Algol-60, with many additions and modifications to 
permit construction and manipulation of complex data structures, includ- 
ing the ability to reference subparts of computer words for efficient use of 
storage. The system also includes an elaborate Macro Preprocessor which 
permits extensive manipulation of the input string of characters before 
compilation, thus providing the user with greatly enhanced expressive 
capabilities. A further feature is that the Préprocessor includes facilities 
for regenerating AED-0 language programs in a meaningful format, which 
exhibits the logical structure of the program. Combined with flexible 
comment" and "remark"! features, this capability makes AED-0 programs 
almost self-documenting, minimizing the need for elaborate flow diagrams 
to describe programs. A source-language debugging system has been 
written, so that the user may trouble-shoot programs at the AED-0 lan- 
guage level, with little concern for the actual machine code generated by the 
compiler. AED-0 development is now complete, and primary attention is 
now focused on its use as a tool to achieve even more sophisticated gen- 


eralized programming and compilation systems. 


The immediate successor to AED-0 will be the AED-1 System now under 
development. A salient feature of AED-] is that it will be highly machine- 
independent. That is, it is being designed so that both the language and 
techniques used to compile machine-code programs will be compatible 

with virtually any type of large-scale computer. Several valuable charac- 
teristics of AED-0 will be contained in its successor, as well as new 
features to give added flexibility and capability. The AED-1 System will 

be able to process AED-0 as well as AED-] language programs. Several 
building-block (subsystem) programs are being written to achieve this goal, 
including a subsystem (AED Jr.) which permits the user to define and check 


out descriptions of arbitrary programming languages. 


With respect to specialized, problem-oriented programming systems, the 
group's first major system in this category is now being prepared. Called 
CADET (Computer-Aided Design Experimental Translator), it permits 
users to carry out computer-aided design applications using verbal or 


graphical language or both. That is, its purpose is to enable construction 
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and manipulation of graphical and analytic structures. CADET will also 
permit analytic functions and physical properties to be associated with 
graphical representations for things. Since programs defining behavioral 
properties of graphical elements may take arbitrary form, they are written 


in standard AED language. 


To assist the initial development,of CADET, which will be a general- 
purpose problem-oriented system, several specific design problems are 
being explored through its use. One of these is concerned with the general 
problem of three-dimensional shape description, and-will use the generalized 
parametric surface patches developed by Professor Coons of the Mechani- 
cal #ngineering Department, as well as the three-dimensional pseudo-pen 
prograim written by Mr. Polansky, and the routines prepared by Mr. Lang 


for ESL's Graphical Display Console. Another application is in the area 


of linear and nonlinear electronic circuits, described in the next section. — 


A further aspect of the entire program in Computer-Aided Design is the 
work being performed on cathode-ray-tube displays for graphical mani- 
pulation of design information and symbols. This work was described in 


the previous section. 


An important goal of the Computer-Aided Design Project is the dissemina- 
tion of information generated by the group to interested users outside M.I. T. 
It was our pleasure to have seven system programmers from industry as 
guests of ESL and Project MAC for the past year. Our guests have parti- 
cipated in developing the programming systems described previously, and 
they plan to apply and further expand these basic concepts upon return to 
their respective organizations. This program of cooperative research and 
development with visiting staff from industry will continue as an important 


aspect of ESL's work in computer-aided design. 
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Computer-Aided Electronic Circuit Design - Michael L. Dertouzos and 


J. Francis Reintjes 


A. SIMULATION OF ELECTRICAL NETWORKS 

A method is being investigated for on-line simulation of electrical net- 
works. This method treats networks consisting of linear and nonlinear 
resistors, as well as inductors and capacitors. A given network can be 
excited with voltage and current sources which may be arbitrary functions 
of time. The voltage response across any node in the network may then 


be called for as output. 


The computer model for the network consists of lists of storage registers 
representing branches connected to other lists representing nodes. The 
method for solution at each instant of time is as follows: The nodes are 
assigned fixed node potentials. The voltage across each branch, and hence 
the current through it, is determined by the two nodes to which it is con- 
nected. Equilibrium requires the sum of the branch currents at each node 
to be zero. If the sum is not zero, new node potentials closer to the 
equilibrium values are computed by perturbing the old ones in a prescribed 
manner. This process is repeated until the equilibrium values are obtained 


within a specified tolerance. 


Thus far, a computer program for use on the Project MAC time-sharing 
system has been written to apply the foregoing method to linear RLC net- 
works excited with step and sinusoidal current and voltage sources. The 
program computes an accurate solution to any given linear network. Com- 
puter time for solution of networks containing up to about 8 nodes and 12 
branches is on the order of a few seconds. An example of a circuit 


analyzed by the foregoing program is shown in Figure 6. 


B, DISPLAY PROGRAM FOR ELECTRONIC CIRCUITS 

A computer program has been written which allows an operator to quickly 

and easily portray planar electronic circuits on the Electronic Systems 

Laboratory Console Display at Project MAC. Circuit elements are formed 

by moving a tracking light pen to the desired position on the screen and j 
actuating an ''element button.'' Circuit characteristics are described ina 


list structure with pointers. Used in conjunction with the analysis program | 
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Figure 6, Response of a Twin-T Filter as Computed and 
{ Plotted by the Circuit-Analysis Program 
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described previously, circuits with multiple elements can be displayed and 


their response to various types of excitation readily determined, 


C. ON-LINE SWITCHING-FUNCTION SYNTHESIS 

The central aim of logical design is the synthesis of any given switching 
function, in terms of given sets of elementary building blocks, for the 
optimization of some performance index in the presence of constraints. 
Although the present state-of-the-art yields algorithmic methods for the 
solution of certain Specific instances of the problem (such as minimization 
of building -block inputs with a two-level, AND-OR realization), no fully- 


algorithmic method exists for the solution of the more general problem, 


Solutions to the general problem are being developed through use of an on- 
line process, where the machine accomplishes those computational tasks 
which can be algorithmically specified, and the user provides those de- 
cisions which he is better qualified to make. The machine portion of the 
system is based on a set of heuristic procedures which, subject to certain 
conditions, guarantee convergence of the process. Results obtained to 
date indicate that this method of synthesis yields more economical struc- 
tures than those spine successive local-optimization Procedures, and does 
not depend on impractical (and usually impossible) exhaustive searches 
through all possible solutions. It is expected that a report on this work 


will be issued during the Fall, 1965. : “fi 


D. THRESHOLD ELEMENT REALIZABILITY 

Previous work has shown that the problem of determining whether an 
arbitrary Boolean function of N variables is realizable with a single threshold 
element can be reduced to the problem of minimizing a quantity which is a 
function of N+1 variables, Several properties of this quantity are being 
derived and studied. It can be shown, for example, that the quantity may 

be interpreted geometrically as a structure of intersecting N+1 dimen- 

sional hyperplanes. Each of these hyperplanes corresponds to some N- 
variable Boolean function, and if the given Boolean function is realizable, 

its associated hyperplane will be "flat" with zero height. The hyperplanes j 
are so situated that they form the boundary of a convex body and hence the 
quantity contains no "pockets"! of local minima. This property enables a 


computer to employ hill-descending techniques in N+] Space which seek a 
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local minimum point. Since, by virtue of the convex nature of the quantity, 
any pOint satisfying local minimum conditions must also satisfy absolute 
minimum conditions, hill-descending techniques will lead to the required 


solution. 


Several hill-descending techniques for performing a minimization process 
are being studied and are yielding dependable results. The question of the 


complete generality of these methods, however, is still an open one. 


E. A LOW-COST PLOTTING BOARD 

A hard-copy graphical-display device that can be operated from low- 
speed, teletype data sources such as those found in time-shared computer 
installations, is being investigated. A novel feature of this device, shown 
in Figure 7, is the conversion of incremental binary data into graphical 
form. This conversion is primarily accomplished by the use of tuned 
dynamic networks at each axis. Inclusion of mechanical networks, use of 
certain specific torquers for actuators, and over-all system organization 
were motivated by a desire to keep cost at an absolute minimum. This 
objective has not adversely influenced performance or speed of operation. 
The latter is fundamentally determined by low teietype information rates 
(100 bits per second) and by upper bounds on the torque-to-inertia ratios 
of reasonably-sized actuators. Theoretical extensions of the foregoing 
basic principle have shown that it will be possible to implement curvature 
control, data smoothing, high-bandwidth plotting and other desirable 


properties. 


F, SIMULATION OF NONLINEAR CIRCUITS 

A different attack on computer simulation of electrical networks, from that 
described in Section A., previously, is being taken in this investigation, 
Primary emphasis is on nonlinedr-network simulation, and network response 
is obtained through succeSsive solutions and modifications of a matrix which 
is a mathematical description of the network being simulated. A computer 
program is now being written in the AED-0 language to establish the net- 
work equations and solve them. Later work will include programming of 
the ESL Display Console for graphically portraying network diagrams and 
responses to various input signals. The system will be able to handle time- 
dependent sources and RLC elements whose characteristics are piecewise 


linear and can be either monotonic or nonmonotonic. 
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Figure 7, A Low-Cost Teletype-Operated Plotting Board 
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Accelerometer System Studies - Alfred K. Susskind and Donald R. Haring 


A novel approach to obtaining digital data from force-rebalance accelero- 
meters is under investigation for the U.S. Air Force under Contract 
AF-33(657)-11311. Because the proposed system contains both sampling 
and quantization, it cannot be treated adequately by analytical techniques. 
Therefore, experimentation is required, which can be economically done 
only by means of simulation. The Project MAC Time-Sharing System has 
been found particularly well suited for this purpose, because of terminal 


equipment such as the ESL Display Console. aa 


By interplay between theory and simulation, good progress has been made 
in understanding and evaluating the proposed scheme. The observed 
system performance is encouraging and the scheme will therefore be 
investigated in greater depth, gradually removing more and more of the 
idealizations now being made and also observing the effect of noise. 
Eventually, a rather complex dynamic system will be simulated, and it 

is expected that sufficient information will be obtained to act as sound guide 


lines in the design of an actual accelerometer system. 
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Simulation of Strapped-Down Naviation Systems - F., B. Hills and 
J. A. Parisot 


In a research program for the Air Force under Contract AF-33(657)-11311, 
we are studying the computations peculiar to "strapped down" navigation 
systems, in which inertial sensors (gyros and accelerometers) are 

rigidly mounted to the airframe rather than being mounted on a stabilized 
platform. The accelerometer outputs must therefore be transformed from 
body axes to the coordinates in which the navigation computations are 
made. We are using the direction cosine transformation. Because of 
vehicle rotations, the direction cosines are not fixed and hence must also 
be computed. This is done by numerically solving nine simultaneous 


differential equations in which gyro data are the independent variables. 


Analytical studies made by the group have shown that the major sources 

of error, insofar as the computations are concerned, are due to the un- 
certainties that arise in vehicle position and attitude because the input data 
must be sampled and quantized. The uncertainties arise primarily be- 
cause acceleration and rotation are not communtative and the components 

of rotation are not commutative. Hence, not only is the intersample 
behavior of the input data important, but even more important is the 

relative intersample behavior between the various input variables. Some 
analytical studies into the size of the uncertainty have been made; however, 
the investigation becomes intractable for all but simple flight paths. There- 


fore, analysis must be augmented by simulation. 


Project MAC's PDP-6 was chosen for these simulation studies because of 
its ability to perform double-precision arithmetic at high speed (46-bit 
accuracy is required), and its suitability for on-line operational control 
(e.g., possibility of connecting external equipment, such as a stick to 
control the flight of a simulated vehicle). Also, macro definitions allowed 
by the PDP-6 assembly program make it possible to set up desired tests 


with a minimum of debugging. 


A diagram of the simulation is shown in Figure 8. In the box labeled vehicle 
simulation, the input data, acceleration (a) and rotation rate (w), are gen- 


erated. At present, these quantities are generated in accordance with the 
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geometry of simple flight paths. In the future, a simple dynamic model 
of the vehicle with manual control may be used to allow a search for 


maximum uncertainty. ‘ 


Two computations of the coordinate conversion computation must be made: 
one for the system under test, the other for a reference computation which 
is considered to be the true solution: The reference computation, like 
that of the system computation, requires sampling and quantizing of the 
input data. Analysis has shown that all numerical formulas compute those 
solutions to the direction cosines that would result if the components of a 
were proportional between samples. Therefore, to make the reference 
computation more accurate than system computation, the samples for the 
reference must be taken more often than those taken for the system. This 
produces a reference solution based on information in the input data not 
available to the system computation, and hence a meaningful uncertainty 
can be computed. The high sample rate of the reference computation 
requires high-speed computation to obtain a close to real-time simulation 
and long word lengths to limit the accumulation of round-off errors to 
acceptable values. These two aspects are the major sources of difficulty 


in programming the simulation. 


The block entitled accelerometer and gyro simulation performs the sampling 
and quantizing done by the inertial sensors. Observe that the outputs are 
not acceleration and rotation rate, but incremental changes in velocity and 


angle, respectively. 


A set of macros have been written, so that tests can be set up and para- 
meters changed with a minimum of effort. Where several tests are to be 
run on the same input data (e.g., changes in word length), provision has 
been made to store pertinent results of the reference computation. Hence, 
the time-consuming reference computation need be performed only once. 
Also, an output routine has been written which produces pages of columnar 
data, each column with a descriptive title. Numbers can be printed in 


decimal or octal radices, and in fixed-point or floating-point forms. 
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Technical Information Project (TIP) 
Prccess Control: Serials and Journals 
Search Procedures 
Measures of Relatedness 
Statistics of Words in Titles 
Educational Use of TIP 


Use of TIP to Update a Data Compilation 
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Technical Information Project (TIP) - Myer M. Kessler, William D. Mathews 


The M.1I.T. Libraries’ Technical Information Project is attempting to 
utilize the time-sharing facilities of Project MAC for literature search 
and retrieval. The system is a prototype, based on twenty-five physics 
journals chosen from a reference matrix started by the Physical Review. 


This literature extends back five years, the average. 


For each journal article, we record: location of the article (journal, 
volume, page), title, authors, institutional affiliation of the authors, cita- 
tions (journal, volume, page); location of the article in Physics Abstracts, 


and subject-index information if available from a published source. 


This information is edited, compressed, and.put on the MAC system as ¢ 
data files. The user may perform search and retrieval operations by 

means of a special language that consists of several dozen words to control 

search programs. The simplest form of engagement may be initiated by 


specifying three commands: 


SEARCH - defines the range of literature to be searched, 


FIND - states the items to be found, 
OUTPUT - defines the nature and content of the computer 
output. 


The SEARCH command may take several forms: 


SEARCH ALL - will search the entire literature 


in store, 


SEARCH ALL NEW will search the last volume of 

each journal, 

SEARCH PHYREV ALL - wil search everything in store of 
a given journal (in this case the 
Physical Review), 


SEARCH PHYREV V. 120 to V- 135 - only volumes 120 through 


135, 
SEARCH PHYREV V. 135 - only volume 135. : 
The program will search as specified above and detect any item described J i: 


under the FIND command, which itself has a variety of possibilities. 
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FIND AUTHOR SMITH - will find all papers that include 
Smith among its authors. 

FIND TITLE CRYOGENICS - will find all papers that contain 
the word "'cryogenics" in the 


title. 


One may similarly define a location, 


FIND LOCATION IOWA STATE - 
or a citation, ‘ 


FIND CITATION 1 131 1165 - to find all papers that cite 
the article in Physical 
Review (code = 1), volume 
131, page 1165. The find 
instruction includes the 
logical manipulations of "and/ 


or/but not." 


Several other sophisticated FIND commands are available, the most in- 


teresting of which is: 
FIND SHARE B journal, volume, page. 


In this command we are looking for papers in the search range that share 
a ditation item with a given paper. This mode of literature search, called 
bibliographic coupling, is a powerful tool for information retrieval. The 
final command, OUTPUT, selects one or more of the options available as 
output. 

j 


/ 


The literature and program facilities of TIP are available to all MAC system 
users, and a monitor program has been written to collect statistics on the 
type and extent of MAC-TIP utilization. The most extensive use of MAC- 
TIP to date has been in connection with a book on Plasma Physics, written 


x by Prof. Sanborn Brown, as noted later in this section. 
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Process Control: Serials and Journals - Patricia M. Sheehan, and 


Marianne C. Henneman 


As libraries grow in size, the problems associated with maintaining con- 
trols grow in complexity. Difficulties are most acute in university 
libraries with scientific commitments, where periodicals represent a 


large proportion of the holdings. 


During 1965, TIP has made a systems study of the traditional manual and 
punch-card methods currently used to control serials and journals in the 
M.1I.T. Libraries. A feasibility report has proposed using serials and 


journals as a microcosm for a computer approach capable of: 


1) extending the features of the MACTIP information 
retrieval system to M.I.T. and other regional users, 

2) offering alternates to supplement the present card 
catalogues, 

3) providing a two-way communication with the computer, 
interjecting reports and trigg*ring actions, and 

4) performing the usual data-processing functions for 


administration, accounting, and statistics. 


The proposed system will be based on a permanent disk file and a series 
of programs combining remote console and direct input means for updating 
and retrieving data. Initial conversion of punch-card data, anda 

limited retrieval program have been completed. Additional programs are 


now being written to purify and expand ‘files to meet defined objectives. 


As one of the largest files in the time-shared system, the library files 
need careful program logic to minimize the effects of continually expanding 
and revising records and items of variable length. A linking technique 
permits simultaneous use of sorted and unsorted data. Additionally, a 
multi-level sorter tries to reconcile file data having key breaks, aftera 
comparatively small number of computer characters, with data having a 
break delayed until the end of a lengthy field. Finally, sorting and re- 
trieval programs expand a variable number of variable-length fields within 
a record and align each field with the corresponding data in other records. 


The pragmatic solution at this point is to have the programs determine 
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tentative maximums, and also monitor and adjust the figures as necessary, 
The extent to which users permit current sources of information to be re- 
placed by computer -produced catalogues and reports will be governed, in 
part, by aesthetics, For this reason, we hope to ultimately have an off- 
line printing arrangement, perhaps using Photon or Similar equipment, 


that will be more versatile than standard computer printers, 


Search Procedures - Evan L. Ivie 
sr ocedures 


The Technical Information Project (TIP) is based on a measure of the 
relatedness between documents that is a correlation coefficient based 

on information theory. This measure was described in the Project MAC 
Progress Report for 1964 (AD-465-088). A quantitative value is assigned 
to the ''closeness" of every pair of documents in a library, based on th: 


way the library is used, 


A program set has been devised which uses this measure to convert a 
request, consisting of a set of interesting Papers and a set of uninteresting 
Papers, into an answer set of related papers. This Procedure now allows 
a user to be very specific, and receive only a few documents in an answer 
set, by indicating no interest in marginally-related documents; or he may 
be very general, and obtain a correspondingly larger answer set, by in- 
dicating interest in almost everything that is presented to him fora 
decision. To ensure that the procedure wil] always converge to a set of 
related documents, it is necessary that the total internal correlation of the 


set being formed always increase, 


a restricted subset of English. Some examples of this language are: 


"Which papers are related to papers by J. R. Jones but not 
Papers using acoustic." 
"Print for decision the titles of articles related to articles 


citing Phys. Rev., Vol. 136, page 22," 


a 12-state automaton With 8-input word types. The file structure used is 
"inverted" in nature. The parameters (words, citations, authors, etc.) 


reside in ordered lists having connectives to other parameters, 
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Measures of Relatedness - John P. Casey 


The volume of journal literature is increasing rapidly. An understanding 
of the structure of this literature may’ provide some interesting informa- 
tion about an expanding communication network, and may help in the 
design of systems to manipulate and order these articles. It is not clear 
what type of structure will fit journal literature, so tentative beginnings 
have been made in three different directions. 

First, a collection of articles may be considered as a graph, with the 
articles as points on the graph. An arc connects two articles if they both 
cite the same article. Questions can now be asked how tightly a subgraph 
is connected together by its arcs, or how closely connected an article is 
toa subgraph. These questions are closely related to the difficulties of 
designing a mechanical system to find a group of articles all related to 


a given reference. 


Secondly, this same network may be considered as a group of nodes con- 
nected by bars, with the subject discussed in a given article diffusing like 
heat along the bars. Coefficients of perineability may be assigned to the 
bars, and rather complex formulas derived for the amount of "heat'! ex- 
changed between an article and others scattered through various journals, 


to determine which are most closely related to it. 


Thirdly, an analysis using the methods of Rashevsky's mathematical 
biology has been started. From a number of simplifying assumptions 
about the distribution of useful ideas in articles and related factors, a 
simple model may be constructed which predicts the likelihood of an 
article being cited and the time span until an article is no longer men- 
tioned in the literature. This may then be checked against statistical 
data available through the TIP system. It is hoped that a model which 
makes reasonably accurate predictions will clarify the general charac- 


teristics desirable in a literature-searching system, 


We are investigating, for a short distance at least, each of these ap- 
proaches. Thus far, there have been no results that would show the 


feasibility, or lack of such, for any of these, 
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Statistics of Words in Titles - Elizabeth A. Dole 
ee ords in titles 


Present studies attempting to describe several characteristics of TIP's 
vocabulary are based on word-count lists compiled from six years of 
Physical Review, Volumes 113-132 (1958-1963), and Volumes 133-136 
(1954) Series A, 


Several criteria are used to assessed total vocabulary, gross word count, 
and average words per title. The basis of each criterion is an exact 
definition by formula or by known arbitrary rule. The first criteria 
developed are those describing what is a "word," and the groups defined 
are successively rejected from the Original unrefined vocabulary. The 


following steps are planned: 


1. A "word" is any string of characters between two Spaces. 


Sample Results: 


; number of individual words = 4,826. 
individual words times number of uses = 56,413. 
number of articles in sample = 7,583, 
average number of words per title = 7.43 


2, Refinement A. A word is any string of characters between 
two spaces that does not start with: 
A. anumber 
B. + or -or= 


C. any punctuation mark 


3. Refinement B. (in progress) A word is any string of 
characters between two Spaces not ruled out by step number two, above, 
that 1s not one of an arbitrarily-s ecified list of words selected for their 
inherent lack of information context, such as connectives, 

In this step, the list of words is selected and tabulated. The 
impact of removing these items from the word list CROSS (A) will be 
calculated to give a new GROSS (B). 


4+. Refinement C. A word is any string of characters be- 
le 


tween two spaces not ruled out in steps two and three, above, that bears 


no basic Similarity to any other word, 
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In this case, words with notably similar bases or roots are 
grouped, and these counts are compounded in the total figures. The 
result should be a reduction in the number of individual words, but not in 
the gross number used for this step. 
5. Refinement D. Remove any remaining combination of 
letters and numbers that is not: 
a) an obvious symbol for a chemical compound, 
b) an obvious isotope symbol. 
6. Refinement E. Remove all chemical compound symbols. 
7. Refinement F. Remove all isotope symbols. 
8. Refinement G. Remove all unclum ed words occurrin 
Scie epee Words occurring 
five times or less. Evaluatethe remaining vocabulary. List all words 
a, 
so removed, otal 
SS Q 


9. Comparison I. Compare a portion (or more) of the results 
above with those gained from an exactly similar list and count from TIP's 


present full file of 25 journals. 


10. Comparison II. Compare search products from equal 
periods of TIP's library with the library of physics references available 
through equivalent editions of Physics Abstracts, and through the Physics 


Abstracts Subject Index Headings List effective at the same time. 


In summary, these vocabulary and retrieval studies should show the 
breakdown, from an information Standpoint, of the vocabulary of title words 
of physics articles as represented in the data base. The studies are dis- 
tinguished by clearly defined stages. These stages are defined explicitly 
enough that, once accepted, they could easily be expressed in logical pro- 
gramming for computer processing. The studies so outlined will afford 
needed information to designers of information files, particularly of in- 


verted title word files, in the field of scientific literature, 
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Educational Use of TIP - Myer M. Kessler, Sanborn C. Brown, and 


Irma Y. Johnson 


A freshman seminar based on TIP is being offered in the fall of 1965. 

This seminar is conducted with the cooperation of Prof. Sanborn C. Brown, 
of the Physics Department, and Mrs. Irma Johnson, Reference Librarian 
at the M.I.T. Libraries. 


Students will investigate the literature of plasma physics using TIP and 
compare it with the results of a search based on traditional library 
methods. The course is experimental and is being observed from two 


points of view: 


1, Can beginning students more profitably study a field of 
physics from current literature references ona computer, instead of 


traditional text books? 


2. Can beginning students compile useful bibliographies and 


compilations of physics literature from a computer-based system? 
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Use of TIP ta Update a Data Compilation - Sanborn C. Brown ren 


I have been investigating the possibility of writing a report which will 

demonstrate using a computer to update a compilation of scientific data. 
Six years agoI published a reference book,” and the present revision is 
based on bibliographic material contained in the physics literature pro- 
grammed in the M.I.T. Library's Technical Information Project (TIP). 


This program is available to all Project MAC users. 


Using this system, the attempt is to create material in an open-ended 
form, so that anyone with access to the computer program can search the 
literature for material which will appear after the text of my book is 
printed. As an aid for doing this, references to individual displays of 
data are given in computer-language form, rather than the usual biblio- 


graphic form. 


To illustrate the method of keeping up-to-date, let me give an example. 
Suppose that the reader is interested in the elastic collision cross section 
between electronics and drypton atoms. He finds a collection of data from 
L. S. Frost and A. V. Phelps (1964) PHYREV V00136 1538. To find out 
what further work has been published, assume that new work would cite 
this paper in its bibliography. The computer program is, therefore, 


searched to see if this has occurred. 


First, the reader must log in to the computer in a manner appropriate to 
his connection with the computer. Once logged in, he calls for the entire 
library file with the command TIP ard then asks the computer to search 
for the citation and deliver the information with the command to search 
ALL and find the appropriate citation, by typing as output the author, title, 


and identification. The command actually typed into the console is: 


S ALL from 1964 
F CITE PHYREV V00136 1538 


G PATI mk 
4 } wat 
* Basic Data of Plasma Physics, John Wiley & Sons, Inc., New York, 1959. ’ } 
i ; 
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On the assumption that some new reference was found, the reader then 
may want to browse further in the associated literature. This is most 
efficiently done by having the computer search for all papers which share 
at least one bibliographic reference with the new paper just found. This 
cén, of course, be done against the entire library file, but since the - 
material of 'Basic Data" was compiled using this technique of literature 
search, computer time can be saved by excluding from the search all those 
Papers which already appear in the text. This can be accomplished by a 
subroutine called ''BROWN!'! created in the process of putting this present 
volume together. The command on the console to get a shared bibliography 


based on the Program for this volume would be: 


S ALL BUT NOT BROWN 
F SHARE B PHYREV V00160 1234 
O PATI 


Dr. Myer M. Kessler, Associate Director of Libraries at M.I.T. and 
Director of Project TIP, studied statistically where most of the physics 
literature was published and developed a list of 25 journals which contained 
70 percent of physics papers (Physics Today, 18, 28, 1965). In develop- 
ing the 'BROWN!'! Program for this report, certain journals were further 
excluded as not being appropriate to experimental gaseous electronics and 
plasma physics. The list of journals in the 'BROWN!! program and their 


program code designation is given in the following list: 


Applied Physics Letters APPLET 
Canadian Journal of Physics PHYCAN 
Helvetica Physica Acta “PHYHEL 
Indian Journal of Physics , INDJPH 

Japanese Journal of Applied Physics ' PHAPJA 
Journal of Applied Physics PHYAPP 
Journal of Chemical Physics JCHEPH 
Journal of the Physical Society of Japan PHYSOJ 

Physica HYSICA 

Physics Letters . PHYLET 
The Physics of Fluids 'PHYFLU 
Physical Review A PHYREV 
Physical Review Letters PHYRET 
Proceedings of the Physical Society (London) PHYPRO 
Soviet Physics - JETP SPJETP 

Soviet Physics - Technical Physics SPTPHY 
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In addition the total TIP program will also search: 


Annals of Physics ANNPHY 
Il Nuovo Cimento NUOCIM 
Nuclear Physics NUCPHY 
Physical Review, Series B PHYREB 
Progress of Theoretical Physics (Kyoto) PROPHJ 
Soviet Physics - Solid State SPSOLS 


Another source of data which are not always to be found in the periodical 
literature are papers published in the proceedings of conferences. A 
number of these were of importance to the material of this report and 
these have been especially incorporated into the TIP program. They are 
Proceedings of the Fifth International Conference on Ionization Phenomena 
in Gases held in Munich in 1961 (MUNICH), Comptes Rendus de la vi= 
Conférence Internationale sur les Phénomenes d'Ionisation dans les Gaz 
held in Paris in 1963 (OPARIS), and Proceedings of the Seventh Inter- 
national Conference on Ionization Phenomena in Gases held in Belgrade in 
1965 (BELGRA). 


If the reader does not wish to exclude the already search-shared biblio- 
graphies in the ''BROWN'' program, both lists of journals will be searched 
when ''BUT NOT BROWN'' is omitted from the console command. 


This study is aimed at serving two purposes: one is to bring the ''Basic 
Data'' book up-to-date; the other is to provide a computer method of find- 
ing new data which appear after the report is issued. The second purpose 
is accomplished through the citation-process, and, therefore, will be 
effective only if there is some article to CITE. Thus, if no published work 
has been found which supercedes that published in the original ''Basic Data" 
book, the previous volume's citation will be repeated with its proper 
reference coding, so that subsequent material may be found by the com- 


puter-search method. ; 


This is possible, because, although the TIP program does not contain 


literature published prior to 1959, the bibliographic coding is complete. If 


the reader so desires, he may call for the computer to present him with J wa 
any complete bibliography he desires by calling for ''B" as an output } 
i 


command. 


u 
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In searching the literature, it is sometimes desirable to determine what 
the common bibliographic references are between two articles found by 
the 'SHARE B!! technique. This is provided for in the TIP program by 


asking the computer to print the linkage between any two papers. 


The program has been written so that it is not necessary for a cited 

reference to be for articles actually available in the program by title, 

author, etc., but only that it be-in the references. Thus, the classic 

paper on elastic collision was published by Brode in the Reviews of ‘ 
Modern Physics in 1932, long before the time interval covered by TIP, 

but use of the LINKAGE request will uncover the proper volume and page 

location of Brode's paper if it has been used as a reference in a recent 


article. 
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Baseband Design of a Unified Carrier - Homer C. Peterson 


Work on a design-optimization program was started using FORTRAN and 
batch-processing compiling. The very nature of the equation and the 
technique of handling were not clearly defined and the infinite limits on 
the summation made it hard to determine what intervals to investigate. 
The first runs were so long that they were dumped from the Lincoln 
computer without progressing far enough to even show a trend to the 
results. With usually 20 hours elapsing between submitting a run and 
seeing the results, it became apparent that the problem would either 


have to be over-simplified or dropped. 


At this time CTSS became available touse. The first big step was to 
recode the program in the MAD language. The program was then modi- 
fied to be handled a section at a time; where hand comparisons could be 
made, a rough plot obtained from the teletype, and programming errors 
readily corrected. Within a week the first outputs were. available, while 
the original had dragged out for nearly a year. Considerable rewrite and 
modification was necessary, as it became clear that certain sections of 
the program did not agree either with the model or with the actual world. 
A new Bessel function subroutine was written, to obtain the desired 
accuracy, only to find that it worked about half the time. The CTSS 
loader didn't always properly load a double-precision hardware programm, 
so the program had to be rewritten to be self-shifting if it was started in 


an odd register. 


For nearly a calendar year (but probably less than two man-months of 
work) the program was in constant change with as many as ten recompila- 
tions occurring in one day. During this time, the program was lost only 
three times due to computer failure. The difference between using a CTSS 
system and batch-processing was a solution with about two man-months 


of work as opposed to abandoning the project. 


4 nnaeinermsneivs 


am 


Bag 


a or a 9 ne et ee Sey ne PR = seme Ss ae 


110 LINCOLN LABORATORY 


On-Line Data Storage and Retrieval - John F. Nolan 


During the past year, the on-line data storage and retrieval system (as 
described in last year's progress report) was coded, checked, and 
operated on a variety of sample files. Although limited in data capacity 
by the disk space allocated to the research problem, the experimental 
system demonstrated successful implementation of advanced techniques 
for on-line access and manipulation of formatted data files. The most 
Significant features of the design are: 

1. Storage of coded information with the data files describing 
their OER an and content. 

2, Representation of this descriptive information by the 
identical conventions used to store the data files themselves; thus all 
system commands for searching, ‘cross-associating, and manipulating 
data files can be used to operate upon the file descriptions as well. 

3. The internal use of list-structures to allow unconstrained 
cross-association of filed data (including, as above, the file descriptions) 
according to the user's needs. 

4, Simple one-word commands and dynamic instructions which 
allow immediate use without any indoctrination or programming knowledge. 


(See Nolan, Appendix C) 


Compilation for a Digital Differential Analyzer - Harold K, Knudsen 


Lincoln Laboratory is constructing a Digital Differential Analyzer (DDA) 
as part of a hybrid computer facility. A compiler is being written on the 
MAC system to make it possible for the nonspecialist to generate solutions 
to ordinary differential equations on the DDA, The inputs to the compiler 
will be: g 

1. The diffe#ential equation written in a FORTRAN-like 
language, 

2. The initial conditions of the variables, and 


3. Estimated maximum absolute values of the variables, 


stil... 
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The output of the compiler will be: 
t, An interconnection matrix which specifies the internal 
connections of the integrators and servos inthe DDA, and 


2. The scaled initial conditions of the problem variables. 


The compiler outputs will be loaded directly into the DDA from the MAC 
computer through a channel consisting of a Model 35 teletype and the 


LINC computer. 


On-Line Experimentation - Joseph W. Bex 


During the past year, CTSS was used for the following activities: 

1, Searching for sequences with desirable autocorrelation 
properties, for use in coding patterns for the Arecibo radar antenna in 
Puerto Rico, 

2. Testing the stability of a new circuit design to be used to 
control the carbon arc for Lincoln Laboratory's new environmental 
facility, 

3. Simulating designs for magnetic orientation systems used 


on Lincoln's experimental communications satellite, LES. 


Exploratory on-line computation, done with the responsible scientist or 
engineer present, was performed for: 


1. Designing antenna arrays for use with the LES satellite, 
involving the adjustment and testing of various parameter values, 

2. Determining the stability or instability of the theoretical 
model for a new solid-state device which would have unusual high-gain 


properties, 
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In addition, the Project MAC computer system began to be used by staff 


of the Bell Telephone Laboratories and of the General Electric Company 


participating in the development of the MULTICS system, 


NON-MIT USERS 1S 
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Model Testing - Merrill M. Flood (University of Michigan) 


The purpose of the model testing project is to develop programs and on- 
line techniques adequate for parameter estimation and goodness-of-fit 
testing on several stochastic process models. These models derive from 
the data of various learning and decision experiments with human subjects 
and animals. We have adapted programs used previously on IBM 7094 
computers at the University of Michigan, at the University of California, 
and at the System Development Corporation, for CTSS use. The pro- 
grams are now suitable for one specific stochastic model and have been 


used successfully on a few small test runs. 


One set of FORTRAN Subroutines, called LOOK, are part of the general 
program, STOl, They are an adaptation of a direct search code developed 
at the Westinghouse Electric Company for the iterative solution of maxi- 
mization problems with constraints. The LOOK sub-routines will be use- 
ful as part of other computer codes and include optimization calculations. 
Eventually, we hope to make LOOK and other similar optimization codes 
available to all interested CTSS users. One FAP Subroutine, called 
PRODUC, within STO1 was written to calculate rapidly and accurately a 
lengthy product of probabilities. Some of our current programming is 
devoted to the development of improved codes for this calculation and 
CTSS is useful in testing such alternative codes for speed and accuracy. 
We are using these optimization codes to investigate, with experimental 


data, a few stochastic process models of rat behavior. 
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Discovery and Learnin Techniques and Programmed Instructions - 


Leonard Uhr (University of Michigan) 


We are working ona program that learns how to handle one or two-dimen- 
sional inputs, such as language strings and sensory patterns. With the 
exception of those parts that build up class structures, sections of the 
program that use previously-learned information to generate the appropriate 
response to an input have been largely completed and debugged. The sec- 
tions that learn have been coded, and are now being debugged. We are also 
continuing tests of the earlier program that handles one-dimensional strings, 
and we are looking into methods for increasing its abilities to handle 


. 


different types of problems, 


In addition, we have coded a program that turns written materials, such as 
textbooks and spec’. ly-prepared instructional materials, into rough drafts 
of teaching-machine programs, conveniently edits such programs on-line, 
and presents the finished instructional program to students. In this pro- 
gram, content material is interspersed with short fill-in questions that 
lead to forward or backward branches as a function of the student's success 
or failure in giving a correct or synonymous answer. We can now generate 
teaching-machine programs in three ways: by compiling from materials 
prepared according toa simple standard format; by generating, during - 
interactions with the student, from a standard description of a problem 
domain; and by converting ordinary books into programs. We should now 
like to combine these methods into a single program that would offer all 


these options, 


oo, 
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Generalized Desk Calculator - Anthony G. Oettinger and 
Adrian Ruyle (Harvard University) 


Project TACT (Technological Aids to Creative Thought) has been develop- 
ing the on-line system of Glen Culler, which functions roughly as a gen- 
eralized desk calculator. With this system, not only can arithmetic opera- 


tions be carried out on single numbers, as they would on a typical desk 


calculator, but operations of more generality — such as exponentiation 
and forward differencing — can be carried out on functions considered as 
such. 


Execution is as quick and easy to initiate as using a desk calculator. 
Operations on complex functions are also possible, as are plots of functions 
on a cathode-ray-tube screen, and display of text. There is also a facility 
for on-line storing of programs" of operator sequences, which can »e 


executed en bloc as single calculator operations. 


The calculator is now developed to the stage where it is a useful tool, for 


about a half-dozen users who will be variously: 


1. Preparing demonstrations, suchas prototypes of course 


lectures in mathematics, and giving calculator exhibitions, 


2. Conducting research, such as kidney physiology and 


frog demography, 


3. Exploring the imaginative possibilities of the calculator 


itself. 
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Aspects of Integer Linear Programming - R6mulo H. Gonz4lez-Zubieta 


Under U, 5S. Army Research Office (Durham) Contract No. DA-31-124- 
ARO-D-209, a primal feasible (all-integer) integer linear-programming 
algorithm has been developed and programmed, together with a related 
procedure for obtaining a first feasible solution. Once a feasible solution 
is found, the algorithm maintains feasibility at each stage, in contrast to 
other algorithms that have been programmed and are currently available, 
These other algorithms do not achieve feasibility until the optimal solution 
is reached. The primal feasible algorithm is based ona particular way of 
applying the cutting planes previously developed by R. E, Gomory, and on 


a specific interpretation of their role, 


The finiteness of convergence has been established for two-dimensional 
problems but not for the general case; however, there appears to be at 


least computational convergence in a considerable fraction of the cases, 


In addition, a Generalized Euclidean Algorithm for finding the greatest 
common divisor for more than two numbers has been defined. The solution 
of systems of linear diophantine equations is presented in terms of integer 


linear programming, 
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Introduction 


This interdepartmental laboratory provides facilities for academic research 
in three categories, designated as general physics, plasma dynamics and 


communication sciences. 


Major support for the research is provided by the Joint Services Electronics 
Program of the Army, Navy and Air Force as well as the Atomic Energy 
Commission, the National Science Foundation, the National Institutes of 


Health and the National Aeronautics and Space Administration. 


As indicated by the following section, a substantial number of RLE research 
projects have received invaluable assistance from Project MAC facilities. 
These faculty and student research activities span a wide range of scien- 


tific and engineering subjects. 


Dynamics of Active Plasma Systems - Abraham Bers 


During the past year we have made extensive use of the MAC computer 

and the ESL display system to solve a large number of difficult and impor - 
tant problems in plasma stability. The complexity of these problems is 
such that without direct man-computer interaction many of these problems 
would not have been attacked. The basic problem of plasma stability to 
small-signal perturbations was described in our progress report last year. 
The most significant programs developed during the last year for our prob- 
lems are summarized in three reports that follow (M. A. Lieberman, C. 


Speck and M. Pennell; and J. D. Mills and A. Bers). 


Our project has brought a large number of our students in contact with the 
Project MAC facilities. The thesis research of these students has ranged 
from computer-display system programming to computer stability analyses 
of hot, gaseous plasmas, and plasmas in solids. From our group, a total 
of about a dozen graduate students, three faculty members, and two post- 
doctoral visitors have made use of the Project MAC facilities during this 


past year. J se 
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Dispersion Relation for Hot Plasmas - Michael A. Lieberman 
eS 


We have written a computer program that finds the zeros of a trans- 
cendental dispersion function Diw,k,...) in the complex w plane. Given 
the complex frequency w, the wave number k, and any other parameters 
that a user may desire, he must provide a subroutine within CTSS that 
computes the value of D. This subroutine may be written in MAD, 


FORTRAN, or FAP. 


In operation, the program continuously steps up k by an increment Ak and, 


at each step, finds a zero, ca of the dispersion function D: 
Dw, k+ndk, ...) = 0. 
=n 


The zero, wo is found through the construction of a grid of values in the 


complex w plane around an initial guess aes as follows: 


Gee* provided by the user 
“gl ~ ato 
* en eed * ei 28s sa A eres 
The program evaluates the real function F = | D( (w, kt ndk,...) ks for every 


point on the grid. Whenever a minimum of F is found at wy the gridis 
refined several times until we is given to three significant figures; then, D 
is checked to verify whether its real and imaginary parts change sign in 
the neighborhood of wii and, if they do, the zero is printed and kis stepped 
up. Thus, the zeros of D in the complex» plane are computed as a func- 


tion of the wave number k. 


Without time sharing, a solution of this problem would be prohibitive; with 
it, a user may specify and alter at will the grid size and spacing in the com- 
plex w plane, the wave number k and its increment Ak, and all other param- 
eters. If at any step a zero of Dis not found, the program requests the 
user to change the grid size, spacing, and location in the complex wplane, 
To help the user, the program will print the values of Dand F at the grid 


points if they are requested. 


Ge SE A EEA MELEE IO) <= eal aa See 


o =. 


RESEARCH LABORATORY OF ELECTRONICS 127 


Plasma Dispersion Relations with Infinite Roots - Carlton E. Speck and 
Martha M. Pennell 


In the study of a particular class of plasma waves, the following equation 


must be solved for its roots, VC = RVC + jCVC: 


oe) oo 


el a? (p) q° (p) Z n°32(p) 3 
1 - (VPC) >) n=1'P! = “n¢ IP! pe, tee) 2% 
Lee (VC - n) P wares (VGe-n) 


We had to determine those values of B42 and VPC which had ranges of p 
that resulted in roots with a negative imaginary part. Since the equation 
has an infinite number of roots and is not expressible as a polynomial, 

we used Newton's method to locate the roots. However, Newton's method 
requires an accurate guess of the root location. To obtain these guesses, 
we wrote subsidiary programs for the approximate portion of the root. 
These approximations derive from an approximate equation in only those 
terms of the infinite sum which are large when lvc} ~n. With such 
''guessing'' programs an operator can find the required approximation for 
the more exact Newton's method and, by varying VPC and B+, can rapidly 
search the VC - p plane until he finds the approximate regions of the nega- i 


tive imaginary part of VC. 
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Stability Analysis of Dispersion Relations - James D. Mills and 


Abraham Bers 


A program has been waateene checked out, and is now being used by mem- 
bers of the Beam-Plasma Group of the Research Laboratory of Electronics 
in the stability analysis of plasma dispersion relations. The program 
makes use of the ESL Display Console for cape outputs. The time- 
sharing system facilitates man-machine communication for changing pa- 


rameters and equations in the analysis. 


The motivation for the prog’am came through the formulation of the Bers- 
Briggs stability criteria. These criteria enable one to analyze the disper- 
sion relation of a plasma (or other) system, to determine and classify the 
kinds of stable and unstable wave responses which are inherent to the sys- 
tem. The technique involves the general problem of mapping a function 

of a complex variable. The dispersion equation of the physical system to 
be analyzed, D(w, k) = 0, relates the two complex variables which are 
ww t jo; and k = ke + jk; One then maps lines of constant wore, into 


the complex k-plane or lines of constant kD or ky into the complex w-plane. 


The program restricts the user to polynomial-type equations of no more 
than degree twenty. The coefficients of the terms of the equation are speci- 
fied simply in a subroutine provided by the user. He then translates this 
subroutine using the MAD command. It is then loaded together with the 
display program which queries the user concerning various parameters. 
When these are specified, the program makes all of the calculations neces- 
sary to the mapping and displays the results on the Display Console. In 
Figure 9 for example, the user specifies the lines indicated on the complex 


w-plane (a) and he sees the complex k-plane (b) plotted on the display. 


The program has been used in various stages of development since April, 
1964. A joint RLE - Project MAC technical report is being prepared for 
publication. (See J. D. Mills thesis entry in Appendix B.) 
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Use of CTSS in a Plasma Physics Experiment - David T. Llewellyn-Jones 


This work combines automatic data acquisition and processing with a labor- 
atory experiment in plasma physics which requires rapid feedback between 
experimenter and his apparatus. The last part of this work will establish 
an interface between a digital voltmeter, connected to the laboratory ex- 
periment, anda teletype, connected to the MAC computer. Since a tele- 
type has been installed beside the experiment, the experimenter has been 
able to simulate actual runs and at the same time, to act and manually 
type the data into the teletype as it is displayed on the voltmeter. In this 
way, speeds of approximately one, three-digit data-point in every 2.5 
seconds are possible. Since the experiment can run at speeds up to about 
one data-point-per-second, realistic simulation is only possible for slow 


runs. 


In simulation runs, the . TAPE. feature of CTSS is very useful. This fea- 
ture immediately creates a permanent data record before the data reaches 
the program, and it has all the useful correction and substitution facilities 
afforded by the INPUT mode. When the data was read directly into the 
program and the necessary correction and substitution facilities had to be 
written into the program itself, the program became cumbersome and slow. 
After installation of a reliable automatic recording system, the data will 

be read directly into the program and the program itself will automatically 


create a.TAPE, of the data as a record. 
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Analysis of Speech - John M. Heinz and Eleanor C. River by 


As part of a study on the constraints which the vocal organs impose ona 
speech signal as a result of their physical configuration, we are trying 
to model the process of speech production from its description by lateral 
cineradiographs up to the resulting acoustic signal. Four stages exist in 
analyzing these lateral cineradiographs and relating them to the corres- 
ponding quasi-static acoustic spectra: 1) obtaining from the film a frame- 
by frame tracing which includes all of the relevant detail, 2) devising a 
reference coordinate system for quantitative measurements which can be 
precisely defined in terms of stable landmarks readily observable on the 
film, 3) devising, from a knowledge of the structural anatomy, a set of 
transformations by which the acoustically-relevant cross-sectional area 
may be inferred from the lateral measurements, and 4) calculating an 
acoustic spectrum from a specification of the cross-sectional area as a 
function of distance along the vocal tract midline. The programs that have 
been written to aid in the construction and testing of the models of speech 


production correspond to parts 3) and 4) above. 


A preliminary version of the midsagittal-to-area function transformation 
program is now in operation. With this program, we can preset the cross- 
dimension-to-area function transformation for a particular talker, either 
by reading in data tables prepared from anatomical measurements of him 
or by specifying certain parameters which determine an analytic descrip- 
tion of the transformation. One of the objectives of the present study is 

the investigation of various possible types of analytic descriptions. Once 
the program has been preset for the talker, cross-dimension data from 

the films are read in. The program then calculates both the cross- 
sectional area as a function of distance along the vocal tract midline and 


the effective vocal-tract length. 


We have written a second program which accepts cross-sectional area func- 
tions derived by the first program as input and, with a Webster's equation 
approximation to an acoustic description of the vocal tract, determines the 
corresponding acoustic output spectra. We have then compared a computed } % 
spectrum with the spectrum of the actual sound produced at the time of the 


x-ray filming. Such a comparison tests the validity of an investigated model. 
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Simulation of the Human Larynx - Thomas H. Crystal 


During the past six months, we have developed a set of programs for sim- 
ulation of the behavior of a human iarynx during the production of speech 
sounds and for control and display of the results of the simulation on a Te- 
mote typewriter console. The programs still need some additional modi- 
fications, but have produced some interesting eeente: 

The larynx itself is represented by two systems. A mechanical system, 
which models the mass and tension of the vocal cords, is used to deter- 
mine, as a function of time, the opening between the vocal cords and air 
flow through the vocal cords. The parameters of this system depend on 
the vocal cord opening as determined by the mechanical system. Differ- 
ence equations that describe these systems are iterated successively to 
produce the simulation. In addition, the program contains provisions for 
the connection of any of a number of vocal tract and sub-glottal models to 
the larynx model. These systems should together provide a complete pic- 


ture of larynx activity. 


A feature of the simulation programming is the ease with which the prog- 
ress of simulation may be controlled and its parameters varied. This 
progress, together with various plotting, printing, and frequency analysis 
routines, is controlled through simple, one-word commands. Thus, a user 
may perform a simulation, examine the various waveforms generated, per- 
form Fourier analyses on any particularly interesting waveforms, change 
the parameters controlling the simulation dependent on the results of the 
previous simulation, and simulate again. Because of extensive control 
facilities programmed into the model, full advantage of time-shared com- 


putation can be obtained in this simulation problem. 
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Articulatory Events of Speech Generation - William L. Henke 
+ pe ecn Veneration 


This project applies computer information-processing techniques to the 
study of the articulatory events of human speech generation. Our data, 
which consist of cineradiograph films of human speakers, together with 
their corresponding phonemic input and acoustic output, are studied and 
modeled mechanically for the purpose of understanding the human speech 
process, in particular the transformation between a discrete phonemic in- 
put and a continuous acoustic Output. As we are dealing with Spacially- 
distributed structures, we require the input/output facilities afforded by 
the ESL Display Console. Since the information-retrieval and modeling 

of these data depend upon techniques of on-line interaction of a user with 

a computer and non-textual types of information, part of our work is the 
development of such techniques. Three classes of programs have been 
developed to the present: 1) cineradiographic graphical data input, dis- 
play, and analysis; 2) implementation and display of a model of tongue 
action; and 3) a very general solution for the acoustic transfer function of 
the vocal tract, with associated displays and inputs. They are currently 
disjointed, but further work will increase the communication between them 
so that they may be used Singly or together, depending upon the needs of the 


aspect of speech being studied. 


One part of our work involves the input, storage, and display of a graphical 
representation of the mid-sagittal plane of a human vocal tract in action. 
The MIT Speech Group has recently made frame-by-frame tracings of pro- 
jections of cineradiographic films of speakers, These'tracings are fortu- 
nately about the size of the ESL Display Console scope face, so that with 
tracing paper affixed to the face of the scope, the light pen can see the 
tracking cross through tracing paper. We enter points and connected line 
segments for various anatomical structures and landmarks, together with 
identifying information, into a dynamic data structure. At the end of the in- 
put for one frame, the data are edited and written out on the disk. The first 
generation display program searches the file for desired frames and si- 
multaneously displays up to three at different intensity levels. The organ- 
ization of this file and the specification and recording of complementary 
acoustic and phonemic data, which wil) be suited to analysis by visual com- 


parison, has prompted our studies on graphical input. 
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We are developing a model of articulatory action which will distinguish be- 
tween operators that are associated with inputs representing linguistic 
categories, and states that are associated with the responding articulatory 
structures. Considering the behavior of the tongue as a responding struc- 
ture, our work develops general subprograms for the display of graphical 
and textual data and constructs a time-domain simulation model which 
stops at any desired time and displays the present state. With this model, 
a user can then compare the state with x-ray tracings and modify some of 


the parameters of the model. 


For acoustic comparison, an algorithm which calculates the frequency - 
domain response of the vocal tract is desirable. We have written one for 
general situations which calculates both the complex spectrum from input 
data of the area function and the loss as a function of frequency, position, 
and area. From the complex spectrum, it further calculates poles of the 
transfer function. Since all this information is displayed ona CRT, the 
volume velocity distribution at any frequency is available. Programs 
developed by other members of the Speech Group to obtain area functions 
from mid-sagittal plane representations of the vocal tract will combine 
with the programs discussed above. 


o 


Grapheme-to- Phoneme Translation of English - Francis F, Lee, 
’ 


It appears that when an English-speaking person ' 


‘sounds out'' a printed 
English word or an English-like non-word, he makes use of several often 
mutually-supplementary processes, namely, recognition of the pattern as 
a whole, identification of certain affixes, in particular the suffixes, iden- 
tification of compound words (e.g. Whitestone as white-stone, not whitest- 
one) and a mapping of graphemes to phonemes employing the graphemical 


context, with the last process occurring at the latest or innermost level. 


The idea of a mechanical procedure, translating a printed English word 
into its phonemic representation using only a single and small exception 
list, runs into immediate trouble. By using multiple exception lists, con- 
sulted when necessary, the total number of exceptions and the list of con- 
text sensitive rules can be substantially reduced. The extent of the neces- 


sity of affix identification is being studied. 
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It appears that the presence of mute e in medial position due to word com- 
pounding, such as spacecraft can only be detected by some sort of table- 


lookup procedure. 


The 18,000 more often used English words (from Thorndike and Lorge) have 
been processed to show the grapheme-to-phoneme cluster relationship in 
graphemic context. This body of data has been used for the creation of 
exception lists. It will be used for the decomposition of complex clusters 


and the formation of rules for the previously-mentioned innermost level. 
Programming Support and Development - Martha M. Pennell 


The Computation Group has provided both programming support and devel- 
opment for users of the time-sharing system. Our program QUIXOT, 
written to alleviate the programming load of the Computation Group, has 
been operational for almost an academic year. A user need only type an 
algebraic expression for the coefficients, and QUIXOT writes, compiles, 
loads, and runs a program to find roots of a polynomial of degree twenty 


or less. 


Encouraged by the success of this program, the Computation Group modi- 
fied QUIXOT so that the user can type in his coefficients as one expression, 
rather than as real and imaginary parts. Because MAD, the source lan- 
guage for this program, does not provide for complex arithmetic, a sub- 
routine call must perform each complex arithmetic operation. Also, we 
incorporated a compiler-writing program into the QUIXOT program, which 
saves the user some initial algebra. An outgrowth of this program is 
COMCAL, a routine to evaluate an expression in complex arithmetic. 
Similar in operation to QUIXOT, it first queries the user for his symbols, 
parameter names, and algebraic expression; then, with the command- 
programming feature of time-sharing, it writes, compiles and executes 


a program to evaluate the expression. 
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RLE Statistics on CTSS - Martha M. Pennell 


I cite examples which might give some idea of the time difference involved 
in a programmer's speed with CTSS versus batch processing. They are 


taken from the work of RLE Computing Group. 


A mathematics major with one summer's FORTRAN experience was hired 

as a scientific programmer. She was given a short demonstration of CTSS 

— compute A from the formula A = B+ 1 where B is input from the tele- 

type — and told to duplicate it where B is an array of N numbers to be 

given at execution time. She did so in about four hours. The next day 

this person was given an RLE problem, and the calling sequences for the 

two necessary subroutines (a Bessel function and roots of a polynomial), 

and she was told to code it in MADTRN. Within three days her program pro- 
duced results. At first it contained errors, but she had set up the subroutine 
linkages correctly and had made no errors inits logic. Her errors lay in 
the algebra for computing the sum. This initial coding was done witha 
minimum of supervision. When I later asked how long she thought it would 
have taken her on batch processing, she said a month, because her errors 
were due to her unfamiliarity with FORTRAN and a series of programming 
blunders, such as keypunching errors and omission of control cards or 


subroutine decks. ‘ 


Another user has five years of experience in FORTRAN and one year with 
MADTRN. The majority of programs she has had to write on CTSS can be 
coded and debugged, on the average, with a two-hour session at the console, 
These same programs, if run and checked out under batch processing here 
at M. I. T., would take at least four to six days. These programs involve a 
main program and possibly two or three subroutines of about 40-50 state- 
ments each. These programs have found roots to polynomials, evaluated 
Bessel functions, performed numerical quadrature, and so on. She 

finds that with programs of this length, three or four tries usually debugs 
them. If you take into account machine errors, bad tapes, lost output, this 


usually means about a week in a batch-processing scheme. 


This person also decided to learn MAD. The problem she chose to code was 


atree sort. Giventhe people inher office, their sex, their marital status, 
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their laboratory affiliation, and their telephone number, the problem was to 
to write a program to sort them, grouped alphabetically, according to one 
or more of their characteristics. The MAD program was about 40-50 
statements long and took six tries before it finally compiled, because the 
programmer was not completely familiar with the punctuation and rules of 
Syntax. Two more tries debugged the program logic. This work was done 


in about two hours of console use, 


Sorting of Personnel Records - Robert L, Rappaport 


The RLE personnel roster includes approximately 700 people, who are 
classified by several categories: position, research group, rooms occu- 
pied, telephone extension, and soon. For example, an associate pro- 
fessor of the Electrical Engineering Department who works in the Plasma 
Physics Group must be listed under each of these files. At present, these 


files are brought up to date and sorted manually. 


We have written a computer program to automate this Organization. The 
program Operates on a masterfile which contains all the essential informa- 
tion for every member of RLE, so that we need make correction only once. 
Sorting can be done on any of the items entered. A second program prints 
an alphabetical directory of all RLE personnel. We have checked out both 
of these programs for a limited number of names and will shortly extend 


them to the entire RLE roster. 
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COMIT - Victor H. Yngve 


During the year we maintained the COMIT system as a standard program- 
ming language available to users of MAC, and we further developed the 
language so that a new system, COMIT II, will soon be released for use 

in MAC. COMIT II includes many new features and conveniences for the 
programmer, while any program that works in COMIT will also work in 
COMIT II. We also developed a program VEDIT, which is written in COMIT 
and is especially useful in the writing and editing of COMIT programs. It 
introduced a number of new features that have found their way into recent 


edit programs. 


Jared Darlington further developed his set of COMIT programs that check 
the validity of arguments expressed in English. These programs are 
unique in that they incorporate several levels of analysis of the input argu- 
ment to determine the lowest level of logical analysis sufficient to prove 
the argument valid. In addition, a program in COMIT which partially 
translates programs from FAP to GE635 GEM has been developed by two 


undergraduate students. 
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Computer-Aided Design - Charles Garman, Richard Bayles, and 


Harold Levin 


During the past year, Mr. Bayles and Mr. Levin were programming the 
section of the CTSS supervisor responsible for maintaining the display 
on the ESL Display Console. Since the design of the software for the 
display system requires a fairly rigid and highly artificial structuring 

of display data, operation of the first version of the supervisory program 


svon showed defects in this design. 


We began in August, 1964, to redesign the interface between the user 

and the display system and succeeded in developing a new supervisory 
program. Among its features were: the recognition of the slave display, 
which had been added to the display console, so that two users can main- 
tain different displays simultaneously; the addition of a system which 
assigns real-time tasks to various analog and digital inputs of the display 
console and to different sections of the image; and the addition of sub- 
routine pictures, which facilitate display of similar graphical objects. 
This project was completed in December; and, with the help of Michael 
Brescia of the ESL Display Group, we have continued to modify and 


improve the system, 


At the end of December Mr. Bayles began programming for the Disk- 
Drum Strategy module that was to be incorporated into the revised File 
Control System of the CTSS supervisor. This work was continued and 
has since been completed by Mr. Garman. He was also involved during 
the previous summer in programming a system that generated and dis- 
played on the ESL console that were surfaces derived from an algorithm 
of Professor S. A. Coons. These surfaces are completely general in 
nature and may model, with remarkable precision, objects as diverse as 
boat hulls, an octant of a sphere, adjacent planes of a regular tetrahedron, 
or a section of an airplane fuselage. In December, he began work on 
modifications and extensions to the CTSS 'ED' command, which is a con- 
text editor for BCD card-image files and he is reprogramming the input- 


output section of 'ED!' to conform to the new file-control system. 
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On-Line Mathematical Anal sis - Roy Kaplow, Stephen L. Strong, and ‘\ 
; John W. Brackett 


With our MAP system a user with little knowledge of computers or pro- 
gramming can solve a wide variety of mathematical problems expressed 
in normal mathematical terminology, MAP contains many interrelated 
programs, only a small fraction of which may be executed in response to 
any one request. The language consists of functional equations and an 
expandable group of commands. Al] communications can be transmitted 


through a standard Project MAC console, 


The system evaluates one-dimensional functional equations that involve 
constants, vectors, library functions and ordinary arithmetic Operators. 
The form used is identical to normal notation, except that multiplications 
only implied in the first form must be specified in the latter and 

functional operations, such as sine and cosine, must be designated by a 
terminating f. MAP automatically generates the q(x) vector at equal 
intervals over any range in x specified by the user. If the range and 
interval in x, the values of the constaut a, or the values of the vectors y (x) 
are not already specified, they are requested by MAP and may then be 


typed in under the control of context editor, 


Commands in MAP correspond to requests for Operations upon functions, 
such as q(x), or upon groups of functions, which cannot Proceed in the 
point-by-point fashion of equations, With such commands at his disposal, 


the user need not remember the parameters required for the Operation, 


As an example, the typing of the single word "integrate will cause MAP 
to: ask the user if he desires to integrate between two constants, the full 
range of an independent variable, cr some combination of constant and 
variable limits; ask the name of the function to be integrated; the name to 
be given to the constant or function answer; and, if necessary, the values 
of the limits. As with equations, if the function to be integrated is not 
already defined, its values will be requested and may be typed in with a 
context editor. The sophisticated user may prescribe all of the necessary 
information, regarding options within the Operation and parameter names 


or values, and avoid unnecessary communication. The following operations 


nan eee a Ss SO rate ie A er ee Coe scien te Bates et ni SR Ada 5 Mm | 


SCHOOL OF ENGINEERING 143 


are now available: integrate, differentiate, convolute, Fourier trans- 


form, change of basis, and least-square fitting. 


In addition to the equations and simple commands discussed above, other 
commands are available which facilitate problem solving and which extend 
the range of MAP. These are briefly described below under five separate 


classifications. 


A. DATA HANDLING 

These commands provide for off-line input of blocks of data, for the 
deletion of data from mass memory, and for the saving of data and pro- 
grams. All information is handled by its symbolic name, and the user 


need not be concerned with computer techniques. 


B. FUNCTION HANDLING 
The ''minmax'' command changes the range and interval for which an 
independent variable or function is defined. The ''select'’ command 


manipulates a portion of an individual function. 


C. OUTPUT 

The "print'' command is used for all output other than system-generated 
messages and questions. The format of printing may be specified by the 
request of auser. When he so desires, a user may request that only 
part of the tabulated range of a function be printed, together with an 


interval in the independent variable. 


D, COMMAND SEQUENCE CREATION AND EXECUTION 

The ''run", ''create'', and ''edit'' requests define a sequence of MAP state- 
ments and store them in mass memory for subsequent use, ''Create" 
allows typing in a sequence of commands under control of the context 
editor. The sequence is assigned a name and executed by "run'! and the 
assigned name, ''Edit'' edits a sequence of commands to correct errors 
or to substitute variable names when the sequence is executed by a "run!! 
command. A user who has not learned a complex language, such as 
FORTRAN or MAD, can thus create a programming structure for his 


specific application and call it by its symbolic name. 
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E, EXECUTION OF MAD PROGRAMS 

With "Execute," a user can write his own command in the MAD language 
and use it within the MAP system, The existence of this command makes 
the system easily expandable and should encourage users to add to the 
system. After the user has written a MAD Program, such as "root," he 
need only type "execute root! to run his Program. As an aid in the 
creation of his own commands, a large number of systeni subroutines 
are available to the facilities of the MAP system. Other subroutines 
accomplish such tasks as: loading any function, previously defined 
within the MAP system, from mass memory into any array specified in 
the user's program; obtaining one or more values of a function from mass 
memory, and storing them as the values of any variable in the MAD pro- 
gram; creating a MAP function from an array in the user's program: and 
executing any CTSS command. All of the necessary subroutines are 
automatically loaded along with the user's Program when the "execute" 


command is given. 


During the past year the MAP system was made available to a number of 
users of the MIT Computation Center. Many of the alterations and 
additions provided during this year are the result of our experience with 
those users. Plans for the forthcoming year include extension of the 
system to functions of more than one independent variable, matrix opera- 


tions, differential equations, and addition of graphical input and output. 
Derivation of Preliminary Ship Lines - Theodore M. Pitidis-Poutous 
SE minary ohip Lines 


The object of this investigation, supported by Bureau of Ships Contract 

No. NObs 90100, was to develop a suitable procedure for the mathematical 
derivation of preliminary ship lines, for a new design, by reference to 

the lines of a parent ship. An algorithm has been formulated to determine 
a new curve or surface, in the form of a polynomial, resembling as 
closely as possible another graphically-specified curve or surface, and 
Satisfying a prescribed set of constraints. The computat‘ons involved are 


facilitated by the use of Special polynomials. 


Various ways of applying this algorithm to the delineation of preliminary 


ship lines have been investigated. In particular, a procedure has been 
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formulated whereby the preliminary lines for a new design can be 
Obtained in two steps. First, the parent form is fitted with a polynomial 


expression, Then the surface equation for the new design is obtained by 


scribed hull form parameters. As an example, the hull form of a trawler 
was first fitted and then modified to conform to new values of prismatic 


coefficient, maximum area coefficient and waterline coefficient, 


on the assumption that the surface could be considered as smooth. It is 
possible to express the hull form in terms of Legendre, Chebyshev, or 
ordinary polynomials, whichever is advantageous in the particular case, 
while the fitting technique allows the number of terms of the expression 


to be chosen according to the required accuracy, 


the preparation of input data was thought to be justifiable, since better 
fitting could be obtained with a smaller number of terms. The approxi- 
mations obtained were deemed acceptable for preliminary and feasibility 
Studies. Undesired Oscillations, which are unavoidable with high-degree 
poiynomials were minimized with the fitting technique and appear to be of 
tolerable magnitude, Polynomial approximations to ship hulls could 
perhaps be improved further only by transforming the hull su rface toa 


more ''respectable"' one before the fitting is undertaken. 


Once a good approximation to the parent hull form was obtained in terms 
of polynomials, the second step involved deformation or modification of 


the parent form according to new form Parameters, 


The use of additional constraints were useful in making radical changes 

to the parent form, It is possible, at least in Principle, to incorporate a 

bulk by applying suitable constraints to the forward sections, and by suit- } 
ably modifying the sectional area and designer's waterline curves, In 


the trawler example, nine such constraints are necessary for each section 
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in the region of the bow, so that the corresponding parent section is 
“Y ignored. Further investigation concerning the use of the constraints 
already employed, or possible additional constraints, will be necessary 


to critically evaluate this procedure. 


These remarks apply to the trawler case examined so far, as full ship 
forms have not been investigated. Preliminary experience with very 

full sections indicate that a few more terms would be required for approxi- 
mations comparable to those obtained for the trawler. Further investi- 
gation is necessary, however, before definite conclusions can be drawn 


with respect to full forms, (See Pitidis-Poutous, Appendix C.) 


Time-Sharing Reactor Code System - Kent F, Hansen, Ian C, Pyle, 


Billy V. Koen and Q. L. Miller 


The purpose of this work is to study the use of time-sharing in Nuclear 
Engineering education and research, Ina question-and-answer session 
at the console, a user derives input data for any nuclear code, in such a 
way that he need have no knowledge of computer programming at all, In 
the past year, we have completed a general-purpose input/edit program 
that obtains input data from the console and numerous nuclear codes. 
The entire collection of Programs is called the TREC system (Time- 


Sharing REactor Codes), 


The TREC system overates as follows: when the user requests a particu- 
lar code, TREC obtains the appropriate files and certain satellite routines 
that read in the input data. The input/edit program asks the user for each 
item of input in a language natural to the problem. The user may review, 
change, or correct his data. After all data is properly given, he may 
request the code to start; TREC then transfers data and control to the 
nuclear code, At the conclusion of the calculation or after an interrupt, 
control is returned to the input/edit program for possible changes in data. 
At the conclusion of the run, the user may select another code, 
+ 

{ The function of the routine TREC is to ask the user what code he wants, 

determine if the code is available or not, and, if it is available, locate 


the symbol table and test routine for the code. The symbol table itself 


' 
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is a list of encoded variable names and desc riptive phrases that define 
each variable. The data is stored in a linear array for delivery to the 
nuclear code itself. The test routines are used to check all input data. 
Occasionally a block of data may be standard information, such as neutron 
delayed -yield fractions for fissionable species, in which case the test 
routine also provides the standard data. Upon completion of the input / 
edit phase, the input data array is written out ina pseudotape for delivery 
to the nuclear code. Within the symbol table, the name of the file con- 


taining the nuclear code is stored for transfer of control. 


Finally, the routine QUIN is loaded. QUIN scans the symbol table and 
prints comments and questions on the console. Thus, the symbol table 
contains the names of all input data along with additional details. The 
test routine is used to test the validity of each piece of input data. If an 
impossible or inappropriate value is given, the test routine prints a 
comment to that effect. Such details as decimal-to-binary conversion 


are also contained in QUIN. 


Several nuclear codes have been prepared and are now routinely available. 
A brief description of each is given below. 

1. KINET solves the infinite medium-kinetics equations with a 
variable number of delayed neutron groups. 

2. DIFFUS is a one-dimensional, multigroup-diffusion code. 

Its program handles up to 16 groups, 5 regions and 5 isotopes per region, 
and also a maximum of 100 space points. 

3. CROSEC is a code that computes group-average, nuclear 
cross-sections. 

4. The Cross Section Library is not a code but rather a large 
block of nuclear cross sections of many isotopes. The data is encoded in 
a form used by CROSEC. 

5. DNPORT is a double Py transport code. This code analytically 
solves the transport equation for one-group slabs. Although it is written 


and debugged, the symbol table is not yet completed. 
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A Stress-Analysis Program - Richard P. Parmelee 


During the past year, I finished and tested a two-dimensional stress- 
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analysis algorithm for thin domes or shells, and then extended the 
algorithm to handle three-dimensional stress analysis of sculptured bodies, 
such as crankshafts. [| completed this three-dimensional program in 


March of 1965 and since then have solved test problems. 


The object of my research to find a method for stress analysis of 
virtually any geometric region, and then implement this method with a 
computer program that will permit a user who knows nothing about com- 
puters to define his problem and obtain accurate answers, Although I 
found an algorithm suitable for analyzing a large class of problems, I 
had problems with the implementation of this algorithm. The availability 
of a time-shared computer made possible the coding and debugging of this 
program in a reasonable time, as well as the guiding and error checking 


of the user input. 


Although I did some graphical display work with the ESL Display Console, 
the input, output, and program control was, for the most part, restricted 


to a teletype. 


Stress-Analysis Coniormal Mapping - Paul A. Wieselmann 


> 


Fundamental to the solution of plane-stress -analysis problems by complex 
variable methods is the generation of a function to conformally map the 
unit circle, or upper infinite half-plane, onto a region described by the 
physical domain of 'the problem. In general, only approximations to 
physical-domain geometry can be obtained, such as an approximation by 


a polygon. 


My work in the past year has attempted to implement analytic methods to 

obtain the mapping function, in particular the Schwarz-Christoffel integral 

My, that maps the unit circle onto the interior of an n-sided polygonal domain, 
( which is restricted to be convex. Riemann's mapping theorem applied to 


the Schwarz-Christoffel integral allows three of the n+z coefficients to be 
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fixed arbitrarily. The program first estimates values for the remaining 
n-1 coefficients, determines an improvement, and continues until it attains 
a predetermined accuracy. At present, this interactive scheme is com- 
pletely free of manual intervention, except in a few special cases which 


require much careful and tedious intervention. 


I shall free the iteration scheme of manual intervention and make complex- 
plane numerical-integration techniques, especially in the vicinity of a 

pole, more accurate in the evaluation of the coefficients. A result of 

this increased accuracy will be an improvement in the iteration scheme, 
since one of its present failures is due in part to an accumulation of error 


in the determination of the length of the boundary. 


Computer-Aided Teaching of Dynamic Systems Behavior - Ronald C., 


Rosenberg 


The three main aspects of our work in the past year were: the continued 
development of the ENPORT system for use as a simulated dynamic - 
system laboratory, the extension of part of its programming algorithm to 
a hand procedure for student use in analysis, and the introduction of a 


simulated blank box as a significant teaching tool. 


The ENPORT system can simulate the behavior of linear, lumped-parameter 
dynamic systems up to fifteenth order. Its input/output request forms have 
been simplified, but have been made more flexible by improved program- 
ming techniques. In addition, a local digital memory to the teletype- 
station-associated analog computer, which are automatically setup by the 
ENPORT program, ‘provide a signal-versus-time oscilloscope display that 


represents the dynamic system behavior. 


The analysis procedure starts with a bond graph model, assig:.s computing 
causality and a sign convention, and selects state variables. ‘fhe system 
equations may then be written from the graph in an organized form and 
reduced to a standard form algorithmically. By means of a-black box, a 
user can specify the input and observe the output in order to discover the 
unknown system,-. The system has been successfully operated as an on-line 


aid to classroom teaching in a class given by Prof. H. M. Paynter. 
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Automatic Network Synthesis - Marc R. Weinberger 


In the past year I attempted to find good search strategies in a parameter 
space. A strategy achieves optimum performance through a choice of 
several possible network structures followed by a search for a best 
parameter combination. It is called optimum if it requires the fewest 
trials or the least time on a computer. I found a strategy which is 
optimum under certain conditions and often substandard in other cases 
where the required conditions are not met. Instead of conducting a single 
series of trials in the entire search space, this procedure conducts a 
multi-stage search in ever-narrowing subspaces. Simple rules indicate 
the allocation of the entire search to successive stages. This procedure 
leads to a large decrease in the number of trials and is more efficient 
than gradient-type methods or other deterministic strategies in several 


nontrivial cases. 


I also studied the effect of noise in the measurements or calculations. 
Since the random search appears to be less affected by these errors than 
gradient-like methods in several instances, it has a broad field of applica- 
bility. The designer must be able to override the automatic search pro- 
gram, whenever he wants to try out his own ideas. Indeed, he learns 
more about his networks during the early stages of design. Since the 
randomness of the chosen trial points often yields good results from the 
very beginning, the designer may want to cut off the automatic routine in 
order to investigate more closely the excellent early returns. Thus, an 
interaction between man and machine, such as the one afforded by time- 


sharing, is a very important part of this method. 
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A Computer Model of Kinesthesis - George Piotrowski 
The purpose of this work is to experiment with the use of the kinesthetic 

in communicating between man and computer. The mode! of kinesthesis 

will be a radially-extensible pivoting arm which has digital shaft encoders 

for pcsition measurement and computer-controlled torque motors to pro- 

vide a force output. The arm will cover a field of about fifteen inches. 

Square and be ab le to exert about five pounds of force. The first version 

of the arm will be as simple as possible and have most of the computa- 

tional functions carried out by software instead of hardware. Thus far, 

I have written a program to simulate the device and also a control program 

for the device; and I am using the two to ascertain the effect of several 

design parameters, notably the sampling rate and the accuracy of measure- 

ment, on overall performance. With the control program, a user can 

select various mechanical systems for simulation. This program has 
been written in a general fashion, so that other algorithms can be added en 


very easily. we 


Computer Solutions for Boundary-Value Problems - Coyt C. Tillman 


I am developing a general, easy-to-use system for the solution of 
dimensional boundary-value problems. The programs written so far area 
command processor, together with associated I/O and free-storage 
routines, boundary shape description routines, and a set of mathematical 
routines for the interpretation of arithmetic argument strings. The 
command processor, though written with boundary-value problems in 
mind, will be more generally usable. It admits to completely format- 
free input and, at the same time, avoids the imposition of the stringent 


punctuation conventions of other format-independent input interpreters. 


A macro facility provides for the definition of new, composite commands, 
the modification of the system's vocabulary, and so on. I am now attend- 
ing to the requirements of the data structure needed to represent the 


difference equations by which boundary value problems are being modelled. 
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Algebraic Expression Compiler - Grover C. Gregory 


I am developing a compiler for algebraic expressions. It is primarily 
intended for use under the supervisory system that Coyt Tillman has des- 


cribed in this section. 


The most recent work has been on the logic of the compiler's second pass. 
This pass is responsible for producing the appropriate machine code, 
using as its input the tables generated during the first pass. The logic 

is similar to that found in the MAD compiler which interprets lists to 
decide how to arrange the machine code. The aim has been to retain the 
flexibility of using such limits while reducing the amount of core they 
require. To that end, a different logic is used to handle the triple inter- 
actions. In addition, the compiler is constructed so that it sees variables 
and unsubscripted functions as identical, a requirement of the system 
under which it will operate. Testing of the whole compiler, as well as 
modifications to allow handling of subscripted variables and Boolean con- 


ditionals, will begin soon. 


Algorithm to Aid Logic Design - Gerald A. Maley 


My main interest at the beginning of this project lay in the development of 
a computer algorithm for the design of minimum Nor-Nand circuits from 
Boolean expressions. Algorithms do exist in this field for the limited 
case where complements of the input variables are available; but, in the 
design of an actual computer, these complemented variables are seldom 
available. When complements are available, only two levels of logic are 
required; but, in the general case, that of no complements, three levels 
are required. With the available computer time, I was able to develop 


the required algorithm for this type of implementation. 


The proposing of an algorithm for a problem such as this is relatively 
simple, but the proving or disproving of such an algorithm can take months. 
On a remote console ] was able to propose an algorithm and run several 
hundred cases in an hour or so. With response times of this order, 1 was 
able to find delinquencies and modify the algorithm in a matter of days, 
rather than months. Iterations of this kind finally produced a simple 


algorithm that I would not have uncovered with normal machine methods. 
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Social Systems Analysis 
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Social Systems Analysis - Ithiel D. Pool 


We are currently working on a generalized editing system for social 
science data. Over past decades the handling of accumulated large files 
of survey records and census records have presented two major problems: 
the adaptation of older files to a computer, and the presence of errors in 
the data that will mislead a computer program. We are developing editing 
routines for on-line revision and correction of these older files. The 
editing system needs to compare data records with the protocol of 
allowable forms and provide feedback with respect to the editing decisions. 
The first version of the editing system is now being used, and, the ex- 
perience of its users is very helpful in the current design phase of the 


next version. 


The social science data surveys and census records on which we focus con- 
sists of both a codebook and records that describe individual people in ways 
outlined in the codebook. The editing system reorganizes the record, and 
a set of programs handles the codebook. The codebook requires that: 1) 
the data be put into computer readable form; 2) parts of the data be tagged 
by people to indicate characteristics of its form and content; and 3) the 
data be checked for internal discrepancies and the user informed of these 


on-line so that errors may be corrected. 


In the CONCOM Project (ARPA Contract No. 920F-9717) we try to simulate 
the flow of messages in a nation. We first represent the audience by creat- 
ing a population of hypothetical individuals that would constitute a fair sample 
of a real population. To describe our hypothetical population of individuals, 
we reconstructed from aggregate statistical data a possible distribution of 
facts about individuals compatible with the reported statistics. In this 
description we used a technique, developed by Frederick Mosteller, which 
synthesizes an n-dimensional table from any number of m-dimensional 

tables that constitute its marginals. For each of the cells into which the 
population is divided, this technique determines its frequency and several 

of its media characteristics. With it we can use not only marginals but } 


also whatever estimates of interactions are available. 
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The next part of the simulation produces estimates of probable media 
exposure by the persons in the hypothetical population. An audience 
estimate is usually gotten by some kind of transformation of a circulation 
or set estimate; for example, the number of radio sets, the total daily 
newspaper circulation, the total magazine circulation, total book publi- 
cation, and soon. All circulation figures then need to be multiplied by 
factors to provide audience estimates. In such total audience estimates, 
we use the Mosteller technique to take account of all marginals about 
which we have information and all interactions about which we have 


information. 


To make our simulation dynamic, we must take into account differences 
in audience habits. For example, a newspaper read by one million people 
each day might conceivably be read by one million people 365 times a 
year or it might conceivably be read by 365 million different people only 
once. The actual distribution lies somewhere in between. We GEKerato 
this distribution of readers by frequency as the concept of cumulation. 

We have developed programs that estimate the underlying distribution of 
probabilities for any available cumulation data and generate cumulation 


curves from assumptions about the structure of the medium. 


Our simulation also accounts for duplication: that is, the incidence of 
shared audience amony different media. Exposure probabilities for each 
individual in the population are thus assigned on the basis of mean proba- 
bilities for the cell of the population to which he belongs, the likely shape 
of the probability distribution as indicated by cumulation facts, and the 
likely non-random duplication among media. In the final part of the 
simulation, a reader of a particular medium will note a particular news 
story in it with a probability p, which need not be constant. In the early 
period of a crisis, when the saliency of the events is low, p may be small; 
later, it may gradually rise. p may also vary with different groups in 
the population. With the simulation, we can thus formulate a set of 
assumptions about the flow of news ina tutal society and reach a con- 
clusion as to who has been reached and how often. At present, we have 
a working model written in FORTRAN which we are running on test data 
and we are translating the entire set of programs into AED-O in order to 


run the simulation with on-line human decision making. 
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With CRISISCOM Project support from the CONCOM Project, we have 
developed a different siumlation of an aspect of crisis commurications, 
namely the CRISISCOM siumlation (partia:ly supported by the. Naval 
Research Laboratory). Over the past year we have tested CRISISCOM 
by representing the Kaiser and the Tsar as they receive about 1500 
messages covering events which occurred during the week before World 
WarlI. On-line, an experimenter manipulates the initial perceptions of 
each decision rnaker, the parameters of his selective perception biases, 
his distortions, and his forgetting mechanisms. He can examine the 
cognitive structure and attention spaces of the decision makers and display 
the spaces in a variety of ways. Allen Kessler is the main programmer 


of this model. 
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Time-Sharing vs. Ordinary Computer Usage 
Electron Scattering Spectra 
Bound-State Wave Functions 

Nucleon-Nucleon Interactions 
Magnet Design Parameters 


Mathematical Research and CTSS 
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“Time-sharing vs. Ordinary Computer Usage - Elizabeth J. Campbell 


To compile actual statistics of CTSS compared with the usual method 
of batch processing is difficult, because the approach to the program 
is usually different. On the console, a programmer tends to be more 
adventuresome and try new approaches to a problem. The long turn- 
around time of batch processing causes a programmer to write a safe, 
conservative program. Comparison is also difficult because when a 
programmer uses batch processing, several programs are going ata 
tirne; while in CTSS, it is possible to concentrate on one program and 
save time, since a train of thought is not interrupted. The following 


examples illustrate this difference. 


In one afternoon, an inexperienced programmer and I worked out a 
routine for Q) the Legendre polynomial of the pee kind, to replace an 
inaccurate section of a program sent to us from Berkeley. We tried 
three new methods, comparing accuracy against tabulated results. 
This work would have taken us two to three days with batch-processing 


procedures. 


In another instance, I spent about three weeks trying to build a Bound- 
State program from pieces of an old program. I decided not to do it on 
time-sharing because there would be many changes due to the differences 
between FORTRAN and MADTRAN. However, when the delay proved too 
great, I] managed to convert the whole program to CTSS in a weekend. I 
decided at the time that I could have done the whole job on CTSS ina 


third to a haif the batch-processing time. 
Electron Scattering Spectra - Elaine Miller 


I have been working on a set of programs that involve problems associ- 
ated with the radiative degradation of electron scattering spectra. One 
of these programs is used as a source for predictions of data. These 
data are employed in a program that introduces all radiative effects into 


otherwise undegraded spectra. The program, when used as the source 


of predictions, computes elastic and inelastic yields throughout triangles 
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selected according to the most reliable predictions. Both straight pole 
fit and Yang assumption are available for the computation of nucleon 


form factors. 


Since, when primary electron energy is increased to 20 BeV, the yield 
often exists only within a very small portion of the selected triangle, 
this small portion is to be expanded into a larger area. Using time- 
sharing, I tested various methods of decreasing mésh size, thereby 
increasing area. Through comparison, elimination, and merging of 
methods, a new, efficient method evolved. For each set of data, time- 
sharing will be used to determine a new triangle, that strictly adheres 
to the conditions imposed upon the original triangle, but uses a finer 


mesh to obtain additional information. 
Bound-State Wave Functions - Elmer Bartels 


Ihave revised and perfected a program, which Elizabeth Campbell 
started with the Bound-State calculation, as the first of a package of 
programs for the use of physicists. A question-and-answer set of 
routines compute the normalized Bound-State wave function and binding 
energy for a Wood-Saxon well. The user may LOADGO the proper set 
of files and then type in well parameters and eigenvalue quantun. numbers 
as they are requested. He has the option to guess the binding energy, or 
the program will supply the guess and subsequently compute the wave 
function and eigenvalue. Another option in the program stores from two 
to four Bound-State wave functions for computing the radial integral 
2 
Le ( R)R,RzRyAR. 
° 
The user may either compute another set of wave functions and binding 
energies with the same well parameters, or he may return to the begin- 
ning of the program and insert all new parameters of the Wood-Saxon 


well without having to LOADGO the files again. 
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Nucleon-Nucleon Interactions - Nancy Spencer 


This use of CTSS is to investigate the applicability of the boundary condition 
model, developed by Professor Feshback and Dr. Lomon, to the proton- 
proton and neutron-proton scattering data in the energy range extending 

to 350 MeV. A basic set of parameters, necessary to specify the 

problem, are read into a program, which calculates such data as scatter- 


ing length, partial cross section, total cross section, and polarization. 


I then compare these theoretical results, predicted by the model, to 
experimental data at various energies and apply a chi-square criterion 

to determine goodness-of-fit. After the program is run off-line, each 
parameter is varied systematically by a specified amount until an optimum 


set of parameters appears. 


Since only one parameter varies at a time, this process does not take 
into account any correlation which may exist between the parameters. 
Thus, with CTSS, a physicist can examine the results as they are com- 
puted and can interrupt the program and change a parameter whenever 
he sees a correlation or feels that one may exist. With such an inter- 
action between a physicist and a computer, the development model, 


which gives a precision fit to the data, can proceed rapidly. 
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Magnet Design Parameters - James E. Spencer 


The main purpose of this problem is to produce a computer program 
capable of calculating all first.and second-order focusing coefficients 
for a serial array of magnets. The design parameters for these 
magnets are arbitrary, so that any geometrical form of boundary can 
be studied. Functional forms of greater than second order (for the 
geometrical form of boundaries) will not affect the first and second- 
order focusing coefficients, but they are included so the program can 
be used as a general ray-tracing subroutine. A realistic field descrip- 
tion is being used, including fringing fields at entrance and exit bound- 
aries which take into consideration the effects of finite widths and 


curvatures on the iso-induction lines. 


Analytic expressions for the fringing fields in the median plane are 
obtained from experimental data by curve fitting. A knowledge of the 
field off the median plane is then obtained using Taylor expansions 
through the third order. The resulting expressions for the field are 
then used with equations of motion to trace rays of any initial conditions 


through the system. 


Once the performance of this program has been determined by making 
specific calculations on known cases, other cases of interest will be 
investigated with a view towards learning the utility of these coefficients 


in the actual design of magnets. 
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Mathematical Research and CTSS - Louis N. Howard 

Fi 
This project has been largely directed toward exploring some of the 
potentialities of CTSS for assisting a mathematician in "everyday"! 
aspects of his research. In general, most pure mathematicians, and 
many applied mathematicians, almost never use a computer. This is 
probably because computers have seemed most appropriate for large 
computations on completely-formulated and numerically-analyzed prob- 
lems; and mathematicians are seldom interested in such problems. 
However, it seems that there are various special aspects of mathemati- 
cal research in which the ability to rapidly make relatively-simple, 
tedious computations, not necessarily of a strictly numerical nature, 
could be of considerable use; for example in seeking counter-examples, 
jesting conjectures, or trying to obtain an overall numerical or geomet- 
rical understanding of some formula or relationship. CTSS, used with 
certain basic programs designed to be fairly self-explanatory, seems 
to offer the possibility of providing such service to mathematicians 


without requiring much knowledge of programming. 


We have experimented with several programs oriented toward this kind 

of work. One example is a very flexible and almost self-explanatory 

program for constructing rough typewriter graphs of contour lines for 

arbitrarily-given real functions of two-variables. This is quite easy 

for the uninitiated to use, and has been found valuable in obtaining a 

qualitative ''feeling'' for the behaviour of such functions. Another 

example is a program for making algebraic computations with matrices 

whose entries are elements of an arbitrary finite field. The field is 

specified by giving its characteristic p, number of elements, and an 

irreducible polynomial over the integers modulo p. A number of matrices 

can be handled simultaneously and operations with them are carried out 

by giving various simple commands. A third example is a program for 

constructing the Lagrange interpolation polynomial through an arbitrary 

set of points and carrying out various operations with this polynomial on 

the basis of simple commands. Other similar programs have been { 
explored for computations in group theory, differential equations (specially 

adapted to time-sharing), and number theory. Finally, a program still | 


under development is intended to provide a simple and flexible means for ~ 
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carrying out operations and sequences of operations that occur in 
elementary numerical analysis. Ifa suitable program can be developed, 


we hope to have students use it fora numerical analysis class. 


While the greatest use of CTSS by mathematicians is likely to be with 
the type of small-scale everyday problems mentioned previously, some 
mathematicians will occasionally be interested in relatively-large-scale 
computational problems. Two such problems in which time-sharing 
appears to be particularly useful have been explored, partially at 
Project MAC and partially at the Computation Center. One is an 
elaborate program for studying the structure of finite groups, developed 
by Dr. C. C. Sims; the other, a set of programs for the study of prob- 
lems in hydrodynamic stability, developed by L. N. Howard and Prof. © 
M. Landahl of the Aeronautics Department. In both of these cases, the 
possibility of guiding computation on the basis of previous results has 


proved to be very useful. 


In general, these experiments suggest that a time-sharing system can 
be effectively and usefully employed in some parts of mathematical 
research; though for general applicability, further simplification of 
programming, and simplifying at least a part of the command system 
of CTSS, should be brought about. The principal limitation of the 
specialized programs mentioned previously is that, in the case of 
complicated, flexible programs which actually do relatively little com- 
putation, like the contour-line program and the numerical-analysis 


program, the amount of time lost in swapping seems excessive. 
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On-Line Simulation System 
The OPS-3 System 
Computer Techniques in Business Problems 
Console-Operated Statistical Routines 
An Automated Stock Exchange 
Computer-Planning System 
Industrial Dynamics 


On-Line Array Processor 
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On-Line Simulation System - Martin Greenberger and Malcolm M. Jones 


The availability of time-sharing systems makes development of an on-line 
simulator both feasible and desirable. From our experience with GPSS, 
SIMSCRIPT, and a preliminary; on-line simulator called OPSIM, we have 
drawn up specifications for a flexible, interactive simulation system. Our 
on-line simulator uses OPS, a’ multi-purpose system with compound 
operators (KOP). One KOP may call another and transmit parameters 

by way of a calling sequence; one may also activate or schedule the occur- 
rence of another by means of special schedule operators. A special com- 
pound operator, called AGENDA, contains a list of future activities, which 
are automatically removed from AGENDA aiter it has been activated. A | 


user may at any time inspect and modify AGENDA. 


The execution of a compound operator may result in the alteration of state 
variables stored in the common system, or it may result in changes to 
AGENDA, Among the operators that change AGENDA are WAIT and DELAY 
WAIT is conditional upon the value of a state variable, and DELAY is un- 
conditional and specifies that a certain period of time elapse before execu- 


tion of the KOP is continued. Thus, activities listed in the AGENDA have 


either a time of execution or a Boolean condition associated with them. f 
Two other operators allow activities scheduled in the AGENDA to be can- i 
celled or rescheduled. In addition, the LOCAL operator lists variables 
to be saved when an activity is interrupted by DELAY or WAIT. 
Among the simulations in progress or completed are: 
1. A simple multi-server queuing model , 
2. Aninventory distribution system \ 
3. A process-control computer system / 
4. A three-level priority interrupt scheduling system 
(similar to the CTSS scheduling algorithm) 
5, A hospital reception desk 
6. An automated stock exchange. 
} 
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The OPS-3 System - Martin Greenberger, James H. Morris, Jr., and 
David N. Ness 


OPS-3 is a multi-purpose on-line system within which a user can write 
and test his own programs regardless of the nature of his particular pro- 
ject. Since the system is comnletely open, a usexv can add his own sub- 
routines and use them in conjunction with system subroutines, as long as 
they obey customary conventions regarding common storage. He is free 
to run his programs within the system during the formative and testing 


stages, and then run them independently during the production stage. 


The system contains operators for symbolic vector and matrix calcula- 
tions, statistical analysis, live data-base reorganization, dynamic loading 
of routines, FORTRAN-like program writing with instant execution, and 
on-line simulations. While sitting at a console, the user may access 
information about the system. The information is stored in regular CTSS 
disc files and may be retrieved from within the OPS-3 system by a single 
key-word mechanism. OPS-3 is now being used in a number of projects. 
These include: ent 

1. an on-line simulation system within which a user may change 
the course of the simulation or examine partial results, 

2. a real-time simulation of the stock exchange that involves both 
simulated and console-initiated orders, 

3. a polynomial manipulation system, 

4. a system for statistical analysis in the social sciences, 

5. a general debugging system for MAD programs in which 
programs are first run interpretively and then compiled into regular BSS 
subroutines. 


OPS-3 is currently available as the CTEST6 command at Project MAC. 
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Computer Techniques in Business Problems - Donald C. Carroll, 


Paul Clermont, Thomas J. Johnson, Lalit Kanodia, and Francis J. Russo 


This project investigates the use of heuristic programming and interactive 
problem-solving techniques in complex business problems, such as 
scheduling and sequencing production in a job shop. We developed a 
simulation model of a hypothetical shop, and, by testing some heuristics 
procedures for the sequencing and routing of jobs under varied and ex- 
tensive conditions, determined the general steady-state behavior of this 
model. The first part of our research studied fully-automatic decision- 
making, in which the heuristics for sequencing and routing require no 
human guidance; the second part is to study the rule of a human pattern- 
recognizer and decision-maker who is closely coupled with the automatic 
system. In the latter research, we are investigating the establishment 

of manning levels for the shop, overtime decisions, and load-leveling. 
These decisions determine the structure under which an automatic decision 


system must operate; consequently, we call them structuring decisions. 


A coinplete set of programs that study non-interactive decision making is 
now available. It includes a highly-parametric input generator (MADTRAN), 
a variety of simulation models that encompass such features as multiple- 
channel work stations, multiple-component orders, alternate-routing 
possibilities, and machine substitution (FAP), and finally an output analysis 
program (MADTRAN). Three pieces of research are complete. One, a 
study of heuristic sequencing rules for single and multiple components 
orders, demonstrated the general superiority of the so-called COVERT 
rule for sequencing. Another, which dealt with heuristics for selecting 
alternate routes for jobs, produced a superior heuristic derived from 
COVERT. The last study, on the effect of feed-back phenomena in heavily- 
loaded shops, led to several generalizations about the control of inputs 

and processing rates to obtain stable operation. A by-product of these 
studies of the minutiae of the decision system is a comprehensive under- 
standing of the behavior of these models with respect to their aggregate 


behavior. 


4 


An extension of these studies to heuristics for machine substitution is 


almost complete; however, the major effort is in interactive 
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decision-making. Toward this end, we are elaborating some of the basic 
programs to include such structural parameters as manning levels by 

shift, and we are developing generalized information retrieval and display 
programs for the ESL Display Console. A rudimentary version of our inter- 


active system should be available in the fall. 


Console-Operated Statistical Routines - James R. Miller and 
Stephen Whitelaw 


A battery of statistical routines have been written, tested, and debugged 
for use with a CTSS console, primarily for the testing of hunches and 
limited hypotheses. Data are typed directly into the console, analyses 
are performed, and results are printed out immediately. Many of our 
routines are non-parametric in nature, and are oriented toward the 


analysis of typical problems in the social sciences. They include: 


1, Standard two-dimensional contingency analyses. 

2. Special three-dimensional contingency analyses. 

3, Rank-order intercorrelation and partial correlation analyses. 

4. Binomial tests. 

5. Single and multiple percentage difference tests. 

6. Homogeneity test on samples classified strictly according to 
nomial categories. 


7. The Mann-Whitney U test of the difference of the sample medians. 


We have also written several parametric routines: 


1. Linear intexweurelation and partial correlation analyses. 
2. Simple and multiple regression analyses. 

3. A T-test of the difference between two sample means. 

4, A one-way analysis of variance. 


5. Atest for normality of sample data. 


We shall expand the range and capability of our current library by improv- 
ing existing routines and by adding new ones. When S. Whitelaw has 
adapted the routines to the OPS system, a user will be able to call routines 
by name from the console, with the name of the data matrix as argument, 


to compound routines, and to save data and results for future use. 
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An Automated Stock Exchange - Stephen Whitelaw 


With the OPS system, methods of automating large securities markets, 


such as the New York Stock Exchange have been explored. This research 


-has included the programming of subroutines or operators to automatically 


handle all functions of a large stock exchange: entry of orders, consumma- 
tion of orders, modification of completed transactions, and automatic 
bookkeeping. Additional operators embody the decision rules of specialists 


who maintain a market by committing capital. 


Several sets of decision rules have been tested on different simulated 
markets. The performance of these rules is measured by price and net 
return criteria. The market model accepts simulated orders from an 
order generator. By treating a common file as an intermediate storage 
area for incoming orders, reports of completed transactions, requests 
for quotations, and output of information, it can also accept orders 
entered from other consoles. The system gives faster response times 
than those found in non-automated markets, and a programmed specialist 


is able to control the market under very adverse conditions. 


Currently the system is being reprogrammed to function in two modes: 


mode 1 provides for continuous simulation and output of consummated 


transactions, and mode 2 for the entry of console orders into the simulation. 


Since the console remains undisturbed by this entry, it can function both 
as a ticker tapéthat reports the results of the simulated market, and also 
as a broker's device that enters orders into the same market. This 
functioning should make an éffective management game and extend the 

use of the automated exchange as a research tool. The simulated market 
will be synchronized with the computer clock to prevent the simulation 


from getting ahead of real time. 
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Computer-Planning System - James C. Emery 


This project studies a management computer planning system. Because 
planning requires a system that can be augmented and modified extensively, 
careful attention has been paid to general techniques that provide for such 
evolution. With such a model, a planner will be able to make sequential, 
hierarchical search for improved plans. First, from an aggregate model 
he will construct a large-scale plan, which is generated through a sequential 
search of alt rnatives. In this search, the human planner proposes alter- 
natives and the computer transforms a proposed plan into its predicted 
consequences. Each probe for an improved plan provides the planner with 
additional information for the next probe. If the process converges, the 
aggregate plan emerging from it is then amplified by more detailed models 
within the constraints provided by large-scale plans. This hierarchical 
process continues until plans are sufficiently detailed for execution. 

os 
The system, specifically designed for budget planning, accepts a fairly 
wide class of functional relationship among variables. It will also deal 
with hierarchically-defined variables, such as manufacturing cost defined 
as the sum of direct labor, direct material, and overhead cost. Since 
functional relationships and the hierarchical definition of variables are 
represented primarily in the form of list structures, the system has the | 
advantages of both flexibility and economy of storage. During the past 
year, Ihave developed the general structure of the system and have also 
begun to specify storage organization and the algorithms for the retrieval 
and manipulation of stored data and to define the set of subroutines needed 


for the generation of displays for budget planning. 
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Industrial Dynamics - Jay W. Forrester 


The Industrial Dynamics Group uses time-sharing and the DYNAMO Com- 
piler for analysis and design of feedback processes in industrial organiza- 
tions. These systems contain multiple-loop feedback structures with 
highly-nonlinear relationships. Sixtieth-order models of systems have 
been simulated which contain several hundred variables. Such mathema- 
tically intractable systems contain important modes of behavior which 
are totally absent in the simpler linear systems. The DYNAMO Compiler 
is very efficient; for, in only a few seconds of computer time, it can 
simulate a system of a hundred or more equations through two hunared 
time steps. The Compiler contains a graphical plotting routine which 


makes time plots of systems responses on remote printing stations. 


One use of the system has explored the dynamics of corporate growth. 
Many of the important growth processes depend on the interactions beiween 
positive feedback loops, negative feedback loops, and the coupling between 


these created by nonlinear elements. Another use has developed illustra- 


tive material for a student seminar on the principles underlying the behavior 


of feedback systems. These principles are generalizations about system 
behavior which have not previously been organized and illustrated. Some 

of them are implicit in the mathematics of linear feedback systems, but 
have not been clearly identified and illustrated; others exist only in muliiple- 
loop nonlinear systems and could not have been previously identified with- 


out the extensive empirical results obtained by computer simulation. 
On-line Array Processor - Mayer E. Wantman 


Operators have been developed on the OPS-3 on-line system that provide 
a simple, flexible way of manipulating and analyzing numerical data within 
the computer. They can be used for general processing arrays, fitting 

of regression models, transformation of variables, rank-ordering of data, 
calculation of roots of polynomials of arbitrary degree, error analysis, 


generalized desk calculation, and MAD-type programming. They are: 


1, Assignment Operator. The operator SET specifies and executes 
MAD-type assignment statements. Operands may be vectors and matrices, 


as well as elements. In addition, SET possesses embedded functions, 
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calls, matrix multiplication and transposition, a power generator, anda 


variety of logical operations. Lines are of the form 


SET ALPHA=(BETA + GAMMA)*DELTA. P. (A + B) 


2. Regression Operator. The Operator FIT creates a multiple 


linear or simple polynomial model derived from given arguments. Sample 
points may be weighted and residuals saved. As an example, multiple 
linear regression of a dependent variable (stored in the vector named 
ALPHA) on independent variables (stored in the vectors named BETA, 
GAMMA, and DELTA) would be performed by the line 


FIT ALPHA TO BETA GAMMA DELTA. 
A cubic fit of ALPHA to BETA would be specified by 


FIT ALPHA TO BETA DEGREE 3. 


Be Root-Solving Operator. The operator ROOT uses the Newton- 


Raphson method of successive approximations for the determination of the 
root of a polynomial equation in one variable. Since the user Specifies 
the starting value, he can obtain all roots by setting different initial values. 


Exponents need not be integral. 


4. Ordering Operator. The operator RANK takes a vector of 
numbers and creates a second vector which is an ordered version of the 


first. A vector which describes the mapping performed may also be 


calculated. 


5. Sorting Operator, The operator SORT sorts by intervals ona 
specified vector. The intervals may be specified either explicitly or by 
range and spacing. If any samples fall outside the sorting range, an 
appropriate message is printed. 
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APPENDIX A 


PROJECT MAC MEMORANDA 


SUBJECT AUTHOR DATE 


On-Line Documentation of the Ite 
Compatible Time-Sharing System 


SQZ/PADBIN (Ref: MAC-M-157)_ R. 
Abstracts of CTSS Console Com- E. 


mands (CC239) (174-2, 1/4/65; 
174-3, 3/7/65 


Access to the MAC System via R. 


Alternate-Use- TWX (''TWX? Prime") 


String Manipulation in the New D. 
Language (AIP memo 71) 


CTSS Library Abstracts (CC240) Bs 
(177-2, 11/23/64; 177-3, 3/22/65 


Winett 7/21/64 


Hussey 7/20/64 


Kliman 8/3/64 


G. Mills 7/24/64 


Bobrow 7/28/64 


Kliman 8/6/64 


Crunch Command in CTSS R, Bayles 8/7/64 
BEFAP Command within CTSS R. Bayles 8/10/ 64 
The ARCHIV Command A. Scherr 10/19/64 
The Structure of On-Line Infor- E, Glaser 10/1/64 
mation Processing Systems J. Dennis 
Preliminary Notes on the Hard E, Glaser 8/28/64 


Matter for the MAC 635 


Proposed Instructions on the M. 


GE 635 for List Processing and 
Push Down Stacks (AIP memo 72) 


On Memory Addressing in the M. 
M. Levin 


GE 635 
MAC and Its Users 


U 
A Note on Asynchronous Parallel H. 
Processing (MSG memo 5) 


Nesting and Recursion of Proce- J 
dure in Segmented Memory E, 
(MSG memo 6) 


Levin 9/18/64 


Minsky 9/22/64 


Neisser 9/29/64 


Witsenhausen 
7/64 


. Dennis 10/ 28/ 64 


C. Van Horn 


Automatic Scheduling of Priority J. Dennis 10/ 64 
Processes 
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MEMORANDUM 
MAC-M-No, SUBJECT AUTHOR DATE 
= Ss PUBJECT pi Re CESS 
189 An Example of Intersphere Com- J. Dennis 9/25/64 
munication and Asynchronous 
Parallel Proces sing: Typewriter 
Console Message Handling by 
Protected Service Routines 
(MSG memo 8) 
190 General Comments about PDP-6 D. Edwards 9/25/64 
Coding and Hardware (PDP-6 
memo 1) 
191 TECO 6 (PDP-6 memo 2) D. Edwards 
192 Unrecognizable Sets of Numbers M. Minsky 11/2/64 
(AIP memo 73) S. Papert 
193-2 TYPSET and RUNOFF, Memor- J. Saltzer 11/6/64 
andum Editor and Type-out 
Commands (CC 244-2) 
194 On the Control Specifications of D. Kuck 11/64 
Computation 
195 ED, A Context Editor for Card R. Dailey 11/20/64 
Image Files (195- 1, 3/15/65) 
(CC245; CC234-1) 
196 Some Aspects of the Structure of D. Kuck 11/13/64 
' Computation R, Kain 
197 Some Time-Sharing Ideas D. Kuck 11/13/64 
198 PLEX-Dump and Relocation in T. Fox 10/26/64 
the AED-O System (9442-M-] 10) 
199 Stack Manipulation Facilities of R. Coe 11/18/64 
the AED-O Language (9442-M-112) 
200 Processes, Spheres of Protection E, Van Horn 
and Independent Computation (no date) 
(MSG 9) 
201 New Operating System for the R. Bayles 11/23/64 
ESL Display Console (9443-M-118 
(Cont'd 12/10/ 64) 
202 Proposal to Improve the Rotation R. Stotz 11/27/64 
Matrix of the ESL Display Console 
(9442-M-114) 
203 A Compacting Free Storage Pro- B. Wolman 3/5/65 
gram (9442-M-113 
204 AED-O Compiler for Batch Pro- H. Spencer 12/2/64 


cessing (9442-M-115 
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MAC-M-No, SUBJECT AUTHOR DATE 
Oe SUBJECT SL ROME RCASG eee 
205 MADBUG: A MAD Debugging R. Fabry — 11/25/64 
System (CC247) 
206 CTSS LISP Notice: Supplement to T. Hart 12/7/64 
A. I. Memo No. 67 (AIP memo 74) 
207 Internal Memos for AED Users C. Feldmann 12/8/64 
(9442-M-116) 
208 ' AED Flash 10, New CTEST2 C, Feldmann 12/10/64 
Command (9442-M-119) 
209 The Naval Mine Warfare Game I. Sutherland 12/18/64 
210 Graphical Output Device A. Rutchka 12/21/64 
211 New TWX-Prime Lines R.G. Mills 12/23/64 
212% Translation of Communication O, Selfridge 6/10/65 
with Automata by Carl Adam Petri 
ee mate. 
213 AED FLASH 11: AED BUG Usage T. Fox 12/28/64 
(9442-M-121) 
215 Television Camera-to-Computer M. Minsky 1/9/65 
Adapter: PDP-6 Device 770 
(AIP memo 75) 
216 New B-Core System for Program- C, Lang 4/30/65 
ming the ESL Display Console 
(9442-M- 122) 
217 Operating Manual for the ESL R. Stotz 3/9/65 
Display Console (9442-M-129) J. Ward 
218 Social Science Data Bank I. D. Poole 1/21/65 
Development 
219 The COMIT Feature in LISP II D. Bobrow 2/18/65 
(AIP memo 76) (LISP II Project 
Memo 2) 
220 The COGO Language D. Roos 2/25/65 
221 The INFO Comimand J. Winett 2/25/65 
222 AED FLASH 12: Multi- Entry D. Ross 2/25+65 
Procedures (9442-M- 124) 
224 An Automatic Plotting Subroutine, T. Stockham, 3/3/65 
“te XYPLOT Jr, 


* In older listing's of MAC-M-memos, this merno was numbered as MAC-M-214, 
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235 


236 
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238 


239 


240 


241 


SUBJECT 


AED FLASH 13: Argument Check- 
ing Procedures for AED (9442-M- 
125) (225-2, 4/15/65) 


AED FLASH 14: Availability of 
AED JR Systems (9442-M- 126) 


The Editing Routine 


A General Cross- Tabulation 
System (COM/COM Sim. memo 23) 
(See Additional Information on 
Cross- Tabs: New Output) 


Proposed Datanet- 30 Program 


Mind, Matter and Models 
(AIP memo 77) 


Description of SNOBOL Commands 
(CC235-2) 


Preliminary Description of 
DAEMON, the Disk Dump and Re- 
load Program Package, and its 
Operation (CC 252) 


Trial Launching of the OPS-3 
System 


Users' Manual on the Use of 
Structural Language 


BLODI Command 

AED FLASH 15: Execution Time 
Procedure Calls (9442-M-131) 

A Generalized File Structure and 


Input/Output System (CC 241-1) 


Some Facilities of the MAC 
Document Room 


KREADQ, SETRDQ, MODRDQ, 
A Set of Format-Free Input 
Subroutines 


Topics in Model Theory 


PDP-6 LISP Input-Output for the 
Dataphone (AIP mezno 79) 
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AUTHOR DATE 


J. Walsh 3/5/65 


D. Ross 3/9/65 


D. Griffel 3/8/65 


W. Selles 4/24/65 


S. Dunten 3/18/65 


M. Minsky 3/26/65 
L. Pouzin 4/9/65 


M. Bailey 4/26/65 


M, Jones 4/26/65 
D. Ness 


R, Logcher 9/15/65 
O. Wright, 5/4/65 
Jr. 

J. Walsh 4/22/65 
%. C. Daley 5/7/65 
R, J. Creasy 

R, M. Graham 
Document 5/10/65 
Room 


(A. Bowen) 


D. Ness 6/16/65 


M. Levin 5/25/65 


W. Martin 6/1/65 
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242 


243 


244 


245 


246 


247 


246 


249 


250 


251 


252 


253 


254 


255 


256 


Zor 


SUBJECT 


PDP-6 LISP Input-Output for the 
Display (AlP memo 80) 


Specifications for a Dataphone- 
Driven Remote Display Console 
for Project MAC 


Proposal tor Commands for 
Dynamic Resource Allocation 
(CC 254) 


DIRECT; A Command Subroutine 
for Controlling Programs 


Computer Experiments in Finite 
Algebra 


AED FLASH 16: Use of AED JR to 
Parse FORTRAN Il (9442-M-134) 


Use of MACDMP (AIP memo 83) 


MAC PDP-6 DECtape rile 
Structure (AIP memo 82) 


PDP-6 TECO, July 1965 Version 


AED FLASH 17: Features of the 
Revised AED-JR (9442-M- 136) 


AED FLASH 18: A Subroutine for 
Generating Coons! Surfaces 
(9442-M- 137) 


FMS Tape Distribution 
(9442-M-127-3) 


Derivation of a Mean Pseudo- 
Associative Search Time 
(MSG memo 13) 


Surfaces for Computer-Aided 
Design of Space Figures 
(9442-M-139) 


Three- Dimension Pseudo Pen 
Sub-routine for use with the ESL 
Display Console (9442-M-138) 


Syntax and Display of Mathemati- 
cal Expressions (AIP memo 85) 


AUTHOR DATE 


W. Martin 6/28/65 


R. H. Stotz 6/24/65 
J. E. Ward 
U. F. Gronemann 


G. Schroeder 4/17/65 


T. Crystal 6/23/65 
E. River 


W. D. Maurer 
6/14/65 


H. Spencer €/ 25/65 


P, Samson = 7/1/65 


P, Samson 7/1/65 


P. Samson = 7/23/65 


J. F. Walsh 7/1/65 
R, Feldmann 


C, Lang 6/30/65 


R, Feldmann7/ 1/65 
D. E, Walker 


E. C. Van Horn 
8/16/65 


S. A. Coons 7/21/65 


R, Polansky 
7/13/65 


W. Martin 7/29/65 
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A Subroutine Trace Program 
forCTSS ; 
CTSS Console Commands 
CTSS Library Subroutines 
TYPSET and RUNOFF, Memo- 
randum Editor and Type-out 


Commands (CC-244-2) 


ED, a Context Editor for FAP 
and MADTRN Program 


AUTHOR 


B. Wolman 


E. Kliman 
E, Kliman 


J. Saltzer 


R. Dailey 
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DATE — 


3/5/65 


3/7/65 
3/22/65 
1/11/65 
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Amstutz, A, E., A Management Oriented Behavo::ul Theory of Inter- 


actions within Consumer Product Markets, M.I. T. Sloan School of 
ee ae Ln ee Ae CLA TSE LS 
Management, Ph.D. Thesis, June 1965 


Armstrong, A, E., A Braille Telecommunication Terminal, M.I. T. 


Department of Mechanical Engineering, S.M. Thesis, June 1965 


Beckreck, L. N,, SEPOL: A Soil Engineering Problem-Oriented 


Language, Department of Civil Engineering, M.S. Thesis, June 1965 


Beusch, J. U., Dynamic Behavior and Control of Communication Networks, 
a! ee LO OV OMMU NLC A LLONEINELWORKS 
M.I.T., Department of Electrical Engineering, Ph.D, Thesis, 
May 1965 


Brach, J. R., A Problem-Oriented Language for Project Management, 


M.1I.T., Department of Civil Engineerinz, S.M. Thesis, January 1965 


Carroll, D, C,, Heuristic Sequencing of Single and Multiple Compound 
Jobs, M.I.T., Sloan School of Management, Ph.D. Thesis, June 1965 


Chang, W. Y., Instability Analysis of Transverse Waves in Solid-State 
Plasmas, M.I.T., Department of Electrical Engineering, B.S, Thesis, 
June 1965 


Chiodi, W., A Mechanism for the Transformation of S, M.I. T, ; Depart- 
ment of Electrical Engineering, S.M. Thesis, June 1964 


Diamond, D. S., A Step toward Coniputer Aided Preparation of Magazine 
Advertising Copy, M.I. T, Sloan School of Management, B.S, Thesis, 


June 1965 


Dvorak, A. Jr., An Input-Output Program for Electronic Circuits Using 
a CRT, M.I.T., Department of Electrical Engineering, Thesis, 


June 1965 
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Efimba, R, E,, Analysis of a Bridge Deck as a Plane Cie View! sae 
——— EB CK Sa Slane Gnid, 
Department of Civil Engineering, M.S, Thesis, June 1965 


Ellis, H, E,, A MAD-based Programming Language, M., I. T. » Depart- 
ee a rogramming Language 
ment of Electrigal Engineering, B.S, Thesis, June 1965 


Falk, G., An Information Retrieval System for Reference Filing, M.1I.T., . : 
ee pine 


Department of Electrical Engineering, B.S, Thesis, June 1965 


Finkelstein, S, C., An Analysis of a Space Beam, M.1.T., Department 


of Civil Engineering, M.S. Thesis, June 1964 


Fiske, R. S., A Real-Time System for Computer Monitoring of the Stock 
emptor Monitoring of the Stock | 


Exchange, M.1, T., Sloan School of Management, B.S. Thesis, 
June 1965 


POR hy Fe 5 Special Input-Output Facilities for the IBM 1620, 


M,.1.T., Department of Civil Engineering, M.S, Thesis, May 1964 
Goeke, R. F,, Introduction to Machine Language Programming ona Time- 
pal be eee LS Gey 
Shared IBM 7094, IMG, 1 Bla Department of Electrical Engineering, 
B.S, Thesis, June 1965 


Gonzalez-Zubieta, R., On Some Aspects of Integer Linear Programmin x 


Ph.D, Thesis, Sloan School of Management, June 1965 (Also 
released by Operations Research Center as Technical Report 
No. 16, June 1965) 


Goodman, R, V., A Preliminary Method of Analysis for a Computer-Aided 


Structural Design Language, Wika hy ALS Department of Civil Engineering, 


S.M. Thesis, June 1965 


Groneman, U. F., Coding Color Pictures, M.I.T., Department of Elec- 


trical Engineering, Ph.D. Thesis, February 1964 


Guttmann, E, G., Investigation of Threshold Logic in High-Speed Displa 
8 es ord =ogic in hig BAAS AEE) 
Systems, M.I.T., Department of Electrical Engineering, S.M. Thesis, 


June 1965 J 7 
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Guzman Barron-Torres, G., Optimum Design of Prestressed Concrete Highway 
Bridges using Digital Computers, M.I.T., Department of Civil 


Engineering, S.M. Thesis, May 1965 


Hadley, R. A., Investigation of Automated Direct Design of Rigid Steel 
Planar Frames, M.I. T., Department of Civil Engineering, M. Ss 


Thesis, February 1965 


Haines, E. C. Jr., The TREET List- Processing Language, M.I.T., 
Department of Electrical Engineering, M.S, Thesis, September 1964 


(Also produced by Mitre Corporation as Report SR-133, April 1965) 


Hamann, W. C., Computer Aided Design of Slender Structural Members, 
M.I.T., Department of Mechanical Engineering, M.S. Thesis, 


May 1964 


Hamilton, M. L., and A. D. Weiss, An Approach to Computer-Aided 
Preliminary Design, M.I. T., Department of Mechanical ee eerie 
and Naval Architecture, S.M. Thesis, August 1964 (Also as 
ESL-TM-228, January 1965) 


Hershdorfer, A. M., Optimal Routing of Urban Traffic, M.1.1., Depart- 
ment of Civil Engineering, Ph.D, Thesis, June 1965 


Hoag, J. D., Instabilities in Transverse Waves Along 6 for a Beam Type 


Distribution, M.I.T., Department of Electrical Engineering, 
B.S. ‘Thesis, October 1964 
& 
, 
F Ishii, H., Instabilities in InSb, M. 1. T., Department of Electrical Engineer- 
\ ing, S.M, Thesis, May 1965 


Johnson, W. F., Use of Traffic Volume Data in Evaluation of Highway 
ae User Costs, M.I. T., Department of Civil Engineering, S.M, Thesis, 


June 1964 


{ Kanodia, L. S., Queue Balancing - Another Approach to the Theory of 
: Scheduling, M. LT. Sloan School of Manag>ment, M.S, Thesis, 
June 1965 
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Kanstroom, E,, An Information Retrieval Program Compatible with the 
M.I. T., Time-Shared Computer System, M.I. T. Sloan School of eens 


Management, B.S. Thesis, June 1965 . 


Kawaguchi, S, K., Preliminary Computation and Data Reduction for STRESS, 


M.I.T,, Department of Civil Engineering, S.M. Thesis, August 1964 


Keller, J., The Critical Path Method, Its Implementation and Effective 
a ee ae ee eee GLE TAL LO DEAN GRILLE CLIME 
Utilization in Construction, M.I.T., Department of Civil Engineer- 
ing, M.S. Thesis, June 1964 


Kim, Numerical Analysis of Open-Channel Transients from Hydro; Power 
Operation, M,1.T., Department of Civil Engineering, S.M. Thesis, 
February 1965 


Kusse, B., Plasma Dispersion Relations and the Stability Criteria, M.I, T, 
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